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ABSTRACT 


This  report  suanarizes  the  work  performed  during  the  last  three  years  on 
an  ARC  sponsored  contract  with  NCSU  for  the  purpose  of  investigating  the 
potential  of  compound  semiconductor  and  hetero junction  devices  for 
millimeter-wave  sources.  Both  two-terminal  IHPATT  and  three-terminal  MESFET 
structures  were  considered.  A  series  of  theoretid  device  models  were 
developed  and  used  to  investigate  device  performance.  Non  equilibrium  effects 
were  considered  with  the  use  of  a  theoretical  Monte  Carlo  transport  model. 
Prototype  IMPACT  devices  were  fabricated  from  GalnAa/InP  grown  by  LFE  and  GaAs 
layers  grown  by  MBE  and  OM-CVS.  The  results  of  this  program  indicate  that 
compound  semiconductors  should  be  useful  for  millimeter-wave  sources.  In 
particular,  both  GaAs  and  InP  IMPACTS  should  be  capable  of  conversion 
efficiency  exceeding  10%  at  100  GBz.  Three-terminal  MESFET  devices  fabricated 
from  InP,  GalnAs  and  GalnAsP  should  be  capable  of  producing  useful  gain 
through  at  least  60-70  (ax.  With  sub-half  micron  gate-lengths,  the  upper 
frequency  of  these  devices  could  possibly  reach  100  GHz . 
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I.  SUMMARY  OF  WORK 

The  purpose  of  this  reseerch  program  ns  to  investigate  the  use  of 
con pound  and  hetero junction  seaiconductor  devices  for  use  as  millimeter-wave 
sources.  Both  two-terminal  IMPACT  devices  and  three- terminal  FET  type  devices 
were  considered.  The  prograa  consisted  of  both  theoretical  and  experimental 
efforts. 

A  series  of  comprehensive  and  accurate  theoretical  device  models  were 
developed  and  used  to  investigate  the  high-frequency  potential  of  devices 
fabricated  from  various  compound  semiconductors.  IMPACT  devices  fabricated 
from  GalnAs/InP,  InP,  GaAs  and  Si  were  investigated  and  predicted  performance 
compared.  Device  performance  was  investigated  from  the  low  microwave  region 
to  over  100  GHz.  In  general,  our  studies  indicate  that  GaAs  is  the  optimum 
material  for  maximizing  device  efficiency  to  a  frequency  of  about  40  GHz,  InP 
works  best  from  about  40  GHz  to  60  GHz  and  Si  produces  the  greatest  efficiency 
for  operation  above  100  GHz.  The  use  of  hetero junction  devices  such  as 
GalnAs/lnP  shows  considerable  potential  for  device  development  and  should  be 
further  investigated.  In  the  structure  we  investigated  the  GalnAs  layer  was 
used  to  reduce  the  contact  resistance  to  the  p-type  InP.  Such  a  structure 
should  allow  the  attractive  properties  of  InP  for  high-frequency  operation  to 
be  exploited.  Our  studies  indicate  that  conversion  efficiency  greater  than 
10Z  should  be  obtainable  at  100  GHz  from  devices  fabricated  from  both  GaAs  and 
InP.  Doping  profile  control,  however,  is  critical  and  growth  technology  such 
as  MBS  and/or  OM-CVD  will  be  required. 

Our  large-signal,  numerical  IMPACT  device  simulations  are  the  first  to 
predict  from  basic  transport  phenomenon  device  conversion  efficiency  for  GaAs 
Read-prof i-le  devices  greater  than  20*  at  Ka-band.  Raytheon  has  reported  such 
conversion  efficiency  for  experimental  Read-profile  devices  and  previous 


L 


2 


theoretical  calculations  have  not  been  able  to  achieve  the  experimentally 
reported  restilts*  Our  results  are  due  to  an  improved  method  for  determining 
semiconductor  material  parameters.  We  use  a  theoretical  Monte  Carlo  materials 
simulation  program  to  determine  the  required  parameters.  This  technique 
provides  a  method  for  investigating  varloua  semiconductors  and  examining 
material  potential.  Using  this  technique  we  have  obtained  excellent  agreement 
between  the  theoretical  results  and  experimentally  reported  device 
performance . 

A  MESFET  device  model  has  been  developed  and  used  to  study  the 
high-frequency  operation  of  three-terminal  compound  semiconductor  devices. 
Various  semiconductors  have  been  investigated  Including  GaAs ,  InP,  GalnAs, 
GelnAsP,  and  AllnAs.  Generally  we  find  that  a  high  saturation  velocity  is 
required  for  high  gain  and  a  high  low-field  mobility  is  required  for  low  noise 
figure.  Although  InP  is  capable  of  high  gain  it  suffers  from  high  noise 
performance.  The  optimum  material  for  high  frequency  FET  type  devices  appears 
to  be  GalnAs.  At  40  GHz  a  GalnAs  transistor  is  capable  of  about  2  db  more 
gain  than  a  comparable  GaAs  device.  In  general,  our  studies  indicate  that 
MESFET  devices  with  half-micron  gate  lengths  should  be  useful  to  about  60-70 
GHz.  Quarter  micron  gate  length  devices  could  possibly  provide  useful  gain  as 
high  as  100  Os.  For  operation  above  70  GHz ,  however,  novel  device  structures 
need  to  be  investigated  and  developed. 

High-frequency,  nonequilibrium  transport  effects  have  been  investigated 
with  theoretical  Monte  Carlo  device  simulations.  These  calculations  have 
helped  provide  insight  into  velocity  overshoot  and  ballistic  transport 
phenomena.  We  find  that  ballistic  or  near-ballistic  transport  effects  are  not 
likely  to  result  in  Improved  device  performance.  Our  studies  of  the 
planar-doped  barrier  transistor,  for  example,  reveal  that  this  device  is  not 
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capable  of  performance  comparable  to  that  achieved  from  standard  submicron 
gate  FETs.  Nonequilibrium  effects,  however,  may  lead  to  novel  device 
structures.  For  example,  the  "ensemble  effect"  revealed  in  our  calculations 
(see  paper  9  in  the  appendix)  can  be  used  to  increase  the  average  carrier 
velocity.  By  using  multiple  acceleration  regions  high  average  velocity  over  a 
significant  distance  can  be  maintained.  This  subject  is  currently  being 
examined  in  more  detail. 

Experimentally,  we  have  grown  epitaxial  semiconductor  layers  of  GaAs  by 
MBE  and  OM-CVD,  and  InP  and  GalnAs  by  OM-CVD  and  LPE.  These  layers  have  been 
used  for  both  characterization  and  device  fabrication.  He  have  fabricated 
GaAs  IMPATT  diodes  from  MBE  and  OM-CVD  material  and  GalnAs/InP  devices  from 
LPE  material.  These  diodes  are  being  tested  in  a  waveguide  oscillator  circuit 
and  their  performance  will  be  compared  to  commercial  GaAs  and  Si  diodes.  We 
have  developed  the  facilities  to  completely  RF  characterize  the  devices. 

Circuit  design  work  has  also  been  supported  on  this  contract.  Our 
efforts  have  concentrated  upon  the  development  of  circuit  design  techniques 
that  will  optimize  component  performance  in  terms  of  gain-bandwidth  tradeoffs 
and  parameter  insensitive  design.  This  work  is  summarized  in  the  publications 
attached  to  this  report. 
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Optimum  Semiconductors  for  High-Frequency 
and  Low  Noise  MESFET  Applications* 

J.  M.  Golio  and  R.  J.  Trew 
Electrical  Engineering  Department 
North  Carolina  State  University 
Raleigh,  NC  27650 

ABSTRACT 

An  analytic  MESFET  device  model  has  been  used  to  study  the  optimum 
velocity-field  characteristics  of  materials  that  are  potentially  useful 
for  millimeter-wave  and  low-noise  MESFET  applications.  Materials  of  current 
interest  have  been  characterized  and  compared.  Results  explain  the  relative 
importance  of  parameters  such  as  low-field  mobility  and  saturated  velocity. 
Differences  between  GaAs  and  Si  performance  are  explained  and  a  number  of 
attractive  compound  semiconductors  for  high-frequency  and  low-noise  device 
fabrications  are  indicated. 


*  This  work  was  partially  supported  by  Rockwell  International  and 
partially  supported  by  the  Army  Research  Office,  under  Contract 
No.  DAA  G29-80-K-0080 . 
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The  GaAs  MESFET  has  been  utilized  In  a  variety  of  low-noise  applications 
at  microwave  and  millimeter-wave  frequencies.  Although  GaAs  performance  has 
been  demonstrated  beyond  30  GHz  with  corresponding  minimum  noise  figures  of 
about  4.5  dB  (1,2),  such  operation  appears  to  be  approaching  a  limit  for 
practical  operation  of  GaAs  MESFET  devices.  For  this  reason,  interest  in 
the  use  of  other  compound  semiconductors  is  growing.  Recently,  special 
interest  has  been  given  to  both  InP  <3, 4, 5, 6)  and  Ga  ^  In  ^  As  (7,8,9). 

The  initial  work  with  these  semiconductors  has  been  encouraging. 

The  purpose  of  this  work  is  to  determine  which  materials  are  most 
likely  to  result  in  optimum  millimeter-wave  and  low-noise  MESFET  performance. 

The  frequency  and  noise  properties  of  these  devices  has  been  examined  as  a 
function  of  electron  low-field  mobility  and  high-field  velocity.  To  accomplish 
this,  a  one-dimensional  analytic  model  similar  to  the  one  reported  by  Pucel 
et  al.  (10)  has  been  employed. 

MODEL  ASSUMPTIONS  AND  TECHNIQUE 

The  model  used  in  this  work  assumes  that  the  electron  transport 

properties  of  a  material  can  be  simulated  by  a  two-piece  velocity-field 

relationship.  The  two-piece  approximation  (Figure  1)  is  defined  from  a 

theoretical  velocity-field  characteristic  determines  oy  Monte  Carlo  techniques. 

For  electric  fields  less  than  an  appropriate  saturation  field,  E  , 

b 

the  electron  velocity  is  described  by  a  linear  expression, 

v- p  E. 
o 

For  electric  fields  above  E  the  electrons  move  at  a  constant,  maximum  velocity, 

v  .  The  velocity-field  curves  for  several  materials  are  compared  in  Figure  2. 
s 

Using  this  information,  a  small-signal  equivalent  circuit  can  be  determined 
and  analysed  to  obtain  RF  performance  predictions,  including  power  gain  and 
noise  figure  as  a  function  of  frequency. 
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For  Che  first  part  of  this  study  the  theoretical  model  was  used  to 
determine  the  relative  importance  of  low-field  mobility  and  high-field 
velocity  on  MESFET  performance.  In  particular,  a  one-micron  gate  length 
geometry  was  assumed  and  the  magnitudes  of  these  semiconductor  properties 
were  varied.  Other  material  and  geometric  properties ' of  the  device  remained 
constant.  The  effects  that  low-field  mobility  and  maximum  velocity  have  upon 
device  figures  of  merit  were  thus  determined. 

In  the  second  part  of  the  study,  materials  were  characterized 
according  to  their  Monte-Carlo  predicted  properties.  The  Monte  Carlo  simulations 
used  for  this  work  have  been  well  documented  and  are  in  excellent  agreement 
with  experimental  measurements  (11,12). 

Low-field  mobility  for  the  various  materials  of  interest  is  obtained 
directly  from  Monte  Carlo  data.  Determining  the  appropriate  maximum  velocity 
for  materials  is  more  involved. 

The  Importance  of  developing  a  systematic,  well-justified  technique  for 
determining  v  has  been  discussed  in  earlier  work  (13)  and  in  this  study, 
such  a  method  has  been  developed.  The  method  involves  numerical  determination 
of  the  carrier  transit  time  under  the  gate  using  an  exact  velocity-field 
relationship.  It  is  then  required  that  the  transit  time  calculated  assuming 
a  two-piece  approximated  velocity-field  relationship  be  equal  to  that  of  the 
exact  analysis.  This  method  yields  a  v  for  GaAs  (doped  to  a  level  N  -  lO^cm  ^ 
in  excellent  agreement  with  the  value  obtained  by  Pucel  (10).  Table  1  gives 
values  for  both  low-field  mobility  and  maximum  velocity  determined  in  this  way 
for  several  Important  materials.  Note  that  both  of  these  material  parameters 
will  be  doping  dependent. 
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The  relative  importance  of  low-field  mobility  and  maximum  velocity 
upon  high-frequency  and  noise  performance  was  investigated  by  independently 
varying  these  parameters  and  computing  the  gain-bandwidth  product  and 
minimum  noise  figure.  The  gain-bandwidth  product,  f^,  is  a  commonly  quoted 
figure  of  merit  for  high-frequency  devices  and  can  be  defined  from  first- 
order  considerations  as 

fT  ■  V'2*C8. 

where  g^  is  the  device  trans conductance  and  C^g  is  the  gate-source 
capacitance.  The  minimum  noise  figure  computation  is  developed  after 
Pucel  et  al.  (10)  and  represents  the  minimum  theoretical  noise  of  the 
intrinsic  device.  Since  such  noise  sources  as  circuit  losses  and  defect 
states  are  ignored,  the  minimum  noise  figures  here  are  somewhat  lower  than 
those  measured  for  lym  GaAs  devices  today. 

RESULTS 

The  variations  of  f^  and  NF^  with  low-field  mobility  is  presented 

in  Figure  3.  The  data  for  the  f^  curves  is  for  devices  biased  with 

while  the  NFm^n  data  is  biased  for  minimum  noise  figure.  Mote  that  as  yc 

is  Increased,  the  predicted  device  gain-bandwidth  product  increases  rapidly 

2 

for  Values  below  about  2000  cm  /v«sec.  For  values  in  excess  of  this  limit, 
little  or  no  improvement  is  observed.  This  information* helps  to  explain  why 
the  frequency  response  for  GaAs  is  much  greater  than  that  for  Si,  while 
further  increases  in  mobility  might  yield  little  improvement  in  that  area. 
Note,  however,  that  increased  mobilities  do  continue  to  significantly 
decrease  the  predicted  minimum  noise  figure 
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The  variations  of  gain-bandwidth  product  and  minimum  noise  figure 
with  increasing  maximum  velocity  is  presented  in  Figure  4.  The  increase  of 
with  vffl  is  almost  linear  showing  an  order  of  magnitude  increase  in  f^  with 


*  corresponding  increase  in  v  .  High  maxi, man  velocity  semiconductors  are  12 

o 

therefore  the  most  desirable  for  improved  high-frequency  performance.  It 
should  be  noted,  however,  that  in  the  range  of  maxinum  velocity  for  practical 
semiconductors  (i.e.  1x10^  -  2x1  (T  cm/sec)  very  little  improvement  in  minimum 
noise  figure  is  made. 

Finally,  the  material  parameters  corresponding  to  the  III-V  compounds 
of  Table  1  were  used  to  simulate  a  0.5  m  gate  length  device.  Maximum  velocity 
and  low-field  mobility  were  determined  as  outlined  above,  while  the  remaining 
necessary  parameters  were  determined  as  described  in  reference  (13) .  The 
devices  were  examined  at  two  different  bias  levels:  biased  for  minimum 
noise  figure  and  biased  for  a  drain  current  of  10  mA.  Figure  5  Illustrates 
the  available  gain  of  the  devices  when  biased  for  L^-lOmA  and  the  noise  figure 
when  biased  for  noise  figure.  Referlng  to  Table  1,  it  can  be  seen 

that  ranking  of  materials  by  superior  gain  predictions  corresponds  exactly 
to  the  ranking  by  highest  maximum  velocity,  while  the  materials  which  show 
superior  noise  properties  are  those  with  highest  low-field  mobility. 

COKCLUSIOHS 

In  conclusion,  a  number  of  semiconductor  compounds  have  been  characterized 
and  compared.  The  relative  importance  of  the  transport  properties  as  described 
by  low-field  mobility  and  saturation  velocity  has  been  examined.  The  results 
Indicate  that  for  optimum  high-frequency  response,  high-field  particle 
velocity  is  more  critical  than  low-field  mobility.  Thus,  a  material  such  as 
InP  may  be  favored.  Similarly,  for  optimum  low-noise  performance,  it  is 
low-field  mobility  which  seems  to  be  most  influential.  Therefore,  a  material 
such  as  Ga  4?In  ^As  may  be  preferred. 
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Table 


:  Low-field  mobility,  y . ,  and  saturation  velocity,  vm, 
materials.  Values  are  appropriate  for  1  um  devices 
level  ND  -  1017ca-3. 


for  several 
doped  to  a 


Material 

2 

U0(ca  /v*sec) 

vm(107cm/sec) 

Si 

500 

1.00 

Ca  As 

4660 

1.30 

InP 

3490 

1.84 

Ga.47In.53Aa 

8900 

1.31 

5050 

i 

0.99 

15 


Figure  1  Carrier  velocity  ea  a  function  of  electric  field  for  GaAs 
doped  to  Njj  ■  io*7car3. 


Figure  2  Static  velocity  -  field  characteristics  for  a  number  of 

semiconductors  of  interest.  The  doping  level  is  Nq  -  1017  cm~3. 

Figure  3  Model  predictions  of  gain-bandwidth  product  and  minimum  noise 
figure  as  a  function  of  low-field  mobility.  Other  material  para¬ 
meters  used  are  appropriate  to  GaAs  doped  to  Nq  ■  10l7  cm**3. 


Figure  4  Model  predictions  of  gain-bandwidth  product  and  noise 

figure  as  a  function  of  maximum  velocity.  Other  material  para¬ 
meters  used  are  appropriate  to  GaAs  doped  to  Nj  ■  1017  cm--*. 


Figure  5  Model  results  of  minimum  noise  figure  and  available  gain  as 

a  function  of  frequency  for  several  materials  of  current  Interest. 
The  noise  figure  data  is  for  devices  biased  to  obtain  m-tn-tmimi 
noise  figure,  and  the  available  gain  information  is  for  devices 
biased  at  Iq  -  10  mA. 
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Minimum  Noise  Figure  (dB)  at  10  GHz 
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Electron  transport  in  planar-doped  barrier  structures  using  an  ensemble 
Monte  Carlo  method*’ 

M.  A.  Littlejohn,  R.  J.  Trew,  J.  R.  Hauser,  and  J.  M.  Golio 

Department  of  Electrical  Engineering,  North  Carolina  State  University,  Raleigh,  North  Carolina  27650 
(Received  6  December  1982;  accepted  11  January  1983) 

An  ensemble  Monte  Carlo  method  is  used  to  simulate  electron  transport  through  narrow  high 
field  regions  which  occur  in  planar-doped  barrier  devices.  These  structures  can  achieve  average 
velocities  substantially  higher  than  the  static  drift  velocity  corresponding  to  average  electric  field. 

Increased  velocities  in  these  structures  arise  primarily  from  two  physical  mechanisms.  These 
mechanisms  are  (1)  velocity  overshoot  and  related  ensemble  hot  electron  effects,  and  (2)  an 
ensemble  effect  due  to  the  collecting  nature  of  the  low  field/high  field  boundary,  which  is  not  due 
to  hot  electrons.  Ballistic-like  transport  can  occur  through  narrow  high  field  spikes.  However, 
ballistic  transport  is  not  the  major  contributor  to  high  ensemble  average  velocities  due  to 
scattering  in  regions  adjacent  to  the  high  field  spike.  Ensemble  average  velocities  at  a  given  point . 
in  a  device  are  influenced  by  the  field  distribution  ahead  of  this  point  as  well  as  that  behind  it 

PACS  numbers:  73.40.Lq,  85.30.Mn 


I.  INTRODUCTION 

The  study  of  transport  in  semiconductors  on  short  time  (tss  1 
ps)  and  distance  (x~0.  1  /rm)  scales  has  identified  velocity 
overshoot1"*  and  ballistic  transport5"9  as  potential  mechan¬ 
isms  for  increasing  carrier  velocity.  Such  studies  receive  im¬ 
petus  from  advances  in  lithography  and  materials  growth 
techniques,  as  well  as  the  desire  to  realize  innovations  and 
improvements  in  high  speed,  high  frequency  devices. 

The  planar-doped  barrier  is  a  majority  carrier  rectifying 
structure  which  is  potentially  capable  of  high  speed  applica¬ 
tions  based  on  hot  electron  transport  phenomena.  I0"‘3  Since 
majority  carrier  transport  occurs  over  dimensions  less  than 
0. 1  p m  and  through  narrow  high  field  spikes,  the  planar- 
doped  barrier  can  possibly  exhibit  both  velocity  overshoot 
and  ballistic  transport  Other  ensemble  effects  not  specifical¬ 
ly  due  to  hot  carriers  can  also  enhance  the  dynamic  velocity 
response. 13  The  purpose  of  this  paper  is  to  discuss  transport 
mechanisms  in  planar-doped  barriers  using  data  obtained 
from  ensemble  Monte  Carlo  simulations.  The  results  suggest 
that  the  dominant  transport  effects  are  velocity  overshoot 
and  ensemble  effects  due  to  the  low  field/high  field  boundar¬ 
ies  and  the  field  distribution. 

II.  TRANSPORT  MODEL 

The  details  of  the  ensemble  Monte  Carlo  method  have 
been  described  previously.*  The  capability  to  include  nonun¬ 
iform  spatially  varying  electric  field  distributions  has  been 
incorporated  into  the  Monte  Carlo  procedure.14-13  In  this 
paper,  we  employ  a  series  of  electric  field  “steps,”  which 
represent  the  ideal  field  distribution  in  a  planar-doped  bar¬ 
rier  transistor." 

The  planar-doped  barrier  is  ana  structure 

grown  by  molecular  beam  epitaxy. 10,13  Tbep*  region  is  thin 
enough  to  be  depleted  of  free  carriers.  Also,  the  lightly  doped 
it  regions  are  negligibly  doped  compared  to  thep*  region. 
The  negative  charge  in  thep'*  layer  is  neutralized  by  positive 
bound  charge  in  the  two  heavily  doped  n*  regions.  If  sn 
impulse  charge  density  is  assumed,  the  resulting  step  electric 


field  distribution  leads  to  a  triangular  electron  potential  en¬ 
ergy  barrier  in  the  conduction  band.  Current  flow  through 
the  barrier  is  asymetrical  with  applied  voltage  when  the  two 
n  regions  have  different  widths. 

Two  barriers  can  be  connected  to  form  a  three  terminal 
planar-doped  barrier  transistor.11,12  The  field  profile  for  an 
ideal  planar-doped  barrier  transistor  is  shown  in  Fig.  1.  Ty¬ 
pical  dimensions  are  specified.  The  base  width  is  0.2  pm,  the 
emitter-base  junction  width  is  0.25  pm,  and  the  collector- 
base  junction  width  is  0.45  pm.  The  area  charge  density  is 
5.5  X  10“*  C/cm2  for  the  emitter-base  junction  p*  acceptor 
plane,  and  3.8  x  10“ *  C/cm2  for  that  of  the  collector-base 
junction.  Figure  1  also  shows  the  field  profile  whic  v  results 
from  an  emitter-base  forward  bias  of0.15  V  and  a  •;»"*' -fctor- 
base  reverse  bias  of  0.3  V.  The  equilibrium  batfri-sr  heights 
are  0.2  eV  for  the  emitter-base  junction  and  0.15  eV  for  the 
collector  base  junction. 

111.  SIMULATION  OF  A  PLANAR-DOPED  BARRIER 
TRANSISTOR 

Field  profiles  of  the  form  shown  in  Fig.  1  are  used  in  the 
simulations.  An  ensemble  of  electrons  (typically  2000- 
10  000)  is  released  at  x  «  0.  This  ensemble  is  transported 
across  the  transistor  to  the  collector  boundary.  During 
transport,  the  electrons  undergo  scattering  according  to  the 
normal  GaAs  scattering  mechanisms.4  The  region  0o t<L  is 
subdivided  into  equally  spaced  regions  of  width  Ax  «*  L  /M, 
where  M  is  typically  100.  Smaller  spacings  are  used  for  the 
very  narrow  field  spikes.  Each  particle  in  the  ensemble  is 
tracked  and  the  ensemble  average  velocity  in  each  Ax  bin  is 
accumulated.  A  particle  reaching  the  collector-base  junction 
with  a  kinetic  energy  lower  than  the  collector-base  barrier 
height  is  removed  from  the  ensemble.  A  particle  can  initially 
traverse  the  collector-base  barrier  and  then  be  badcscattered 
into  the  base.  If  such  a  particle  approaches  the  collector-base 
barrier  again  with  less  kinetic  energy  than  the  barrier  height, 
it  is  then  removed  from  the  ensemble.  Thus,  the  ensemble 
average  velocity  of  particles  which  reach  the  collector  is 
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Fie.  1.  Electric  field  distribution  vs  distance  in  an  ideal  planar-doped  bar¬ 
rier  transistor  with  and  without  bias  voltages. 


computed.  Unless  otherwise  stated,  all  calculations  present¬ 
ed  will  be  for  GaAs  with  material  parameters  used  previous¬ 
ly.4 


Figure  2  (curve  (a)]  shows  ensemble  average  velocity  for  a 
planar-doped  barrier  transistor  with  1000  A  base  width  and 
300  A  widths  for  the  emitter  and  collector  high  field  spike 
regions.  For  convenience,  the  indicated  field  polarity  is  posi¬ 
tive  downward  in  this  and  subsequent  figures.  Also,  the 
emitter-base  and  collector-base  barrier  heights  are  0.21  and 
0.05  eV,  respectively  and  then*  regions  illustrated  in  Fig.  1 
are  assumed  to  have  zero  thickness  for  these  simulations.  In 
Fig.  2,  there  are  four  regions  of  interest  The  first  is  the  re¬ 
gion  prior  to  the  emitter-base  high  field  spike  where  the  ve¬ 
locity  rapidly  rises  to  above  2  X 107  cm/s.  The  second  is  the 
region  within  the  emitter-base  high  field  spike  where  the  ve¬ 
locity  rapidly  rises  to  between  (3-4)  x  107  cm/s  and  then  falls 
to  below  2  X 107  cm/s.  The  third  region  is  the  base  where  the 


Tabu  !.  Calculated  values  of  transit  time  (r) ,/»  (Ini'1,  tad  ratio  of  the 
number  of  camera  collected  to  the  number  lost  in  the  beet  [Nc/N,]  tor  two 
values  of  collector -beat  bartici  height  idee)- 

Base  dec  ■» 0. 1 1  eV  dec-0.16eV 


width 

(A) 

t1p»l_ 

(2wrT’ 

(GHi| 

f/c/N, 

tfpel 

(2rrT' 

(OHa) 

Nc/N, 

2000 

mm. 

6.3 

25J 

3.0 

1000 

5.6 

28.4 

6.7 

250 

It 

5.1 

3U 
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velocity  varies  slightly,  and  the  fourth  region  is  the  base- 
collector  junction  where  the  velocity  is  well  above  2X 107 
cm/s  for  over  0.1  /tin.  The  details  of  these  regions  will  be 
discussed  in  the  next  section.  Also,  curve  (c)  in  Fig.  2  shows 
the  ensemble  velocity  for  a  transistor  with  a  base  width  of 
250  A  (all  other  parameters  remain  fixed).  The  results  are 
similar  to  curve  (a),  with  the  velocity  in  the  base  increasing  to 
above  2.5  X 107  cm/s  in  this  case. 

Using  data  of  the  type  shown  in  Fig.  2,  the  transit  time  r 
across  the  device  has  been  calculated  as  a  function  of  base 
width.  Values  of  r  are  shown  in  Table  I  for  two  different 
values  of  collector-base  barrier  height.  Also  shown  in  Table  I 
is  the  value  of  /  =  (2jrr)*‘.  While  /is  not  the  true  cut-off 
frequency  for  the  device,  it  represents  a  lower  bound  to  the 
cut-off  frequency.  The  values  of  r  in  Table  I  properly  reflect 
the  emitter-  and  base-transit  times.  However,  the  actual  col¬ 
lector  transit  time  will  be  shorter.  Thus,  GaAs  planar-doped 
barrier  transistors  should  operate  above  30  GHz  with  base 
widths  less  than  0. 1  /tm.  Table  I  also  shows  corresponding 
values  of  Nc/Nt ,  where  Nc  is  the  number  of  carriers  reach¬ 
ing  the  collector  and  Nt  is  the  number  with  insufficient  ener¬ 
gy  to  cross  the  base-collector  barrier.  The  values  of  Nc/Nt 
can  be  equated  with  the  d-c  common  emitter  current  gain 
based  on  a  simple  transport  model.17  In  the  simulations,  cur¬ 
rent  gain  greater  .than  ten  can  be  achieved  in  planar-doped 


Fie.  2.  Ensemble  average  velocity 
[curve  (a)]  end  electric  field  distribu¬ 
tion  [curve  (b)]  vs  distance  fora  QaAa 
planar-doped  barrier  transistor  with  e 
base  width  of  1000  K.  Curve  (c)  abowa 
ensemble  average  velocity  for  a  device 
with  a  base  width  of 250  A 


J.  Vac.  Set  Tectmot  B,  Vot  1,  Mo.  2,  Apr  .-Juno  IMS 


Ml 


UMajoim  tat-. 


23 


Fra.  3.  Ensemble  evenfe  velocity 
[cum*  (*).  (b).  tad  (e)]  and  electric 
Add  distribution  [com  (d)]  n  die- 
tuot  for  a  OaanliiauAi  planar- 
dopad  barriar  translator.  Results  for 
three  values  of  intervalley  energy  se¬ 
paration  (dn)  are  shown. 


barrier  transistors  with  base  widths  less  than  about  1000  A. 

Figure  3  shows  the  ensemble  average  velocity  for  a  struc¬ 
ture  with  material  parameters  the  same  as  Ga^,  In^jj  As14 
for  three  different  values  of  intervalley  energy  separation 
dfrt .  The  device  parameters  are  the  same  as  those  in  Fig.  2 
with  base  width  *  1000  A.  Until  the  final  field  step  in  the 
collector  region,  the  transport  in  this  material  is  not  drasti¬ 
cally  different  from  that  of  GaAs.  In  the  collector,  the  veloc¬ 
ity  transition  and  final  velocity  value  is  greatly  influenced  by 
At rL.  This  illustrates  that  transport  in  this  region  is  in¬ 
fluenced  by  velocity  overshoot.17  Also,  increased  speed  and 
frequency  performance  can  be  expected  in  pianar-doped  bar¬ 
rier  transistors  fabricated  from  such  central  valley  dominat¬ 
ed  materials  as  GalnAs. 11 

IV.  DISCUSSION  OF  TRANSPORT  MECHANISMS 

In  the  region  to  the  left  of  the  emitter-base  high  field  spike 
there  is  a  large  increase  in  the  ensemble  average  velocity. 
This  velocity  increase  has  been  shown  to  be  due  to  the  col¬ 
lecting  nature  of  the  emitter-base  high  field  boundary. 11  A 
carrier  which  crosses  into  the  region  x>0.2  pm  has  a  low 
probability  of  being  scattered  back  into  the  low  field  region 
x  <0.2  pm.  The  probability  for  such  backscattering  de¬ 
creases  as  the  barrier  height  of  the  emitter-base  high  field 
spike  region  increases.  This  is  illustrated  in  Fig.  4,  which 
shows  the  ensemble  average  velocity  for  high  field  spikes 
with  barrier  heights  of  0.06,  0.12,  and  0.3  eV.  The  Monte 
Carlo  data  points  for  an  ideal  collecting  barrier  an  also 
shown.  These  data  points  wen  generated  in  a  separate 
Monte  Carlo  simulation  for  which  all  particles  reaching 
x  ••  0.2  pm  wen  removed  The  0.06  eV  barrier  is  not  quite 
an  ideal  collecting  boundary,  while  the  0. 12  and  0.3  eV  bar¬ 
riers  an  essentially  ideal.  Our  simulations  have  shown  that  a 
barrier  height  of  0.1  eV  or  greater  acts  as  an  ideal  collecting 
boundary  in  pianar-doped  barrier  structures.  For  barrier 
heights  greater  than  0. 1  eV,  then  an  very  few  particles  at  the 
boundary  which  have  negative  velocities. 

The  data  shown  in  Fig.  4  also  relates  to  transport  in  the 
base  region  of  the  pianar-doped  barrier  transistor.  The  en¬ 


semble  velocities  reach  a  maximum  of  between  2.7  X 107  and 
3.8  X  107  cm/s  within  the  high  field  region,  depending  on  the 
field  value.  Also,  the  ensemble  velocities  decrease  rapidly 
before  the  high-to-low  field  transition  occurs.  Outside  the 
high  field  region,  the  ensemble  velocities  for  step  field  values 
of  20  and  40  k V/cm  stay  above  the  steady  state  velocity  ( Va ) 
for  the  1  kV/cm  background  field.  This  behavior  is  unex¬ 
pected  based  on  the  expected  quick  acceleration  propsoed 
for  the  pianar-doped  barrier  transistor."  In  this  picture,  the 
emitter-base-high  field  region  is  kept  small  enough  (~150- 
300  A)  so  that  collisions  are  unlikely.  These  electrons  are 
then  halliitically  launched  into  the  base  region  where  they 
are  expected  to  travel  up  to  0.3  pm  before  losing  a  large 
fraction  of  momentum.  "  H  While  Fig.  4  shows  that  the  ve¬ 
locity  in  the  base  is  higher  than  Vu  for  fields  of  20  and  40  kV/ 
cm,  and  that  it  can  stay  high  for  a  substantial  distance  (Le.,  a 
few  hundred  angstroms),  the  ensemble  velocities  are  lower 


Fra.  4.  Ensemble  ever***  velocity  *t  distance  for  high  Add  (pitas.  The 
electric  Add  Iwtsnswiss  Air  tta  (pitas  in  20  kV/em  [cam  (all.  40  kV/cm 
[cam  (b)).  sad  IOOkV/Cta[carvc|c|lovsra1tatarsof300ATtavchissof 
Vm  shown  tndioesee  the  steady  stats  avenge  velocity  in  tta  t  kV/cm  bock- 
froaad  Arid. 
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Fkj.  3.  Ensemble  avenge  kinetic  energy  v»  distance  for  a  Add  ipike  of  70 
kV/cm  over  }Q0  k  with  and  without  inoervalley  scattering. 


by  a  factor  of  2  than  would  be  anticipated  based  on  the  quick 
acceleration  argument. 

As  proposed  in  the  quick  acceleration  model,  our  simula¬ 
tions  show  that  the  ensemble  can  be  accelerated  baJlistically 
through  the  emitter-base  high  field  region.  This  is  illustrated 
in  Fig.  5,  which  shows  ensemble  average  kinetic  energy  ver¬ 
sus  distance  through  a  70  kV/cm,  300  A  high  field  spike. 
When  the  £(111)  conduction  band  minimum  is  removed 
from  the  simulation,  the  ensemble  average  kinetic  energy  at 
the  edge  of  the  high  field  spike  is  0.21  eV.  This  is  the  energy 
gained  ballisbcally  from  an  electron  traveling  300  A  in  a  70 
kV/cm  field.  The  ensemble  actually  travels  balhstically  for  a 
distance  of  about  330  A  since  the  energy  does  not  change 
abruptly  when  the  field  changes.  However,  when  the  upper 
(satellite)  valley  is  included,  the  ensemble  is  not  accelerated 
to  the  ballistic  energy,  but  reaches  a  maximum  of  about  0. 16 
eV.  In  both  cases  the  kinetic  energy  decay  occurs  over  a 
relatively  long  distance,  indicating  a  large  energy  relaxation 
time. 

In  Figs.  4  and  3,  the  energy  which  can  be  ballisbcally  im¬ 
parted  to  an  electron  from  the  field  spikes  are  less  than  0.33 


eV,  the  energy  separatum  between  the  central  (F)  and  satel¬ 
lite  {L )  valleys  (deri )  in  GaAs.4  Questions  then  arise  as  to 
why  intervalley  scattering  affects  the  velocity  and  energy 
shown  in  Figs.  4  and  3,  respectively.  The  energy  shown  in 
Fig.  3  is  total  ensemble  average  energy  for  all  particles  in 
both  r  and  L  valleys.  We  have  previously  shown  that  the 
distribution  of  carriers  in  the  T  valley  extends  to  much  high¬ 
er  energies  than  the  ensemble  average  kinetic  energy  for 
fields  well  above  the  threshold  field.1*  This  energy  redis¬ 
tribution  occurs  very  rapidly  (t«0.S  ps)  when  electrons  are 
injected  into  a  high  field  region.  Thus,  for  the  cases  illustrat¬ 
ed  in  Figs.  4  and  5,  the  energy  distribution  in  the  F  valley 
extends  to  energies  well  above  Aft.  .  intervalley  scatter¬ 
ing  is  important  in  spite  of  the  fact  that  the  ensemble  average 
kinetic  energy  is  well  below  . 

For  the  fields  and  dimensions  considered,  the  dominant 
scattering  mechanisms  in  the  GaAs  planar-doped  barrier 
transistor  are  polar  optical  scattering,  ionized  impurity  scat¬ 
tering,  and  intervalley  scattering.1’  Intervalley  scattering  is 
highly  isotropic,  while  polar  optical  scattering  and  ionized 
impurity  scattering  are  primarily  small  angle  scattering 
events. ,9J0  If  isotropic  scattering  becomes  a  significant  frac¬ 
tion  of  the  scattering  events  in  the  base  there  will  be  a  large 
number  of  electrons  which  are  backscattercd  toward  and 
into  the  emitter-base  high  field  spike  region.  Since  most  of 
these  carriers  do  not  cross  the  barrier  back  into  the  emitter 
region,  their  negative  velocities  lower  the  ensemble  average 
velocity  in  both  the  base  region  and  in  the  emitter-base  high 
field  spike  region.  This  is  apparently  the  mechanism  for  ve¬ 
locity  reduction  in  these  regions  (Fig.  4).  The  increase  in 
intervalley  scattering  is  enhanced  because  the  energy  relaxa¬ 
tion  time  is  large.  Since  the  ensemble  average  energy  in  the 
base  region  remains  high  over  a  substantial  distance,  the 
scattering  rate  for  intervalley  scattering  remains  high  until 
the  energy  decays  sufficiently.  The  scattering  rate  decreases 
as  the  energy  decreases,  which  further  prolongs  the  energy 
relaxation.10 

Figure  6  further  illustrates  these  points.  The  ensemble 
average  velocity  is  shown  for  a  20  kV/cm  step  field  occurring 
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at  x  =•  0.1  fj m  and  for  a  field  spike  of  20  kV/cm  at  x  =*  0.1 
tim  with  a  width  of  0.1  fim.  For  x  >  0. 1  fit n,  the  ensemble 
average  velocity  for  the  step  field  is  larger  than  the  ensemble 
velocity  for  the  field  spike.  The  potential  barrier  for  the  field 
step  is  0.2  eV  across  the  0. 1  nm  distance,  which  is  more  than 
adequate  to  prevent  backscattered  particles  from  crossing 
the  barrier  into  the  region  x<0.1/sm.  However,  the  poten¬ 
tial  barrier  decreases  linearly  for  the  field  step,  and  is  0.4  eV 
at  *■  0.3/*mand0.6eV for*:*  0.4/sm.  Thus,  carrienin  the 
region  x>0.2  fim  ate  more  readily  backscattered  into  the 
region  0.1<jc<0.^tm  for  the  field  spike  than  for  the  step 
field.  Also,  Fig.  7  shows  the  ensemble  average  energy  for 
these  two  field  profiles.  The  ensemble  avenge  energy  is 
greater  for  the  step  field  in  the  region  x>  0.2  because  the 
carriers  gain  additional  energy  from  the  higher  field.  The 
decay  to  steady  state  is  faster  because  the  amount  of  energy 
loos  required  for  the  step  field  distribution  to  reach  steady 
state  is  less  than  for  the  field  spike. 

Thus,  isotropic  scattering  plays  a  major  role  in  real  apace 
transfer  devices  by  reducing  the  ensemble  average  velocities 
in  the  base  region  and  in  the  emitter-base  high  field  spike 
region.  The  effects  of  isotropic  scattering  on  transport  in  real 
space  transfer  devices  and  other  quick  acceleration  schemes 
cannot  be  eliminated.  However,  to  a  first  approximation,  its 
effects  on  velocity  reduction  in  the  base  region  can  be  mini¬ 
mized  by  limiting  the  energy  which  can  be  imparted  halHtti- 
cahy  in  the  emitter-base  high  field  spike  region  to  less  than 
about  one-half  the  mtervalley  energy  separation.  Thus, 
qEd  <  0.3d rc ,  where  £  is  the  emitter-base  accelerating  field 
and  d  is  the  width  of  the  accelerating  field  region.  This  was 
the  case  for  each  value  of  An  in  Fig.  3,  and  is  the  reason 
there  that  the  emitter-  and  base-region  transport  ate  not  af- 
faded  by  d  rt  It  is  important  to  realize  that  transport  of  a 
distribution  of  carriers  in  a  device  (ensemble  transport)  is 
influenced  by  the  field  distribution  both  behind  and  ahead  of 
the  mam  diplacement  of  the  distribution.  This  is  a  unique 


ensemble  effect  which  is  not  predicted  using  single  particle 
ballistic  or  scattering  models. 


V.  CONCLUSIONS 

The  GaAs  planar-doped  barrier  transistor  is  a  three  termi¬ 
nal  device  capable  of  operating  at  frequencies  above  30  GHz. 
Ballistic,  transport,  velocity  overshoot,  and  ensemble  effects 
all  play  a  role  in  the  physics  of  this  device.  Based  on  the 
results  presented  here,  factors  which  should  be  incorporated 
in  device  design  are  (a)  the  emitter-base  barrier  height  should 
be  greater  than  about  0. 1  eV,  (b)  the  base  width  should  be  less 
than  about  1000  A,  (c)  the  product  of  field  times  distance  in 
the  emitter-base  field  spike  should  be  less  than  one-half  the 
potential  due  to  intervalley  energy  separation,  and  (d)  mate¬ 
rials  with  high  intervalley  Energy  separations  should  be  used 
to  increase  the  velocity  in  the  collector-base  junction. 
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A  Broad- Band  Amplifier  Output  Network  Design 

ALFY  N.  RIDDLE  and  R.  J.  TREW,  member,  IEEE 


jUma-A*  — lyto  4tdgn  —tori  for  a  lossy  [So  ro—inThu 
ntfwock  is  prosoMod  and  tto  adroaH|o»  of  loopy  —works  wo  tote— *4. 
Exafin  of  nm-uifo  nopMIsn  mtog  FITs  sod  Mpohr  n —Hoots  wo 
pmsMod  to  stow  too  f— fbtoty  of  tots  pwtlndw  otrworfc  ia  low  powor 
a— Ml*r  dosffns,  Ttoso  actors  ofxsia  ptm  of  1&4±R5  dB  wfto  a 
23-dB  —ri—  oobt  fig—  fa  too  4.9-4.0-GHx  fri—acy  range  aad 
145=  1.2  «  wit*  »  ai  os  ha  at  iapw  VSWR  of  1. flit  soar  too  1.9-2.9-GHz 
frequency  range*  raspnethefy. 

I.  INTRODUCTION 

Lossy  gain-compensating  output  networks  can  provide  lower 
input  reflection  coefficients,  t  lower  amplifier  noise  figure,  and  a 
more  predictable  amplifier  design  (1).  The  resistive  nature  of  this 
type  of  network  may  also  improve  amplifier  stability  and  distor- 
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tion  by  reducing  standing  waves  within  the  amplifier.  Although 
lossy  broad-band  gain-compensating  networks  are  often  used 
[I]-[4],  explicit,  analytic  design  techniques  for  these  networks 
have  not  been  reported. 

This  paper  presents  an  output  circuit  design  based  upon  a  * 
matching  network  combined  with  a  bandpass/bandstop  diplexer. 
As  a  result  this  network  contains  boau  th»  drain  supply  induc¬ 
tance  and  the  dc  blocking  capacitor,  which  am  needed  in  any 
output  network,  as  integral  elements.  Explicit  fonmaici  for  the 
element  values  of  this  network  are  derived  and  presented  shortly. 
This  method  is  different  from  previous  methods  [1]  because  this 
technique  allows  the  device  load  to  approach  SO  0  as  frequency  is 
decreased.  However,  this  network  is  similar  to  that  used  in  [2] 
and  [3)  in  that  the  supply  voltage  ia  inserted  through  a  quarter- 
wavelength  shunt  stub  and  a  series  resistor.  This  method  results 
in  greater  stability  and  usability  for  the  amplifier. 

A  bandpass/bandstop  diplexer  is  more  useful  than  a  simple 
low-pass/higb-pass  diplexer  because  it  provides  an  exact  match 
at  one  frequency,  and  an  arbitrary  amount  of  attenuation  (limited 
only  by  network  dement  Q't)  at  any  frequency.  Diplexing  net¬ 
works  may  be  used  in  ether  input  or  output  networks  depending 
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on  noise  figure,  power  output,  stability,  and  other  amplifier 
constraints.  This  paper  considers  the  characteristics  of  a  resonant 
diplexer  in  the  output  network. 

U.  Netvoux  Design 


The  output  network  is  designed  in  two  parts.  First,  a  w  network 
matches  the  transistor  output  impedance  to  50  Q  at  the  highest 
frequency  of  interest  fH.  The  insertion  loss  of  this  network  ILr,  at 
the  lowest  frequency  of  interest  fL,  is  calculated  and  subtracted 
from  the  total  loss  desired  for  gain  compensation  ILT  to  give  the 
loss  required  of  the  diplexer  ILd.  This  may  be  expressed  as 

/ M/l)=  H-s(/l)+  IL,(fL)=MAG{fL)-MAG{f„)  (1) 

where  ILr(fH)=0=  ILd(fH)—  IL^/h),  and  MAG(f)  is  the 
frutiffwim  available  gain  of  the  device  and  input  network  at  the 
specified  frequency.  In  this  paper,  insertion  loss  is  in  decibels, 
and  is  with  respect  to  the  maximum  power  available  from  a  50-0 
source.  Tins  simple  two-point  method  of  compensation  works 
wen  for  the  gradual  gain  versus  frequency  slopes  of  single-stage 
amplifiers.  The  diplexer  center  frequency  is  /„,  and  the  diplexer 
bandwidth  is  determined  by  the  required  loss  ILd,  at  fL. 

A  standard  parallel  RC  network  equivalent  circuit  is  derived 
for  the  transistor  output  impedance.  This  circuit  may  be  derived 
in  a  conventional  manner  from  Sn(fH)  in  the  unilateral  case, 
from  the  conjugate  of  the  output  match  reflection  coefficient  for 
mifitwmfB  noise  at  /„,  at  from  the  conjugate  of  the  complete 
match  reflection  coefficient  at  f„.  In  the  *  matching  network 
shown  in  Fig.  1,  R0  and  C,  are  derived  from  the  transistor  output 
admittance.  It  is  worth  noting  that  L2  consists,  at  least  in  part,  of 
the  hooding  wire  to  the  transistor  chip  from  the  microstrip 
circuiL  The  equations  for  this  matching  network  are  derived  by 
separating  the  network  at  node  A  in  Fig.  1  and  writing  expres¬ 
sions  fev  the  impedances  looking  towards  R0  and  RL.  Real  and 
imaginary  (one  direction  being  corqugated)  parts  are  equated  to 
obtain  a  match.  Thus,  we  have 


Qo- 

uh*  o£i 

(2) 

,  1,/J 

UW*  gCj  * 

[> 

*  +  ,H] 

(3) 

C,« 

Ql 

uhRl 

(3a) 

Lj  = 

RjC, 

(4) 

i+e« 

l  +  tfi. 
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Fig  3.  Distributed  fora  of  output  network  lor  laser  gain  coatipcnaaiion. 


where  X2  ~  uHL2.  Since  the  diplexer  proems  a  constant  50  Q  to 
this  network,  this  insertion  loss  is  affected  mainly  by  changes  in 
the  elements  of  the  transistor  output  equivalent  circuit.  The  FET 
design  shown  later  suffered  very  little  from  an  85-percent  change 
in  iu  output  conductance  across  the  operating  frequency  range. 
In  these  situations  an  average  output  conductance  weighted  to¬ 
ward  the  top  of  the  frequency  range  may  be  useful  The  mismatch 
of  the  v  network  at  jL  has  not  been  found  to  disturb  the  diplexer' s 
characteristics.  The  desired  diplexer  loss  ILd,  at  fL  determines  the 
diplexer  bandwidth  by  the  relation 


gw* - — — ^ - 

//rA(10,t'/,#-l)‘/J 


(7) 


where-*  *2* /„. 

The  insertion  loss  for  this  network  may  be  derived  by  describ¬ 
ing  the  network  as  an  A  BCD  matrix  seen  from  RL.  The  matrix 
elements  may  be  converted  to  5-panmeters  via  the  relation  [5] 


_ _ 

A(ju)R9+  D(jm)RL  +  8(j*)+C(jv)R0RL 


(5) 


and  expressed  as  an  insertion  loss  in  decibels  by 

/^.,(y«)»l01og|5JI(y«)|, 


which  cooies  from  the  expression 


IL4  *  lOlog- 


1  + 


»"\  /m/l  /  J 


(7a) 


where  ILd  is  a  positive  number  and  the  bandwidth  BW  is 
normalized.  This  equation  is  a  rearranged  form  of  the  insertion 
loss  for  a  simple  resonant  network.  The  equations  presented 
below  for  the  diplexer  in  Fig  2  may  be  found  in  any  network 
synthesis  text  [6] 
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Pit.  7.  HFET-1 101  4.0-6.0-CHz  lowooiM  duracttrittia.  (« 

<b)|S„|.«:>|Sn|. 


actual  deagni  with  optimized  dement  values  in  pares  thesis.  Figs. 
5  an d  7  show  the  characteristics  of  a  two-stage  cascade  of  each 
amplifier  before  and  after  optimization,  with  the  solid  lines 

The  MRF901  amplifier  was  designed  for  minimum  input  re¬ 
flection  coefficient  and  maximum  Oat  gain  over  the  1.0  to  2.0-CHi 
frequency  range.  The  equivalent  circuit  for  this  device  was  based 
on  a  complete  coagugaie  match  at  1.0  and  2.0  GHz.  this 

transistor's  input  circuit  appeared  as  a  dominantly  inductive 
series- resonant  circuit,  an  iterative  Smith  Chart  procedure  was 
need  to  detign  the  input  network.  The  output  network  was 
de  tignad  according  to  die  aquations  of  the  previous  section.  Since 


the  MAG  of  the  MRF901  was  8.0  dB  at  10  GHz  and  13.8  dB  at 
1.0  GHz,  5.8  dB  of  insertion  lose  was  required  of  the  output 
network  at  1.0  GHz.  This  insertion  lorn  was  oompoaed  of  1.05  dB 
of  w  network  loss  and  4.75  dB  of  diplacer  loss.  The  diplater  loss 
of  4.75  dB  determined  e  bandwidth  of  1.065.  Figs.  4  end  5 
indicate  that  this  design  was  easily  cascaded  and  required  no 
optimization.  Lossless  gain-compensating  output  networks  were 
tried  with  this  combination  of  device  and  input  network  with  lees 
predictable  results.  Figs.  5(a),  (b),  and  (c)  show  of 

£  1.4  dB  of  gam  variation  around  16.8-dB  gain,  1.78:1  input 
VSWX,  and  61:1  output  VSWR. 

The  HFET-1 101  amplifier  was  designed  for  minimum  aoitc 
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figure  in  the  4.0  to  6.0-GHz  frequency  range  with  a  maximum  flat 
gain.  The  analytic  approach  used  in  [2]  was  the  basis  of  the  input 
network  design.  This  input  network  was  designed  to  match  the 
conjugate  of  the  optimum  source  impedance  for  low  noise  at  5.0 
GHz.  Interpolation  was  used  to  derive  this  impedance  from  the 
published  data  at  4.0  and  £.0  GHz.  The  network  element  values 
were  derived  from  an  «=3,  0.01-dB  ripple,  Chebyshev  low-pass 
prototype  [8].  The  first  impedance  inverter  (Ku)  was  omitted 
from  the  final  design  because  its  50.75  Q  impedance  contributed 
little  to  the  overall  match. 

The  transistor  output  equivalent  circuit  was  derived  from  the 
conjugate  of  the  reflection  coefficient  for  a  matched  output  with 
minimum  noise.  The  unilateral  gain  approximation  as  described 
in  [9]  was  used  to  find  the  maximum  amplifier  gain  at  4.0  and  6.0 
GHz.  The  source  reflection  coefficient  derived  from  the  conjugate 
of  the  transistor  input  equivalent  circuit  at  4.0  and  6.0  GHz  was 
used  in  this  unilateral  gain  approximation.  This  required  3.55  dB 
of  loss  at  4.0  GHz.  A  maximum  gain  of  8.45  dB  at  6.0  GHz  was 
expected  for  the  amplifier.  The  diplexer  network  needed  to 
provide  2.85  dB  of  the  3.55-dB  loss,  and  so  a  normalized  band¬ 
width  of  0.864  was  used.  Although  the  cascade  of  two  identical 
single-stage  amplifiers  provided  acceptable  results  without  adjust¬ 
ment  of  any  of  the  circuit  elements,  optimization  improved  the 
gain  flatness  and  output  reflection  coefficient.  The  noise  figures 
of  2.1  dB  and  4.0  GHz  and  235  dB  at  6.0  GHz  were  changed  to 
1.9  dB  at  4.0  GHz  and  248  dB  at  6.0  GHz  during  optimization. 
The  relatively  high  VSWR  of  this  amplifier's  input,  as  shown  in 
Pig.  7(b),  is  due  to  providing  an  optimum  source  impedance  for 
ooise  minimization. 

The  formulas  presented  in  this  paper  are  useful  whenever  the 
equivalent  circuit  of  the  device  is  a  parallel  RC  network.  There¬ 
fore,  bipolar  transistor  amplifiers  operating  well  below  the  device's 
fT  may  contain  the  circuit  just  presented  as  an  input  matching 
network. 


(61  F.  F.  Kuo.  Network  Analysis  ami  Synthesis.  New  York:  Wiley.  1966.  p. 
406. 

(7)  George  L.  Ragan.  Microwave  Transmission  Grants.  New  York:  Dover. 
1969.  p.  617. 

(8)  Anatol  L  Zverev.  Handbook  of  Filler  Synthesis.  New  York:  Wiley.  1967. 
p.  315. 

(9)  S- Parameter  Design.  Hewlett-Packard  Company  HP  Application  Note 
154.  Apr.  1972.  p.  25. 


IV.  Conclusion 

A  simple  output  network  was  shown  to  simultaneously  provide 
gain-compensation  and  a  predictable  amplifier  design.  The  feasi¬ 
bility  of  this  design  method  was  demonstrated  by  cascading  two 
identical  single-stage  amplifiers  and  calculating  the  total  ampli¬ 
fier  5-parameters  before  and  after  optimization.  Although  the 
input  and  output  networks  were  designed  by  treating  the  ampli¬ 
fier  as  if  these  networks  did  not  interact,  the  actual  results  agreed 
well  with  simple  theory.  Explicit  formulas  for  the  design  of 
lumped  and  distributed  output  networks  were  presented. 

Acknowledgment 

The  excellent  work  of  typist  Ms.  N.  Tyson  is  greatly  appreci¬ 
ated  by  the  authors. 

References 

(I]  Nod  Marshall.  "Optunmng  multi-stage  amptihen  (or  low-noue."  Micro - 
wooes,  pp.  62-64.  Apr.  1974.  pp.  60-64.  May  1974. 

(21  C.  A.  Uachfi  and  R.  L.  Tinman.  "Design  sad  performance  of  microwave 
amplifiers  with  GaAs  Scboitky-Galc  field-effect  traneuton."  IEEE  Trans. 
Microwave  Theory  Tech.,  vol.  MTT-22.  pp.  510-517.  May  1974. 

01  D.  P.  Hontbuckle  and  L.  J.  Kublman.  Jr.,  "Broad-Band  medium-power 
ampUficaiten  in  the  2-124-GHz  range  with  GaAs  MESFETs."  IEEE 
Trans.  Microwave  Theory  Tech.,  vol  MTT-24.  pp.  338-  342,  June  1976. 

(4)  Charles  A.  Uechti,  “Microwave  Odd-effect  trwssistors- 1976,"  IEEE  Trans. 
Microwave  Theory  Tech„  vol  MTT-24,  pp.  279-29B.  June  1976. 

(5]  M.  A.  R.  Guneton.  Practical  Moat *  Algebra.  New  York:  Elsevier.  1970. 
p.  in. 


PARAMETER  INSENSITIVE  MATCHING  CIRCUITS 
FOR  LOU  COST  INTEGRATED  CIRCUITS* 


32 


A.  N.  Riddle  and  R.  J.  Traw 
Electrical  Engineering  Daparcaanc 
North  Carolina  Stata  University 
Raleigh,  North  Carolina  27650 


(919)  737-2336 


ABSTRACT 


A  new  theory  for  increasing  tha  tolerance  of 
latching  networks  to  load  variaelona  is  presented. 
This  theory  is  baaed  on  Hatching  the  anglaa  of  tha 
changes  in  Sn(f)  dua  to  frequency  and  load  variations, 
and  its  use  can  double  allowed  device  tolarancas  in 
many  cases.  As  a  byproduct,  this  theory  also  shows  why 
Chebyshev  Hatching  filters  have  poor  tolerances.  Exam- 
ples  are  given  showing  am  auch  as  an  order  of  aagnltude 
lHproveaenc  in  circuit  tolerances. 


INTRODUCTION 

The  manufacture  of  aicrovave  amplifiers,  sixers 
and  oscillators  mist  Include  component  and  device  tol¬ 
erances  in  che  design  procedure -if  lower  cost  and  more 
reliable  circuits  are  to  be  produced..  A  new  theory  for 
aaxialxlng  che  tolerance  of  Hatching  necworks  to  load 
variations  is  presented  here.  The  tolerance  of  Su(f) 
to  variations  in  the  load,  or  any  ocher  component,  is 
maximised  when  che  variaelona  causa  a  change  in  Stl(f) 
that  Is  perpendicular  to  the  direction  of  Su<f). 
Knowledge  of  che  angles  of  the  differential  changes  in 
Sn(f)  with  frequency  and  load  variations  allow  che 
design  of  Hatching  networks  with  a  aaxiHum  tolerance 
over  a  band  of  frequencies. 

The  concern  here  is  to  mexiHise  che  colerance  of  a 
circuit  to  variations  in  che  active  device  used.  For 
che  lossless  networka  discussed  hare,  it  has  been  found 
that  che  aagnlcude  and  frequency  variation  of  the  input 
reflection  coefficient  determines  che  senslclvlty  to 
load  variations.  Slnple  formulae  are  presented  for 
both  differential  and  large  change  sensitivities  in 
terms  of  S-paraawters.  These  formulae,  along  with  the 
nacure  of  Hatching  networka  allows  us  co  establish 
worst  case  tolerances,  and  show  how  matching  networks 
capable  of  doubling  che  cyplcal  worst  case  load  tol¬ 
erance  may  be  designed.  Also,  tha  rapid  changes  with 
frequency  due  to  che  ripples  of  a  Chebyshev  network 
will  be  shown  to  decrease  les  colerance  relative  co 
what  la  possible  with  Hatching  networks  using  a  flat 
mismatch  equal  to  che  peak  of  the  Chebyshev.  Simula¬ 
tions  of  narrowband  end  broadband  feedback  FET  ampli¬ 
fiers  illustrate  chase  conclusions. 

THEORT 


By  using  an  ABCD  aatrlx  representation  of  s 
Hatching  network,  we  aay  see  Che  transformation  from 
che  load  laaiccance  co  tha  input  reflection  coefficient 
as  two  bilinear  transformations.  This  aeans  thee  since 
changes  in  che  real  and  iaaglnery  parts  of  the  load  are 
orthogonal,  resulting  changes  in  the  input  reflection 
coefficient  will  be  orthogonal.  Since  we  usually  want 
to  increase  che  tolerance  of  tha  circuit  to  reactive 
variations  in  che  load,  we  aay  exploit  this  ortho- 
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gonallty  co  Increase  our  colerance.  The  following 
equation  will  prove  useful,  and  is  derived  assuming 

2t<0  -  «  +  J  *  Cf): 

SSu(f).  ,  S2l(f)  S12(f)  (1) 

(f)  2  zL  (f) 

Figure  I  shows  the  center  of  che  reflection  coefficient 
plane  wich  a  circle  enacriblng  che  maximum  tolerable 
reflection.  I*-..  A  point,  Su(f0)  is  located  Inside 
chls  circle  and  represents  s  point  on  the  curve  of 
input  reflection  coefficient  versus  frequency  for  the 
amplifier.  Increases  and  decreases  in  che  device  reac¬ 
tance.  typically  the  input  capacitance,  will  cause  the 
point  S ^ ( f 0)  co  nova.  This  differential  movement  is 
given  exactly  by  aquation  (1)  for  any  network  descrlb- 
able  via  S-paraaeters. 


Figure  l.  Reflection  plane.  The  circle  of  aaxlmua 

tolerable  reflection,  r,,,*,  and  che  circuit 
reflection  coefficient  at  f0,  Su(f0)  are 
shown.  The  dashed  and  alternating  lines 
show  Che  minimum  and  maximum  absolute 

changes  in  S^Cf,,),  respectively. 

The  alternating  line  of  Figure  1  shows  the  maximum 
symmetric  tolerance  abouc  the  point  Su(f0).  Tha 
dashed  line  of  Figure  1  shows  an  asyaamtrlc  tolerance 
which  for  practical  purposes  should  be  limited  to  its 
maximum  syweerlc  variation.  This  shows  how  variations 
in  the  device  capacity  will  have  a  minimum  tolerance  if 
they  cause  radial  excursions  in  Sn(f0)  and  a  maxi¬ 
mum  tolerance  if  they  cause  excursions  of  S j. 2 ( f 0 )  along 
the  alternating  line  of  Figure  1.  As  e  reminder,  the 
tolerance  In  the  device  capacity  would  be  that  percent¬ 
age  change  which  caused  S|,i(fo)  to  lie  on  the  rMX 
circle. 

The  next  step  toward  a  more  tolerant  network  lies 
in  the  expression  of  the  change  in  S22(f)  with  fre¬ 
quency  as  given  by  the  equation 


1C  *■  «■ 

l  Si 


ssuCO  .  .  s21(f>sl2(f)  .  »x 
ST~  iiiTTT  Tt 
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where  A,  1.  C.  and  D  are  Che  paraaeters  of  che  ABCS 
■atria  for  che  input  matching  network.  Equation  (2) 
ahowe  that  for  all  low  0  or  resistive  Batching  networks 
che  absolute  direction  of  8Sn(f)/8f  1*  exactly  the 

ease  as  that  of  8Sn(f)/8X.  depending  only  on  the  sign 
of  8X/8 f.  For  chess  neeworks  react ire  variations  In 
Che  load  have  their  aexlaua  tolerance  when  the  Sn(f) 
le  a  constant  over  the  passband.  l.e.  the  Matching 
filter  has  a  flat  aitaaceh  response.  While  practical 
networka  are  not  resistive  and  rarely  low  0.  equation 

(2)  leads  to  che  conclusion  that  matching  networks 
should  provide  a  flat  mismatch  ac  least  over  che  center 
of  the  band  and  snould  not  have  large  ripples  such  as  a 
Chebyahev  filter.  These  conclusions  will  be  elaborated 
on  via  examples  on  Chebyahev  and  Butterworth  notching 
filters,  as  wall  as  a  broadband  FET  feedback  aapllfler. 

Another  result  of  this  study  is  given  by  equation 

(3)  which  gives  the  aexlaua  percentage  tolerance  of  a 
series  or  parallel  load  capacitance  assuming  chat 
SH(f0)  •  0,  che  matching  network  is  lossless,  the  net¬ 
work  is  also  reciprocal,  and  che  percentage  change  In 
the  capacitance  Is  smell: 

X  tolerance  -  ™  #  (3) 

Of  • 

where  Of  Is  the  Q  of  the  load  at  fa. 


Figure  2.  Idealised  bandpass  aacchlng  filter  for  tast¬ 
ing  tolerances  to  load,  *L(f )  and  circuit 
variations.  Sea  Table  I  and  Figures  3  and  4. 

RESULTS 

The  network  of  Figure  2  wee  used  to  exaaina  Che 
tolerance  of  Ideal  Chebyahev  and  Butterworth  responses 
to  variations  In  che  load  capacitance.  The  load  is  a 
narrowband  model  of  the  Input  impedance  of  a  230  p  wide 
FET  used  In  the  3.7  to  4.2  GHz  range.  Note  chat  Ideal 
responses  were  obtained  by  neutralising  the  load  reac¬ 
tances  with  negative  components  and  using  an  ideal 
transformer.  Figures  3  and  4  are  plots  of  third  order 
Chabyshev  and  Butterworth  bandpass  filters  with  a 
normal  and  2Z  high  FET  capacitance.  These  plots  are 
only  from  3.73  to  4  GHz  for  clarity,  aad  have  the 
normal  capacitance,  .3  pF,  as  the  circled  points  and 
the  high  value  as  trlaogular  points.  These  filters 
were  designed  for  peak  reflections  less  than  .26  from 
3.7  to  4.2  ®a. 

The  problem  with  the  Chabyshev  filter  Is  not  only 
that  the  reactance  variation  causes  a  change  in  Sn(f) 
that  le  not  quite  tangent  to  Sij(f)  at  the  ripple 
peaks,  but  that  tbs  change  In  S}](f)  due  to  reactive 
load  variations  becomes  tangent  to  the  curve  of  Sn(f) 
at  higher  frequencies.  This  means  that  the  rapid  sag- 
nltude  change  in  Su(f)  causes  reactive  variations  In 
In  the  load  to  cause  near  radial  variations  In  Sn(f) 
and  thus  a  lower  tolerance.  Chabyshev  responses  are 
only  optimal  when  passband  and  stopband  responses  are 
considered.  As  can  be  seen  from  Figures  3  and  4,  the 
al snatched  Butterworth  response  has  a  greater  tolerance 
In  the  aldband  even  though  Its  sasllast  reflection  is 
equal  to  the  peek  reflection  of  the  Chabyshev.  The 
Butterworth  filter's  tolerance  Is  reduced  at  the  band 


Figure  3.  Plot  of  Chebyahev  filter  response  (circles) 
and  response  with  load  capacitance .  Cf ,  In¬ 
creased  by  2Z  (triangles).  Also  shown  is 
the  asxlaua  allowable  reflection,  C*,*  * 
.333.  Only  3.73  to  4  GHz  of  the  3.7  to  4.2 
GHz  response  is  shown  for  clarity.  Note  re¬ 
sponse  variation  to  load  capacitance  becoaea 
tangent  to  the  response  curve  as  frequency 
is  Increased  (dashed  lines).  See  circuit  In 
Figure  2. 


-JSO 


Figure  4.  Plot  of  Butterworth  filter  response  (cir¬ 
cles)  and  response  with  load  capacitance. 
Cf,  increased  by  2X  (triangles).  Also  shown 
is  the  aaxlmua  allowable  reflection,  rmax  - 
.333.  Only  3.75  to  4  GHz  of  the  3.7  to  4.2 
GHz  response  Is  shown  for  clarity.  Note  re¬ 
sponse  variation  to  load  capacitance  becomes 
tangent  to  response  curve  in  the  midband 
(dashed  lines).  See  circuit  In  Figure  2. 


edges  because  variations  In  the  load  reactance  cause 
changes  in  Sn(f)  which  begin  to  be  normal  to  the  var¬ 
iation  In  Su(f)  with  fraquency.  Since  the  Butterworth 
filter  has  near  ideal  tolerance  properties  in  its  aid- 
band,  modifying  Its  band  edge  response  will  lead 
toward  optlaally  load  tolerant  filters. 

The  subject  of  losd  tolerant  filters  In  lossless 
reciprocal  structures  demands  selecting  filter  proper¬ 
ties  that  yield  the  highest  tolerance  to  varlatons  In 
reactive  or  resistive  variations  io  the  load.  Recipro¬ 
cal  structures  aust  rely  on  loss  or  angular  relations 
for  tolerant  design,  as  seen  vis  the  821(f)  812(f) 
terms  In  equations  (1)  and  (2).  The  Si2(f)  tern  In 
equation  (1)  shows  how  circulators  and  Isolators  allow 
tolerance  to  losd  variations  In  nonreciprocal  networks. 
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TABLE  I 


BANDPASS  HATCHING  FILTER  TOLERANCES 


Allowed  Element  Tolerances  (Z) 
for  Tmx  -  .333  «nd  BH  *  3. 7-4. 2  GH* 


Circuit  Parameters 


(aea  Figure  2) 


Device  Circuit 


Filter 

Q 

P  MX 

Cf 

1 

2 

Chsbyshev 

6.23 

.278 

1.2Z 

>.5Z  L2 

>.5X  C2 

Buccervorch 

4 

.28  (.25  midband) 

.81 

>1Z  C3 

>11  L3 

Butcerworth 

4 

.21  (.164  mldband) 

2.1Z 

>3Z  Ci 

>4Z  L2 

Buccervorch 

.25 

.278 

2.2Z 

>6X  n 

>8X  Lf 

Ideal  (eqcn  3) 

0 

0.0 

3.9Z 

TABLE  II 

MATCHING  FILTER  TOLERANCES 

Allowed  Element  Tolerances  (X) 

Broadband  PET  Amplifier  (1-9  GHz) 

for  Tgex  -  .333  and  S2i>1.75 


Filter 


Hatching  Filter 
Ct(pF)  L2(nH) 


Device 

Cf 


Circuit 


Tuned 

Flattest 


.37  .92  >.5Z 

.2  .7  >6Z 


>.5X  L2  >1Z  g. 

>11Z  g„  >I3Z  Rre 


A  second  point  of  clarification  concerns  the  straight 
lines  drawn  in  Figure  1.  Of  course.  the  bilinear 
property  of  networks  causes  all  variations  to  have  a 
finite  curvature.  This  is  most  familiar  to  chose  using 
the  MAP  function  of  COMPACT.  The  result  of  this  Is  to 
cause  Che  calculations  based  on  the  scraighc  lines  of 
Figure  l  to  be  approximate. 

Various  bandpass  networks  represented  by  Figure  2 
were  tested  for  their  tolerance  characteristics  over 
che  3.7  co  A. 2  GHz  band  with  Tut  -  .333.  Table  I  pro¬ 
vides  information  not  only  about  load  variations,  but 
also  about  tolerances  co  filter  components  and  che 
results  of  equation  (3).  The  tolerance  information  Is 
given  for  a  maximum  Sj2(f)  -  .333  with  only  one  compo¬ 
nent  being  varied.  The  additional  tolerances  In  Table 
I  are  for  che  first  and  second  most  sensitive  network 
elements,  as  labeled,  and  their  actual  tolerance  is 
lass  chan  che  next  higher  Integer. 

Table  II  gives  che  coleranca  results  for  a  broad¬ 
band  FRT  feedback  amplifier.  This  amplifier  used 
a  730  u  wide  FET  of  Cf  •  .6pF  and  gg  ■  .063  mhos  with 
an  LC  match. ng  filter  (C[ ,  l2)  and  a  feedback  resistor 
<*FB>  of  180  ohms.  The  FET  modal  Included  source  in¬ 
ductance  and  other  pertinent  components,  fl*o  different 
different  matching  filters  were  used  with  chit  ampli¬ 
fier,  and  have  responses  shown  In  Figure  3.  The  first 
design  was  'tuned’  co  provide  zero  reflection  near 
7  GHz.  and  results  In  a  familiar  humped  response.  The 
second  response  was  designed  co  give  a  finite  but 
flatter  VSWR  across  che  entire  band.  This  lectsr 
design  Is  less  straightforward,  but  results  In  a  much 
greater  tolerance  to  both  che  FET  and  the  circuit 
elements. 
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Figure  3.  Input  VSVR  of  1-9  GHz  FET  feedback  ampli¬ 
fier.  Solid  line  Is  response  with  a  tuned 
matching  filter,  -  .248.  Dashed  line 

is  with  a  flat  matching  filter,  puy.  *  .2. 

CONCLUSION 

In  addition  to  the  presentation  of  several  equa¬ 
tions  useful  for  designing  more  colersoc  circuits.  It 
has  bean  shown  that  the  peaks  and  rapid  magnitude 
changes  of  a  Chsbyshev  filter  create  a  low  tolerance 
response.  The  use  of  nacchlng  filters  which  have  a 
flac  mismatch  across  most  of  che  band  was  shown  co  pro¬ 
vide  a  large  improvement  in  device  and  circuit  tol¬ 
erances  In  all  but  the  large  reflection  with  high  0 
cases.  High  coleranca  reciprocal  networks  are  achieved 
by  designing  che  circuit  response  so  element  variations 
causa  response  variations  which  are  perpendicular  to 
che  response  vector.  The  response  should  also  be  as 
flae  as  possible.  The  synthesis  of  maximally  tolerant 
filters  will  aid  che  production  of  microwave  circuits 
by  simplifying  the  cutting  procedure. 
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ABSTRACT 

A  study  of  Ion- Implant ad  MESFET  performance  as 
a  function  of  the  Implantation  energy  and  fluency 
and  including  the  effects  of  deep-laval  trap 
concentrations  in  the  substrate  has  been 
conducted.  Carrier  concentrations  as  a  function 
of  depth  are  determined  through  the  use  of  LSS 
theory  and  a  profiling  modal.  An  analytic  device 
model,  which  computes  both  DC  and  RF 
characteristics.  Is  then  employed  to  predict 
MESFET  perf ormsnces .  The  study  Includes  the 
effects  of  depth  dependent  transport  properties 
mnA  has  indicated  a  number  of  design  rules  for  the 
fabrication  of  optimized  ion-implanted  devices. 


INTRODUCTION 

The  performance  of  FETs  fabricated  by  ion- 
implantation  depends  greatly  on  carrier 
concentrations  and  velocity-field  characteristics 
as  a  function  of  depth  into  the  active  device 
layers.  The  presence  of  deep-level  traps  in  the 
semiconductor  contributes  to  the  complexity  of 
problems  associated  with  the  characterisation  of 
ion- implanted  devices. 

This  work  is  targeted  mainly  at  determining 
design  rules  for  high-frequency  ion-implanted 
MESFETs.  The  study  utilises  theoretical  models  of 
both  material  properties  end  device 
characteristics.  This  information  is  combined 
with  experimental  material  characterisation  to 
provide  improved  quantitative  accuracy  of  tha 
model. 

Initially,  the  determination  of  typical  ranges 
for  concentrations  of  deep-level  traps,  N ^x) 
(normally  due  to  chromium  in  the  GaAa  substrate) 
was  made  from  a  novel  measurement  technique  using 
both  differential  capacitance  and  conductance  DLTS 
data.  Shallow-level  donor  concentrations,  Ns(x), 
were  then  determined  from  LSS  theory  as  a  function 
of  implantation  parameters.  These  concentrations 
as  a  function  of  depth  wars  than  used  in  a 
profiling  model  to  determine  the  free-carrler 
profile,  n(x) ,  for  the  material.  Carrier 
transport  properties  were  also  determined  from  the 


•This  work  was  supported  by  Rockwell  International 
Electronics  Research  Center  and  tha  U.  S.  Army 
Research  Office,  Research  Triangle  Park,  NC  on 
contract  DAAS29-80-K-4080. 


trap  and  donor  profiles.  This  was  dona  through 
the  use  of  Monts  Carlo  particle  simulations  and  a 
nodal  to  account  for  the  effects  of  compensation 
in  tha  semiconductor - 

An  analytic  model  which  utilises  all  of  thla 
information  is  then  used  to  assess  performance 
potential.  The  distinction  between  free-carrler 
and  donor  profiles,  the  effect  of  deep- level 
traps,  and  the  depth  dependence  of  mobility  and 
velocity  are  shown  to  be  important  considerations 
which  have  been  Ignored  in  previous  models.  The 
DC  characteristics  and  smell-signal  S-parametara 
along  with  figures  of  merit  are  computed  by  the 
model. 

MATERIAL  CHARACTERIZATION 

Figure  1  shows  free-carrler,  background  donor, 
and  deep-level  trap  profiles  typical  of  ion- 
implanted  semiconductor  material.  Deep- level  trap 
centers  and  free-carrler  diffusion  from  highly 
doped  to  lowly  doped  regions  will  cause  tha  free- 
carrler  profile  to  differ  from  that  of  the  ionised 


Figure  1  Typical  concantrelon  profiles  for  ion- 
implanted  material. 

donors.  Our  studies  show  that  as  witch  as  an  order 
of  magnitude  difference  can  exist  between  these 
two  profiles  for  ion-implanted  materials  typically 
used  in  the  fabrication  of  GaAs  MESFETs.  The 
deep-level  sites  will  also  have  a  degrading  effect 
on  carrier  transport  properties. 

Determination  of  the  various  profiles  preeanc 
in  a  sample  is  s  difficult  problem.  Traditional 
C-V  analysis  measures  the  free-carrler 
concentration  if  traps  ere  not  present  (1).  When 
traps  era  present,  however,  knowledge  of  the  trap 
concentration  profile  is  naeeessary  to  extract  the 
free-carrler  information  12).  Likewise,  XTS  data 
can  give  trap  concentration  Information  if  the 
frae-carrler  profile  is  known.  For  this  work. 
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that*  two  oat nod •  were  combined  to  determine  the 
fraa-carrlar,  shallow* Laval  doaor  and  deep-level 
trap  coocaocratloaa  for  a  typical  torn- implanted 
davlca.  Tha  datalls  of  this  technique  ara 
disc  us  sad  alaawhara  (3).  Needed  for  tha  technique 
Is  tha  low-flald  aohility  of  tha  aatarlal  as  a 
function  of  doaor  daaslcy  and  background 
compensation.  This  la  obtaload  using  tha 
thaoraclcal  rasults  of  Ualukiewlcs  it  al.  [4]  la 
conjunction  with  Monta  Carlo  velocity-flald 
pradlctloaa.  Our  own  ttoata  Carlo  rasulta  (3]  wara 
uaad  to  determine  tha  aohility  aa  a  function  of 
background  doaor  daaslty  with  no  trapa  praaant. 
and  tha  Ualufclewlcx  values  wara  chan  aoraalixad  to 
cha  Nonce  Carlo  nuabars.  Tha  noraallxad  data  was 
finally  curva  fit  to  obtain  an  aapirlcal 
axpraaaion  for  aobllity  aa  a  function  of 
background  donor  daaslcy  and  coapansation  ratio. 
Tha  rasulclng  expression  Is 


ko  -  — — -  •  (1-  «)b  (» 

I  +(IaOn)  c 
Ho 

whera 

8380  (ea2/v*sec) , 

No  -  23.2353, 
c  *  23.0, 

8  *  MT/Nd. 

A(log  NdJZ-  B(  log  +  C 


b  « 

0.114992 

with 

A  -  0.025, 

8  «  0.817278, 

C  -  6.252838 


for  %  >  1021  (a-3), 
for  Sj,  <  1021  (a'3). 


and  whara  N_  Is  glvaa  In  (m-3).  Expression  (1)  la 
ploctad  against  cha  aoraalixad  Ualukiewlcs  valuas 


rigura  2  Low-flald  aobllity  as  a  function  of 
doping  and  daap-laval  coapansation.  Tha 
solid  Unas  ara  coaputad  froa  aquation 
(1)  for  8*0.0,  0.2,  0.*,  0.6,  0.8.  Tha 
data  points  ara  froa  cha  aoraalixad 
chsoratlcal  coaputaclona  (4]. 

la  Figure  2.  As  can  ba  scan  froa  tha  plot,  tha 
agraaasnt  la  quits  good. 


k  ona  alcron  gaca  langch  DUTR  along  with  a 
dlffaraacial  capaclcanca  pat  cam  wara  fabrlcatad 
on  Silicon  laplantad  Cr-dopad  SsAa  substrata. 
Cooductanca  OLTS  and  C-V  nans ur an ants  wara 
parfocaad  on  thasa  devices  and  a  dan  inane  daap- 
laval  trap  ataca  war  ldanelfisd  0.736  aV  balow  tha 
conduction  band.  Tha  nsasuransnes  wara  uaad  as 
Input  data  to  a  conputar  simulation  which  conputaa 
tha  daslrod  profllas  aa  nanelonad  abova. 

Tha  final  r a suit lag  frsa-carrlsr ,  shallow-laval 
donor,  and  daap-laval  crap  conesntratloos  aa  a 
function  of  dapeh  Into  cha  aatarlal  ara  shown  In 


Flgura  3  Rasulclng  concancratlon  profllas  and 
low-flald  aobllity  profile  for  one 
device . 

Flgura  3.  Notice  that  deap  Into  tha  channel  there 
la  sons  scatter  of  tha  shallow-laval  donor  data. 
This  bagina  to  occur  whan  cha  crap  concancratlon 
and  tha  shallow-laval  donor  concentration  ara  of 
tha  same  order  of  magnitude.  The  uncertainties  in 
the  exact  shallow-level  concentration  at  this 
depth  Into  tha  channel  ara  noc  critical  co  cha 
profile  predictions.  This  la  true  since  cha 
aagnlcude  of  all  tha  profllas  of  Interest  ara 
snail  at  this  dapeh  when  coopered  to  their 
magnitudes  near  tha  Implantation  peak. 

The  low-flald  mobility  profile  obtained  froa 
this  analysis  la  also  shown  In  Figure  3.  Tha 
curva  can  be  compared  co  Che  results  of  Das  and 
Rim  [6],  end  is  In  good  qualitative  agreement. 

THS  DggXCj  MODEL 


The  device  modal  used  here  la  a  one-dimensional 
model  with  a  small-signal  analysis.  Onm- 
dlmanslonal  models  offer  a  number  of  advantages 
for  a  study  of  this  type  over  mqra  elaborate  two- 
dimensional  models.  For  example,  cha  fora  of  tha 
rasulta  obtained  froa  one- dimensional  analysis  la 
more  useful  In  terms  of  many  device  optimisation 
and  design  problems.  Also,  Information  needed  to 
examine  eircuit/dsvlce  Interface  phenomena  can 
easily  ba  obtaload  froa  one-dlaaasional  rasults. 
Of  critical  Importance  hare,  however,  is  the  cost 
difference  between  one-  and  two-dimensional 
simulations..  In  thasa  studies,  over  one-hundred 
different  device  geometries  wars  examined  at  can 
to  twenty  different  bias  levels.  To  obtain  this 
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Intonation  Iron  two-dimensional  simulation!  Mould 
have  boon  at  laaat  an  otdar  of  magnitude  aora 
costly  in  cans  of  both  tlaa  and  money. 

The  aodal  asaunas  that  tha  alactron  cranaport 
propartlas  of  a  aacarlal  can  bo  slaulatad  by  a 
two-piece  valoclty-fiald  ralatlonahip.  Tha  two- 
piaca  approximation  Is  dafinad  from  a  thaoratleal 
valoclty-fiald  charactaristic  datanlnad  by  Itonta 
Carlo  techniques.  For  alactric  flalds  lass  chan 
an  appropriata  saturation  fiald,  Em  tha  alactron 
velocity  is  dascrlbad  by  s  linear  axprasaion, 


For  electric  fields  above  Eg  tha  electrons  nova  at 
a  constant,  aaxlmun  velocity,  v«.  Us In*  this 
Information,  a  small-signal  equivalent  circuit  as 
shown  in  Figure  4  can  be  datanlnad  aiwi  analysed 
to  obtain  RF  performance  prodletlons,  including 
power  gains  as  a  function  of  frequency. 

Low-field  mobility  as  a  function  of  depth  for 
tha  model  la  obtained  directly  from  aquation  (1). 
Determining  the  appropriata  maximum  velocity  for 
the  material  is  more  involved,  however. 

Tha  importance  of  developing  a  systematic, 
well-justified  technique  for  determining  maximum 
velocity,  v„,  was  discussed  in  aarliar  work  {7,8] 
and  in  this  study,  such  a  mac  hod  haa  been 
developed.  Tha  method  involves  numerical 
determination  of  the  carrier  transit  tlma  under 
the  gate  using  an  exact  velocity-field 
relationship.  The  aaatmMd  field  distribution  for 
this  calculation  la  obtained  from  two-dimensional 
simulation  results.  It  is  than  required  that  tha 
transit  time  calculated  assuming  a  two-piaca 
approximated  velocity-field  relationship  be  equal 
to  that  of  tha  exact  analysis.  This  method  yields 
a  v„  for  Case  (doped  to  a  level  lD  -  1017aa“3  )  in 
excellent  agreement  with  the  value  obtained  by 
Fucal  at  al.  {•].  The  resulting  data  was  curve 
fit  and  expressed  as 


v,  *  v0  -  A  log  [(1-y)2  +  *y] 


y  •  (Unix)/*-)2'5, 

■d-  1.5  x  10«, 

A  -  0.0282, 

■  -  0.4, 
v^  1.40. 

Equation  (3)  gives  va  in  105  a/see  when  Eq(x)  Is 
expressed  in  a"3 .  The  expression  is  also  assumed 
to  have  the  sane  dependence  on  compensation  ratio, 
6,  as  equation  (1).  Thus,  the  factor  (1  •  8)b  is 
multiplied  with  equation  (1)  to  obtain  ^  in  the 
presence  of  traps.  The  exponent,  b,  is  defined  in 
equation  (1). 

Equations  (1)  sad  (3)  in  conjunction  with 
knowledge  of  the  three  profiles,  n(x),  Np(x)  sad 
Nf(x),  allow  for  the  derivation  of  a  device  nodal 
which  Includes  the  effects  of  varying  transport 
properties  as  a  function  of  epl-dspth. 

The  current-voltage  predictions  of  the  aodal 
for  a  one  aicron  ion- implanted  device  are  compared 


Figure  4  The  equivalent  circuit  for  an  FET  used 
la  the  analysis. 


Mam  VOLTLQS  I  VOLTS) 

Figure  5  Comparison  of  model  predicted  and 
measured  1-V  characteristics  for  s  one 
aicron  lon-laplanted  device. 

in  Figure  5  with  the  measured  I-V  curves.  The 
agreement  is  excellent.  It  should  be  noted  that 
without  the  Inclusion  of  the  effects  of  traps  on 
carrier  transport,  this  agreement  could  not  be 
obtained.  The  traps  have  a  tendency  to  'soften' 
the  pinch-off  characteristics  of  the  device  (  le. 
whan  traps  are  Included  in  the  simulation,  the 
slope  dip/ dV g  is  not  as  great  near  pinch-off). 
For  all  of  the  devices  studied  in  this  work,  this 
'softening'  effect  was  required  to  obtain  best 
agreaaaat  with  measured  I-V  characteristics. 

EESPLTS  AMD  COUCLPSIOKg 

The  results  that  follow  ware  obtained  following 
the  modeling  steps  outlined  in  the  Introduction. 
Three  parameters  were  varied  Independently.  They 
are  1)  implantation  energy,  £  ,  2)  peak  doping 

density,  gj^,.  and  3)  trap  concentration,  NT(x). 
The  implant  species  was  assumed  to  be  Si  in  GaAs 
and  the  activation  was  assumed  to  be  100X  for  all 
devices.  The  trap  concentration  was  assumed  to  be 
constant  as  a  function  of  depth  for  these  studies. 
Mote  that  tha  peak  doping  density  can  be  converted 
to  e  corresponding  ion  fluency  through  the  simple 
relationship 


♦  ■  Op  Syj 


Op  is  tha  standard  deviation  of  the 
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orojectad  »*<»*•• 

For  cho  implant  anargy  and  peak  doping  studies, 
the  crap  lava!  hi  left  eoaacaac  ac  Jfy-2.0  x  10 15 
ca"3 ,  This  iiabac  was  chosen  to  be  in  ganaral 
agreement  with  cha  raaulca  shown  la  Flgura  3.  Tha 
iaplanc  energy  waa  varied  from  50  Co  150  kaV 
whila  Cha  paak  doping  cook  values  bacwaan  8  t  10^ 
and  4  x  1017  aa~3. 

Figuraa  6a  and  6b  illuacraca  cha  effects  chasa 
paraaacara  hava  on  cha  gain- band tfldch  produce,  f  t. 


so  too  (SO 

MflUT  SNOSHV.  {Kill 


Flgura  6a  PradlcCad  gain-bandwidch  produce,  fy , 
»a.  Implantation  anargy,  £  .  Tha  paak 
doping,  H,„ .  la  uaad  aa  an  Independent 
paraaacar. 


so  too  iso 

IMPLANT  C NOSOV.  (Sotfl 


Flgura  6b  PradlcCad  aaxlaia  frequency  of 
oaclllacion,  f  aajr  vs.  lnplancaclon 
anargy,  £  .  Tha  paak  doping,  Noax> 
uaad  aa  aa  indapandanc  paraaacar. 

and  eha  aaxlaun  fraquancy  of  oaclllacion,  f .... 
Tha  valua  for  fy  la  eonpucad  froa  flrae  ordar 
conaidaraclona  co  ba  given  by 

■  *m/2»Cgs.  (7) 

Tha  quandey  f  mtx  la  obealnad  by  noting  eha 

fraquancy  ac  which  Mason's  unilateral  gala  (aa 

predicted  by  cha  nodal)  goas  co  unity.  Tha 

raaulca  claarly  Indicate  tha  superiority  of  low- 
anargy  laplants  for  high-frequency  oparatlon. 

dll  of  tha  devices  conaldarad  la  coaplllag 
Figuraa  6a  and  6b  wore  conparad  at  a  blaa  of  XjflO 
sa.  Tha  dashed  line  falling  off  rapidly  at  e!n 
low- anargy  sod  of  cha  curves  Indicates  that  for 


anarglas  lower  chan  this,  Toss<10  t>A. 

For  the  crap  study,  an  Iaplanc  anargy  of  70  kaV 
and  a  paak  doping  density  of  2  x  10i7  c£ 3  was 
assuaad.  Tha  background  crap  concancracion  was 
varied  froa  0  co  1016  as-3. 

Figure  7  Illustrates  Che  effects  Craps  hava  on 
cha  zero  gaea  blaa  currant,  IDSS>  and  cha  pinch- 
off  potential  •  *  oo'  Flgura  8  shows  cha  plnch-off 
"softening*  efface  aeneioned  In  the  previous 


HUP  CONCSNTNATNM  tt»*l 


Figure  7  PradlcCad  affects  of  varying  crap 
concancrstlona  on  DC  characcarlstics. 


NORUALIZSO  OATS  VOL  TAOS  H,/V„) 


Flgura  8  Mo  nail  sad  drain  curranc  predictions  vs. 

nornalissd  gate  voltage  with  and  without 
traps. 

section  for  one  particular  devleo. 

Deep- lav el  traps  also  hava  sane  effect  oo  the 
IP  performance  of  tha  device  as  Illustrated  In 
Figuraa  9  and  10.  Figure  9  shown  claarly  tha 
decrease  of  f  T«od  fM,  associated  with  Increasing 
trap  concentrations . 

In  Flgura  10  aote  chat  for  low  bias  currents  fy 
Increases  whan  few  traps  are  present  while  It 
decreases  for  higher  trap  concantratlons.  This  is 
saslly  sxplsined  In  tarns  of  the  degrading  effects 
traps  hove  on  nobility  and  velocity.  As  the  gate 
bias  restricts  current  flow,  a  larger  fraction  of 
tha  carrlors  are  forced  deeper  into  the  channel. 
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Figure  9  Predicted  gain- bead  width  produet,  fj,  end 
maximum  frequency  of  oec illation,  f 
ve.  background  trap  concentration. 


Figure  10  Predicted  gain-bandwidth  product,  fj ,  ea 

a  function  of  normalised  drain  current. 

The  background  trap  concentration  la 

uaad  aa  a  variable  paraaeter. 

This  correaponda  to  the  sore  lightly  doped  regiona 
of  the  device.  If  the  compensation  ratio  la 
fairly  email  (  la.  few  trapa)  then  f roe  equatlona 
(1)  and  (3)  the  tranaport  propertlaa  are  auperlor, 
and  ff  increaaaa.  If,  on  the  other  hand,  the  trap 
level  la  on  the  aaee  order  of  eagnituda  aa  the 
a hallow-lev el  donor  concentration,  then  the 
coepenaatlon  ratio  approochea  one.  Thla 
correaponda  to  axtrene  degradation  of  nobility  and 
velocity  and,  thua,  foreaa  fy  to  decreaae. 

In  conclualon,  deep  level  trapa  in  ion- 
iaplanted  devicea  degrade  carrier  tranaport 
propertlaa  in  the  aanlconductor  eatarlal.  The 
degradation  la  more  aavare  near  the  tall  of  the 
implant  profile  than  near  the  peak.  Thua,  tha 
tranaport  propertlaa  of  the  device  will  be  depth 
—  or  blae  —  dependent.  A  modeling  technique 
which  account a  for  thla  dependence  hae  bean  uaed 
to  etudy  device  propertlaa  aa  a  function  of 
fabrication  peraaetera  and  deep-level  trap 
concentratlona.  The  reaulta  indicate  that  low 
energy  iaplanta  abould  poeaae  auperlor  high- 


frequency  propertlaa  and  that  lowering  tha  trap 
lavala  in  the  material  abould  improve  device 
performance. 
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ELECTRON  CONCENTRATION  AND  ALLOY 
COMPOSITION  DEPENDENCE  OF  HALL 
FACTOR  IN  Ga,ln1-,A*rP1.r 


Indexing  terms:  Semiconductor  inices  and  materials.  Hall 
•fftcu  lli-V  semiconductors 

]*•  Hall  factor  is  calculated  by  the  iterative  method  over  the 
whole  compost  boo  rup  at  Ga>,.,As  P.  ,  Uttice- 
^  ei*tr00  concentrations  used  are 
1  *  10  cm  sod  1  x  10 17  cm*’,  with  carrier  ootnpensa- 
ooo  ratios  (SB  +  NA)/n  1.  2.  sad  J,  which  ate  commonly 
obeerwd  m  poem  GalnAsP  alloys.  The  Hall  factor  decreases 
with  laareaaing  alloy  content  y  far  relatively  low  electron  coo- 
QgCTtxHU,  but  this  is  not  the  case  in  compensated  alloys 
because  of  varying  predominance  of  scattering  mechanisms 
with  composition  and  ionised  impurity  concentration. 

There  is  rapidly  growing  interest  in  the  transport  properties  of 
GalnAsP  alloys  lattice- matched  to  InP.1"1®  Improved  FET 
performance  is  expected  for  GalnAsP  because  of  its  high,, 
electron  motoUily1-*-*-1®  and  higher  electron  drift 
velocity  •  than  those  of  GaAs,  In  the  analyses  of  these 
transport  properties  the  measurement  of  low-field  mobilities  is 
one  of  the  essential  evaluations.  The  mobility  is  usually  derived 
from  a  combination  of  the  resistivity  and  the  Hall -effect 
measurement,  and  the  thus  derived  mobility  is  the  Hall  mobi¬ 
lity.  On  the  other  hand,  the  slope  of  the  velocity/field  curve  at  a 
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low  electric  iield  strength,  which  is  equal  to  the  drift  mobility, 
is  very  often  normalised  by,  or  compared  with,  the  measured 
Hall  mobility.  However,  the  Hall  mobility  ttH  is  not  generally 
equal  to  the  drift  mobility.  The  ratio  is  the  Hall  factor 
rH.  This  factor  also  correlates  the  carrier  concentration  n  and 
the  Hall  coefficient  R„  through  the  equation  n  —  rK/eRH 
where  t  is  the  electronic  charge.  In  the  limit  of  high  magnetic 
fields,  and/or  carrier  degeneracy,  r„  approaches  unity,  m»  *■  Po 
and  the  carrier  concentration  can  be  derived  directly  from  a 
measured  Hall  coefficient  However,  it  is  not  practical  to  attain 
the  high  magnetic  field  limit,  and  so  for  many  applications 
nondegenerate  semiconductors  are  of  interest 

There  are  a  few  experimer  ts  of  the  determination  of  the  Hall 
factor  r„  in  GalnAsP  from  comparisons  of  Hall  mobilities  at 
low  and  high  magnetic  field. 5,11  Although  measured  values  are 
close  to  unity,  one  paper’  reported  that  it  was  equal  to  or  less 
than  1,  and  others’2  observed  it  equal  to  or  greater  than  1. 
Since  the  Hall  factor  has  a  dependence  on  scattering  mechan¬ 
isms,  its  value  varies  with  the  carrier  concentration  and  the 
composition  of  GalnAsP  alloys.  There  are  no  reports  of  the 
Hall  factor  covering  the  whole  alloy  composition  of  GalnAsP 
which  show  the  carrier  concentration  dependence. 

In  this  letter,  we  show  the  calculated  Hall  factors  in  Galn¬ 
AsP  over  the  whole  alloy  composition  lattice-matched  to  InP 
in  a  moderate  electron  concentration  range  between  1  x  10“ 
cm" 1  and  l  x  1017  cm'1  at  300  K. 

The  calculation  method  is  an  iterative  technique  in  which 
the  Boltzmann  equation  is  solved  by  numerical  iteration.1 3-1  ’ 
Both  the  drift  mobility  and  the  Hall  mobility  are  calculated  to 
any  arbitrary  degree  of  convergence  without  using  the  relaxa¬ 
tion  time  approximation,1*  Matthiessen’s  rule,  or  complicated 
mathematical  expressions.17'1*  All  the  calculations  in  this 
letter  are  carried  out  within  0-5%  error.  The  allowed  scattering 
mechanisms  are  polar-optical  phonon  scattering,  acoustic 
phonon  scattering,  piezoelectric  scattering,  ionised-impurity 
scattering  and  alloy  scattering.  The  quaternary  alloy  scattering 
rate  proposed  by  Littlejohn  er  a/.20  is  used  with  the  Phillips 
electronegativity  difference  for  the  scattering  potentials.20-21 
The  material  parameters  for  the  quaternary  alloy  were  eva¬ 
luated  by  the  linear  interpolation  between  InP  and 
Ino.5jGao.47As.  Some  of  the  parameters  of  InP  and 


n*.  I  Alloy  amfotUUm  dependence  of  Hall  factor  In 
Ga,lalm,As,Pl.tal  100  K 

Electron  coaesatndoni  are  1  x  10“  cm*’  (solid  liaet)  and 
I  x  1017  cm*'  (dashed  tinea).  Carrier  compeasatioe  ratio* 
(Af»  +  Nty*  are  I,  2  and  5,  impec lively 
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Ino.5jGao.47As  used  in  the  calculations  are  listed  in  Table  1. 
The  electron  distribution  function  is  the  Fermi- Dirac  function 
because  the  electron  distribution  in  Ino.5jGao.47As  may 
degenerate  at  a  relatively  low  electron  concentration 
(it -2  x  JO17  cm* 2  at  300  K). 

Figs,  la-c  show  the  calculated  Hall  factor  as  a  function  of 
the  alloy  composition  y  in  Ga.l^-jAs^P,  _r  Carrier  compen¬ 
sation  ratios  {N„  +  1.  3  and  5,  respectively.  Here, 

Afo,  Na  and  n  are  the  donor  estxxn  (ration,  the  acceptor 
concentration,  and  the  free-electron  concentration,  respec¬ 
tively.  In  the  model  both  the  donor  and  acceptor  are  assumed 
to  be  fully  ionised,  so  that  n  is  equal  to  (ND  -  SA)  and 
{No  +  Na)  is  the  total  ionised  impurity  concentration.  As  a 
general  tendency  in  a  purer  alloy  semiconductor,  the  Hall 
factor  decreases  with  increasing  alloy  composition  y.  This 
comes  from  the  decreasing  phonon  energy  with  composition  y, 
as  can  be  seen  from  Table  1.  The  value  of  the  Hall  factor  for 


Table  l  SOME  OF  THE  MATERIAL  PARAMETERS 
OF  InP  AND  Ino.5jGao.47As  USED  IN 
THE  CALCULATIONSt 


InP  I  n0. 5  3  Ga,>.  As 

Effective  mass  ratio  m*/m„ 

008 

0041 

Static  dielectric  constant  e. 

12-35 

13-77 

Optical  dielectric  constant  e«, 

9-52 

11-38 

Optical  phonon  Debye  temperature  (K) 

504 

396 

Energy  gap  at  0  K  (eV) 

1-42 

0812 

Acoustic  deformation  potential  (eV) 

6-5 

589 

Piezoelectric  constant  (C/m2) 

0-035 

0099 

t  Linear  interpolations  between  the  two  wen  used  for  material 
parameters  of  GalnAsP 

The  effective  masses  an  found  in  Reference  24.  Dielectric  constants, 
energy  gap  at  0  K.  acoustic  deformation  potential  and  piezoelectric 
constant  of  InP  are  found  in  Reference  13.  Optical  phonon  «««nw«  an 
reported  by  Reference  23.  Other  parameters  of  In«.,,Ga*47As  an 
linearly  interpolated  between  GaA$  and  InAs 


polar-optical  phonon  scattering  is  a  maximum  at  a  tempera¬ 
ture  close  to  the  Debye  temperature  of  the  polar  optical 
phonon.14'22  As  the  ambient  temperature  increases  r„  falls. 
For  GalnAsP  the  Debye  temperature  of  the  polar  optical 
phonoo  varies  from  504  K.  to  396  K  with  increasing  y,  so  that 
at  room  temperature  the  value  r„  decreases  with  increasing  y. 
However,  in  a  combination  of  several  scattering  mechanisms 
this  general  tendency  is  modified.  Besides  polar  optical  phonon 
scattering,  alloy  scattering  and  ionised-impurity  scattering  are 
the  influential  scattering  mechanisms  in  a  moderately  doped 
alloy  at  room  temperature.  The  combination  of  alloy  scatter¬ 
ing  flattens  the  temperature  dependence  of  the  Hall  factor,  but 
the  value  at  300  K  is  little  changed.21  On  the  other  hand,  the 
combination  of  ionised  impurity  scattering  and  polar -optical 
phonon  scattering  tends  to  make  the  total  scattering  rate  less 
momentum-dependent,  so  that  the  value  of  the  Hall  factor 
decreases,  When  the  scattering  rate  is  independent  of  carrier 
momentum,  the  Hail  factor  is  equal  to  1.  Because  of  the  vary¬ 
ing  importance  of  the  scattering  mechanisms  with  the  alloy 
composition  y,  the  Hall  factor  does  not  necessarily  decrease 
monotonically  with  increasing  y  in  a  moderately  doped  alloy. 
In  the  Figure  the  Hall  factor  has  the  lowest  value  1-02  near 
y  -  02  at  b  -  l  x  1017  cm'1  and  R  -  2. 

At  a  higher  carrier  compensation  ratio  the  Hall  factor  in¬ 
creases  because  of  the  predominance  of  ionised-impurity  scat¬ 
tering  which  has  a  Hall  factor  ss  large  as  1-93. 

In  summary,  the  Hall  factor  in  Gi,ln1-,AsrP1-„  depen¬ 
dent  on  electron  concentration  and  alloy  composition,  has 
been  calculated  by  the  iterative  method.  In  a  moderate  electron 
concentration  range  the  Hall  factor  is  greater  than  1  over  the 
whole  composition.  The  results  should  be  uteftil  in  the  calcula¬ 
tion  of  drift  mobility  and  electron  concentration  from  the  re¬ 
sistivity  and  the  Hall  effect  measurements  in  which  the  Hall 
mobility  and  the  Hall  coefficient  are  obtained. 
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ABSTRACT 


A  new  odd  order  impedance  matching-  network  with  reduced  aenaltlvlty  to  active  device  capacitance  variations 
la  presented.  A  synthesis  procedure  for  these  networks  Is  presented  and  experimentally  verified  with  the  con¬ 
struction  of  a  alcrowave  amplifier.  These  networks  are  useful  in  the  development  of  low  coat  alcrowave  Integrated 
circuits  since  they  reduce  the  harmful  effects  of  device  variations. 


Introduction 

In  this  paper  a  general  analytic  method  is  pre¬ 
sented  chat  extends  network  synthesis  to  Include  both 
even  and  odd  order  low-pass  Impedance  matching  fil¬ 
ters.  Previous  network  synthesis  procedures  for  la- 
edance  matching  filters  were  based  on  the  use  of  even 
order  filter  structures.  These  odd  order  networks 
are  Important  since  they  are  relatively  Insensitive 
to  varlatlooe  In  the  capacity  associated  with  a  load 
Impedance  lower  chan  the  generator  Impedance.  The 
lasemeltlvity  of  these  nacworka  to  capacitance  varle- 
tlooe  can  minimise  or  eliminate  the  eoetly  tweaking 
of  clreulte  caused  by  variations  in  the  active  dev- 
Icee.  Two  claaaas  of  polynomials  are  presented  so 
that  bach  broed-and  narrow-band  circuits  any  be  de¬ 
signed. 


Discussion 

The  development  of  even  order  matching  filters 
(11  has  proved  vary  useful  la  amplifier  design.  It 
will  he  shown  In  this  paper  chat  odd  order  filters 
are  unique  In  that  they  provide  an  insensitive  cir¬ 
cuit  element  lnetaad  of  broadening  the  bandwidth  of 
the  Impedance  match. 

The  design  procedure  Is  based  upon  consideration 
of  the  fundamental  relatlooahtp 

l*H<«s>lZ  *  l*2l(w)|Z  -  1.  (1) 

with 


(2) 

«P) 


*2l<»>  ■  L  .  O) 

7TiTp) 

sad 


«  *  CR|-«^)Z  <*> 

where  p  -  s  ♦  Jw.  The  tens  Rj  and  R,  are  the  gene¬ 
rator  and  load  resistances,  respectively.  The  char¬ 
acteristic  function,  D(p)  may  be  defined  by 

K»>  -  (l*bp)M(p2),  (J) 


where  the  parameter  b  is  proportional  to  the  square 
of  the  parasitic  capacity  associated  with  the  device 
co  be  matched.  This  parameter  la  set  to  aero  for 
even  order  networks.  The  even  polynomial,  M(pS),  may 
be  defined  according  to  the  traoafotmncloa  presented 
in  (t|  for  broad-band  matching,  or  It  may  use  the  co¬ 
efficients  of  the  binomial  eapanalon  for  narrow-band 
matching.  During  synthesis,  the  polynomial  Mp1)  de¬ 
termines  the  bandwidth  and  ripple  of  the  Impedance 
match,,  and  appronlmate  formulas  for  b  In  terms  of  the 
load  capacity  nay  be  derived  which  thus  define 
*21<P).  Once  Sj[(p)  is  defined,  the  derivation  of  a 
.network  Is  straightforward. 


Theoretical  Results 

The  relative  Insensitivity  of  the  odd  order  low 
pass  matching  filters  co  varlatlooe  In  a  shunt  capac¬ 
itor  on  Che  low  Impedance  side  of  the  filter  may  be 
demonstrated  both  analytically  and  experimentally. 
For  the  analytic  case,  the  circuit  of  Fig.  la  was  de¬ 
signed  from  a  Cebysev  prototype  with  b-t  and  a  frac¬ 
tional  bandwidth  of  31Z.  The  VSUR  plots  in  Fig.  lb 
were  then  obtained  from  the  Input  of  the  filter  as 
the  capacitor  on  the  low  Impedance  aide  was  varied 
t)0X.  Fig.  h  shows  an  even  order  (s»i)  Cabyeev 
matching  filter  with  the  same  parameters  as  the  cir¬ 
cuit  In  Fig.  la  except  that  b»0.  By  comparing  Figs, 
lb  and  ,2b  one  may  see  that  the  synthesised  filters 
give  nearly  Identical  response  (solid  curves),  but 
chat  the  t30X  variations  In  the  component  nearest  the 
load  causes  much  greater  variations'  In  the  TM  of 
the  even  order  filter.  For  example,  the  midhead  9 SUE 
Is  designed  as  l.l,  but  reaches  1.67S.04  under  varia¬ 
tion*  In  the  even  order  filter  while  ealy  becoming 
1. Sit. OS  under  variations  la  Che  add  order  filter. 
The  odd  order  matching  filter  thus  uses  the  addition¬ 
al  element  far  Insensitivity  Instead  of  bandwidth. 


mantel  Results 

As  amplifier  selsg  a  Hater ala  igf  901  bipolar 
transistor  has  been  destgaed  malag  a  fifth  order 
Cebysev  matching  network  for  Che  Input  circuit,  and 
the  lossy  gala  equalisation  circuit  described  in  (21 
for  the  output  circuit.  Fig.  Si  shows  the  amplifier 
designed  for  operation  from  690  Mis  to  900  Mis  at  a 
bias  currant  of  9  uA  and  9v  T^g.  The  parameters  of 
the  Input  circuit  design  were  b*S,  Eg  -  IS  ohms.  fa  - 
100  Mis  and  A  •  .3119.  This  circuit  design  pthduced 
s  theoretical  and  experimental  gain  of  19.91.R  db 
from  690  to  900  Mis. 


•This  work  was  supported  by  the  0.  S.  Army  under 
cooCtacr  0AAC29-60-E- 0060 . 
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(a) 
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.  1  (a)  Odd  order  (a-3)  Cabysev  low-paoa 
matching  filter  with  a  bandwidth 
parameter.  A,  of  .31. 

(b)  The  Input  VSVK  of  thla  network 
ae  the  1.64F  capacitor,  C-,  la 
varied  ±30Z.  The  solid  daahed 
and  alternating  lloea  are  for  . 

C  -  1.64F,  2.46F  and  .82F 
respectively. 


m  ^ 

Tt-»  T'-y 

(e) 


(b) 


Mg.  2  (e)  Even  order  (n*4)  Cabysev  low-paaa 

■etching  filter  with  a  bandwidth 
parameter.  A,  of  .31.  Thla  la 
the  filter  of  Fig.  1(a)  with  b-0. 
(b)  The  Input  VSWX  of  ehla  network  aa 
the  .331  inductor,  L, ,  la  varlad 
tSOZ.  The  solid,  dashed  and 
alternating  llaee  are  for  L.-.33K. 
.3231  and  .173H  reapectlvely. 


rx 

The  lneeneltlvlty  of  the  odd  order  matching  net* 
work  to  changea  In  the  devlce'e  input  capacity  la  Il¬ 
lustrated  In  Flga.  3b, c  which  coapares  the  odd  order 
■etching  network  with  an  Ideal  transformer.  Fig.  3b 
shows  tha  vswt  of  the  aapllflar  with  an  odd  order 
■etching  network  In  theory  and  experlaent  aa  well  aa 
with  aa  Ideal  transformer  of  1.38:1  turns  ratio.  The 
curves  are  represented  by  solid,  dashed  and  alternat¬ 
ing  lines,  respectively.  While  Fig.  3b  shows  reason¬ 
able  VSWK  response  In  all  cases.  It  la  Fig.  3c  which 
shows  the  need  for  odd  order  Batching  filters.  Fig. 
3c  shows  the  change  In  VSVK  when  a  5  pF  capacitor, 
giving  approxlaately  30Z  sore  Input  capacity,  la 
placed  In  parallel  with  the  transistor's  base-ealtter 
junction.  Tha  change  In  VSVK  for  the  Ideal  trans¬ 
former  notched  aapllflar  la  obviously  greater,  and 
thus  sore  sensitive,  than  the  changea  la  both  the 
theoretical  and  experimental  odd  order  Batching  fil¬ 
ter.  Thla  la  because  the  shunt  capacitor  on  the. low 
lnpedance  side  of  the  filter  does  not  contribute  to 
the  lnpedance  Increase  needed  for  Batching. 


(e) 


Fig.  3  (a)  MRF-901  amplifier  with  distributed  fifth 

order  matching  filter.  Kota  the  100  pF 
blocking  capacitor  and  *0*  long  trans¬ 
mission  line  are  far  the  bias  network 
only. 

(b)  VSV*  Of  MKF-901  amplifier  la  theory 
(solid  lino),  experiment  (daahed  line) 
and  when  the  Input  network  la  replaced 
by  a  1.38:1  Ideal  transformer  (alternat¬ 
ing  11m). 

(c)  Change  la  VSVK  of  KEF-tOl  aapllflar  as  a 
3  pF  capacitor  lo  placed  serosa  its  base- 
emitter  junction.  Sana  11m  equivalences 
an  la  (b). 
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Those  networks  can  also  ba  uaad  Co  lap rove  the 
Insensitivity  of  broadband  F.E.T.  feedback  amplifi¬ 
ers.  Such  an  aapllf ler  was  designed  using  an  Avancek 
AT-8111  F.E.T.  chip  blaaad  aC  approximately  JJ  uA  and 
4v  VDS.  The  circuit  shown  In  Fig.  4a  producer  51.3 
db  gain  from  1  to  10  CHx  and  an  Input  VSWk  shown  by 
tha  solid  curve  in  Fig.  4b.  The  daahed,  alternating 
and  double  alternating  curves  In  Fig.  4b  represent 
variations  of  the  circuit  In  Fig.  4a  with  no  matching 
network',  a  Smith  chart  matching  network  design,  and  a 
fifth  order  Gebysav  matching  filter;  respectively. 
Mote  how  the  circuit  exhibits  a  rising  VSWK  versus 
frequency  with  no  matching  network  (dashed  curve), 
and  a  VSWK  which  goes  to  1:1  and  then  rises  sharply 
above  8.2  CHx  for  the  tuned  charactc riot lea  of  eha 
Saleh  chert  matching  neework  (alternating  curve). 
The  analytic  aaechlng  networks  of  third  and  fifth  or¬ 
der  both  exhibit  smooth  VSVK  curves  oscillating 
around  1.3:1  and  never  exceeding  1.82:1.  The  capaci¬ 
tor.  In  the  matching  network  oeareee  eha  devlea  la  ac¬ 
tually  the  Intrinsic  device  capacity,  and  thus  la  ac¬ 
counted  for  In  the  deelgn  but  not  present  In  the  cir¬ 
cuit  layout.  The  Smith  chart  design  method  Inherent¬ 
ly  accounts  for  the  device  Input  capacley,  and  so 
forms  an  odd  ordar  matching  network  albalc  of  nonop- 
tlaal  design.  The  curves  of  Fig.  4c  show  the  change 
In  VSVK  of  aach  circuit  when  eha  F.E.T.  Input  capaci¬ 
tance  la  raised. by  30X  (.193  pF).  Nora  how  the  cir¬ 
cuit  with  no  aaechlng  network,  exhibits  a  rising 
change  In  VSWK  versus,  frequency,  while  all  the  odd 
-  order  aaechlng  networks  exhibit  a  vary  low  change  In 
VSWK  up  to  aoma  break  frequency.  Thus,  these  match¬ 
ing  networks  trade  lneenslelvlty  at  moat  frequencies 
for  an  Increased  sensitivity  at  vary  high  frequen¬ 
cies.  Tha  Smith  chart  matchlag  netwrok  (alternating 

curve)  ban  tha  lowest  break  frequency  and  highest 
sensitivity  at  10  CHx  as  shown  In  Fig.  4c.  and  so  Is 
not  the  optimal  matching  naework.  Tha  fifth  ordar 
network  (double  alternating  curve)  shows  a  very  low 
VSVK  change  up  eo  9  CHx  which  makes  It  tha  optimum 
natwork  for  Increases  In  capacley.  Tha  poor  perform¬ 
ance  of  tha  fifth  order  natwork  far  reductions  la 
capacity  may  make  this  natwork  leas  than  optimum  In 
soma  situations.  Tha  analytically  designed  third 
ordar  network  gives  an  overall  improvement  over  tha 
Smith  chart  aaechlng  network,  and  no  aaechlng  natwork 
la  all  capacitance  variations.  Thus,  this  method  of 
analytic  aaechlng  network  daslgn  la  shown  to  yield 
Improved  and  relatively  lnaanaltlva  circuit  designs. 


Conclusion 

la  conclusion,  the  theory  of  low  peas  matching 
filters  la  extended  to  Include  odd  order  networks. 
The  aaalytlc  theory  presented  here  should  be  of  par¬ 
ticular  use  In  designing  Internal  aaechlng  networks 
for  bipolar  translators  sad  F.E.T. 's.  These  odd  or¬ 
der  antuorke  offer  reduced  sensitivity  to  parasitic 
capacity  variations  1m  active  devices,  thus  potenti¬ 
ally  reducing  the  coot  of  monolithic  circuits  by 
eliminating  turning. 


■•••  .let 


(c) 

Fig.  4  (a)  Avaatek  AT-S111  broad-band  F.E.T.  ampli¬ 

fier  with  a  distributed  third  order 
matching  filter. 

(b)  VSVK  of  AT-8111  amplifier  for  Input  net¬ 
works  consisting  of  third  order  matching 
filter  (solid  line),  no  matching  network 
(dashed  lino).  Smith  chart  designed 
matching  network  (alternating  line),  and 
fifth  order  matching  filter  (double 
alternating  line). 

(c)  Change  la  VSVK  of  AT-8111  snpllflar  as  a 
.193  pF  eapacltor  Is  placed  across  lta 
gate-source  leads.  Seme  lime  equiva¬ 
lences  as  in  (b). 


Inferences 

(1]  C.  1.  Matt  heel.  Tables  of  Chabyaav  Impedsnee- 
Tramsformlitg  networks  of  low- Pass  Filter  Form,* 
Free.  IEEE,  vel.  32,  pp.  939-943.  August  19*4. 

(2|  A.  H.  Kiddle  and  t.  J.  Trow,  ’A  Bread-land  Ampli¬ 
fier  Output  Hptwark  Design,*  XEES  Trans.  Mlcro- 
wave  Theory  Tech.,  vel,  HIT- JO,  pp.  192-191, 
Feb.  1M2. 


Acknowledgment 

Aa  always,  the  authors  are  indebted  to  the  ex¬ 
cellent  work  of  Ns.  I.  K.  Tyson,  the  typist. 


i 

i 


461 
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GalnAsP  quaternary  alloys  are  known  to  have  two  or  three  longitudinaL  optical  (LO  (-phonon 
modes.  The  two-LO-phonon  model  was  employed  for  the  electron  distribution  and  drift  mobility 
calculations.  Perturbation  distribution  has  sharp  drops  due  to  the  two-LO  phonons  and  these 
drops  cause  a  further  downward  bowing  of  the  mobility  along  the  arsenic  composition  compared 
with  a  one- LO- phonon  model  which  has  mosdy  been  used. 

PACS  numbers:  72.20.Fr,  63.20.Kr 


There  has  been  a  great  interest  in  GalnAsP  alloys  lat¬ 
tice  matched  to  InP  for  their  applications  to  optoelectronic 
devices 1-1  and  potential  applications  to  high-speed  logic  de¬ 
vices.4-*  Better  performance  in  GalnAsP  than  GaAs  of 
high-frequency  field-effect  transistors  and  high-speed  logic 
devices  is  expected  because  of  its  higher  electron  mobil¬ 
ity0-7**  and  higher  drift  velocity4-10*11  than  GaAs.  Though 
the  experimental  verification  of  these  electron  transport  pro¬ 
perties  is  being  carried  out,  the  mechanism  of  the  electron 
transport  in  this  alloy  is  not  very  clear.  There  are  two  main 
problems.  One  is  alloy  scattering  due  to  the  random  arrange¬ 
ment  of  atoms  in  the  alloy,  and  the  second  is  the  multi¬ 
phonon  modes  observed  in  this  alloy. 12-11  The  problem  of 
alloy  scattering  in  this  alloy  has  been  discussed  by  several 
authors.  °-7-*-1 1-14  However,  they  assumed  the  one-longitu¬ 
dinal  optical  (LOj-phonon  model  in  this  alloy  except  for  Ref. 
14,  and  this  assumption  contradicts  the  experimental  obser¬ 
vation.  We  have  reported  the  theoretical  calculation  of  the 
electron  Hall  mobility  using  the  two-LO-phonon  model  and 
have  shown  a  good  agreement  between  the  calculation  and 
experiments.14 

In  this  letter,  we  report  the  effects  of  the  two  LO- 
phooon  modes  on  the  electron  distribution,  especially  on  the 

Tr— enddnas  Department  of  Electrical  Ensjnawjn»KyoioUaiT«r- 
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perturbation  distribution  in  GalnAsP,  and  clearly  illustrat 
the  reason  the  electron  mobility  calculated  for  the  two-LO 
phonon  model  is  lower  than  that  for  the  one-LO-phonon 
modeL 

The  calculation  method  of  the  mobility  is  the  iterative 
technique. 15-10  The  electron  distribution  function /(k)  unde 
the  influence  of  a  small  electric  field  can  be  expressed  as 
/( It)  rmfjjk )  +  cas0g[k ),  where  fj^k )  is  the  equilibrium  distri 
bution,  cost?  is  the  cosine  between  the  momentum  k  and  elec 
trie  field  F,  and  g[k )  is  the  perturbation  distribution.  The 
Boltzmann  equation  is  solved  for/(k).  After  substituting/^ 
withy^k )  +  cos9g\k ),  the  Boltzmann  equation  is  multiplier 
by  the  first-order  Legendre  polynomial  cosd  and  integrate: 
over  9.  The  remaining  equation  for  the  one-LO-phonon 
model  is 

s&k)-s^kn+s^kn-  II 

where  S0  is  the  total  scattering  rate  and  consists  of  the  scat 
tering-out  rate  by  the  LO-phonon  absorption  (£,£.),  and  b> 
the  LO-pbonon  emission  ^» ),  and  the  differences  betwest 
the  scattering-in  rate  and  scattering-out  rate  for  all  elastic 
scattering  processes  of  piezoelectric  scattering  (5^ ),  acousii 
phonon  scattering  ($„),  ionized  impurity  scattering  (£„),  ant 
alloy  scattering  ($„).  Thus  50»  S£,  +S^  +Sm  +SK 
+  5,  +S*.  These  scattering  rates  for  the  elastic  process  ar> 
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equal  to  the  inverse  of  the  relaxation  times.  S  (5  ^ )  is  the 
scattering-in  rate  by  the  LO-phonon  emission  (absorption). 
The  electron  in  the  state  of  k  *(k  “)  which  is  the  momentum 
evaluated  at  the  energy  corresponding  to  E  plus  (minus)  the 
energy  of  the  LO  phonon  may  be  scattered  into  the 
momentum  state  k  by  the  phonon  emission  (absorption).  The 
expressions  for  all  these  scattering  rates  are  found  in  Refs.  13 
and  17  for  alloy  scattering. 

For  the  two-LO-phonon  model  the  scattering  process 
by  the  LO  phonon  can  be  separated  into  two,  i.e.,  the  scatter¬ 
ing  by  the  InP-like  phonon  and  by  the  Gain  As- like 
phonon.12  If-S*..  (5,^.1  represents  the  scattering-out  rate 
of  the  electron  in  the  momentum  state  *  by  the  absorption  of 
the  InP-like  (GalnAs-like)  LO  phonon,  the  scattering  rate 
includes  a  weighting  factor  of  (1  —  x)(l  — y)  (y). l2- 14  The  cou¬ 
pling  strength  between  the  electron  and  the  LO  phonon  is 
assumed  to  be  proportional  to  the  concentration  ratio  of  InP 
and  GalnAs  in  the  quaternary  alloy  Ga, In,  _,As,P,_,. 
The  Raman  scattering  intensity  by  the  InP-like  and  the 
GalnAs-like  LO  phonon  was  observed  to  be  proportional  to 
the  concentration.12  Similarly,  the  scattering-out  rates  by 
the  emission  of  the  LO  phonons  (S^.S  )  and  the  scat¬ 

tering-in  rates  by  the  emission  or  absorption  of  the  LO  phon¬ 
ons  ($£i»  S'*.  S  ^o)  are  weighted  by  the  corre¬ 

sponding  concentration  ratios.  Thus,  the  Boltzmann 
equation  for  the  two-LO-phonon  model  leads  to 

£ag(&)  +  S  a  )  +  ^ipo|g(Ar  j~) 

+S£o*k  (2) 


where  5g  » IS  ^.»  +  S  +  S^q  +  S  apoo 

+  Sn  +  S^.  The  mode  frequency  of  the  InP-like  LO  phonon 
varies  linearly  from  350  cm-1  and  that  of  the  GalnAs-like 
LO  phonon  increases  from  200  to  275  cm" 1  with  increasing 
y.12  The  drift  mobility  is  calculated  from  the  integration 


1  V  f-Jfc 2  ..  .. 


where  « is  the  free-electron  concentration, 

,  mdE/dk 
~#k 

is  the  augmented  density  of  states, 12  and  Fis  the  crystal 
volume. 

The  calculated  drift  mobilities  at  300  Kfor  the  one-LO- 
phonon  model  and  the  twoLO-phonou  model  are  compared 
in  Fig.  1.  The  electron  concentration  a  is  1 X 10“  cm-2  and 
the  carrier  compensation  ratio  (Af0  +  NA)/n  is  1.  Here,  N0 
and  Na  are  the  donor  concentration  and  acceptor  concentra¬ 
tion  and  are  assumed  to  be  fully  ionized.  These  curves  should 
be  compared  with  the  previous  calculation  which  predicted 
an  upward  bowing  of  the  dirft  mobility  along  the  alloy  com¬ 
position  y. 2  This  was  caused  by  the  lower  effective  mass 
which  was  — 1  by  the  interpolation  method.  The  ex¬ 
perimentally  determined  linearly  varying  effective  mass 
gives  the  downward  bowing  of  the  mobility  curve  which  is 
the  general  observation  in  experiments.  A  farther  bowing  is 
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FIG.  1. Electron  driftmobtlitmcaiculatedbytheitentmtecbniquefortlie 
one- LO-phonon  model  (dashed  line)  and  for  the  two-LO-phonon  model 
(solid  line).  The  electron  concentration  n  is  1 X 10“  cm"1  and  the  carrier 
compensation  ratio  (NB  +  ffA]/n  is  I. 

observed  in  the  two-LO-phonon  model  and  it  is  very  close  to 
the  experimental  results. 

Since  the  difference  in  the  mobilities  for  the  two  models 
is  largest  sty  ■  0.2~0.3,y  *•  0.3,  i.e.,  Gaa^In^ASo  jPa, 
is  chosen  to  discuss  the  effects  of  the  two  LO-phonon  modes 
on  electron  distribution.  Some  of  the  material  parameters  at 
y  «  0.3  used  in  the  calculation  are  listed  in  Table  I.  These 
values  are  estimated  from  the  linear  interpolation  between 
those  of  InP  and  Gag  4T  In* „  As.  In  Fig.  2  the  perturbation 
distribution  g\k )  for  the  one- LO-phonon  model  (dashed  line) 
and  for  the  two-LO-phonon  model  (solid  line)  is  shown  as  a 
function  of  the  normalized  momentum 

r  /  ff2k2/2jw*2V/2 
V  k,T  ) 

at  300  K.  When  the  electron  energy  is  equal  to  or  greater 
than  the  phonon  energy,  the  electron  will  lose  its  energy  by 
emitting  the  phonons.  In  the  one-LO-phonon  model  a  sharp 
drop  at  X  «■  1 .26  is  observed  due  to  the  LO-phonon  emission 
of  the  energy  40.64  meV.  In  the  two-LO-phonon  model 
drops  at  X*  1.10  (GalnAs-like  phonon,  28.33  meV)  end 
X  »  1.30  (InP-like  phonon,  42.60  meV)  are  clearly  illustrat¬ 
ed.  The  fall  in  the  curves  *t  around  X  »  1.18  is  the  perturba¬ 
tion  due  to  the  emission  of  two  LO  phonons.  Since  the  drift 


TABLE  L  Material  parameters  at  y  —  0.3  um4  in  the  calculation. 


Eftctivt  maa  ratio  mVm, 

0.0683 

Optical  ditlactric  ooaMaai «. 

10.01 

Static  tidaenic  coaaunt  a. 

1178 

Energy  gap  at  0  K  (aV) 

1.24 

Optical  pboaoa  tatty  (K) 

oaa-LO-phoaoa  moO*l 

471.6 

two-LO-phoooa  modal  laP-Bka 

4*4.4 

OaInA»4ika 

340-3 

Aopi,  Ptwa.  Lett,  Vol.  40,  No.  9,  i  May  1 9#2 


Ti 


tt. 


837 


48 


M03 


PI O.  I  Pwuitaiw  distribution  m  »  function  of  thn  momentum.  The  ooe- 
LO- phonon  model  i>  the  dashed  line  and  the  two-LO-phooon  is  the  solid 
lint.  The  sharp  dropa  in  the  distribution  an  due  to  the  phonon  emission.  In 
the  two-LO-phonon  model  it  has  two  dropa  corresponding  to  InP-like  and 
GalaAsdika  phonons. 

mobility  is  proportional  to  the  integral  of  the  perturbation 
distribution  as  is  clear  from  Eq.  (3),  the  difference  of  these 
curves  gives  us  the  drift  mobility  differences;  3343 
cm2V-ls~'  for  the  one-LO-phonon  model  and  3132 
cm2V-,s" 1  for  the  two-LO-phonon  model.  For  the  calcula¬ 
tion  of  these  mobilities  the  electron  concentration  1 X  101S 
cm-3  and  the  carrier  compensation  ratio  l  were  used.  The 
degree  of  a  drop  depends  on  the  strength  of  the  electron- 
phonon  coupling.  It  will  be  straightforward  to  extend  the 
model  to  the  three  or  more  LO  phonons  as  observed  by  the 
far-infrared  reflectivity  study,13  although  the  coupling, 
strength  of  the  electron  to  these  phonons  is  still  in  question. 

In  summary,  the  electron  drift  mobility  and  the  elec¬ 
tron  distribution  in  Ga,In,  As,P,  lattice  matched  to 
InP  have  been  calculated  by  the  iterative  technique.  The 
electron  mobility  is  lower  for  the  two-LO-phonon  model 
than  for  the  one-LO-phonon  model  over  the  whole  composi¬ 


te  Appi.  Ptiys.  Lao..  Vot.  ao,  No.  0. 1  Mpy  1M2 


non  y.  The  perturbation  distribution  in  the  two-LO-phooon 
model  shows  two  sharp  drops  at  the  energies  equal  to  the  two 
LOphonon  energies  which  result  with  the  lower  electron 1 
mobility  compared  with  the  one-LO-phonon  model 

This  work  was  partially  supported  by  Army  Research 
Office — Durham  under  Grant  Number  DAAG29-80-K- 
0080,  and  partially  by  the  Office  of  Naval  Research,  Ariing 
ton,  Virginia. 
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ABSTRACT 

This  paper  presents  a  study  of  the  spatial  ensemble  velocities  of 
electrons  in  GaAs  which  result  from  abrupt  changes  in  electric  fields 
occurring  internal  to  a  device.  Calculations  are  presented  for  single  field 
spikes  as  well  as  the  complex  field  configuration  which  occurs  in  planar  doped 
barrier  transiscor  type  devices.  The  calculations  provide  estimates  of 
currenc  gain  and  carrier  transit  times  for  PDB  transistor  structures.  Hign 
current  gains  are  found  to  require  very  narrow  base  widths.  ~- 


INTRODUCTION 

The  ability  to  grow  very  thin  semiconductor  layers  with  accurately  con¬ 
trolled  doping  has  brought  about  new  device  concepts  in  recent  years  [1,2].  To 
improve  high  frequency  performance  many  of  the  new  device  concepts  attempc  to 
use  non  steady-state  electron  dynamics  such  as  ballistic  transport  or  velocity 
overshoot  which  occur  over  short  distances  when  charge  carriers  are  subjected 
to  large  electric  fields.  The  GaAs  planar  doped  barrier  (PDB)  illustrated  in 
Fig.  I  is  a  rectifying  junction  [2]  which,  when  fabricated  with  an  accurately 
controlled  doping  profile,  should  demonstrate  hot  electron  characteristics. 

By  combining  two  barriers,  a  PDB  transistor  structure  can  be  formed  as 
shown  in  Fig.  2.  The  basic  idea  is  that  carriers  injected  by  forward  bias 
over  the  larger  of  the  two  barriers  will  become  hot  carriers  in  the  heavily 
doped  base.  Most  of  tnese  carriers  should  have  high  velocity  and  sufficient 
energy  to  traverse  the  smaller  collector  barrier.  Interest  in  such  device 
structures  has  provided  the  motivation  for  the  present  ensemble  Monte  Carlo 
study  of  carriers  subjected  co  abrupt  changes  in  electric  field. 


•This  work  was  supported  by  a  research  grant  from  the  Office  of  Naval 
Research,  Arlington,  VA  and  contract  No.  DAAG29-80-K-0080  from  the 
Army  Research  Office,  Durham,  NC. 
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ABKUPT  FIELD  CHANGES  AND  HIGH  FIELD  SPIKES 
Ensemble  velocity  profiles  near  an  abrupt  change  in  the  electric  field 
are  shown  in  Fig.  3.  Three  examples  are  shown  with  the  field  for  x  >  0  being 
20  kV/cm  and  for  x  <  0  taking  on  the  three  values  of  2,  5  and  8  kV/cm.  The 
curves  in  Fig.  3  show  the  general  behavior  of  velocity  overshoot  expected  as 
carriers  enter  a  high  field  region.  The  high  velocity  transient,  however,  be¬ 
gins  a  considerable  distance  before  the  carriers  enter  the  high  field  region 
with  the  velocity  exceeding  2x107  cm/sec.  Since  the  '.trs  cannot  have  their 
energy  changed  by  the  high  field  before  exposure  to  ch  »h  field,  this  ini¬ 
tial  velocity  increase  cannot  be  due  to  hot  carrier  effects. 

A  decalled  study  of  this  pre-overshoot  effect  has  demonstrated  that  it  is 
not  a  hoc  carrier  effect,  but  an  ensemble  effect  on  the  carriers  for  x<0  due 
to  the  carrier  boundary  condition  established  by  the  high  field  region  begin¬ 
ning  at  x-0.  Basically,  Che  low-field/high-field  interface  acts  as  an  ideal 
carrier  collecting  boundary  to  Che  low  field  region.  Thus,  any  carrier  reach¬ 
ing  the  high  field  region  has  an  extremely  low  probability  of  being  back  scat¬ 
tered  into  Che  low  field  region.  The  velocity  distribution  for  x  slightly 
less  chan  zero  must  have  essentially  no  particles  with  negacive  velocities. 

The  above  discussion  provides  a  simple  explanation  for  the  ensemble  aver¬ 
age  velocity  near  Che  low-field/high  field  interface.  Consider  a  drifted  Max¬ 
wellian  velocity  distribution  function  with  a  mean  velocity  of  v0.  a  calcula¬ 
tion  of  che  ensemble  average  velocity  of  only  the  positive  velocity  particles 
yields 


<v>+ 


exp(-£o/kT) 

- +  Vq 

1  ♦  erf  (/S0/kT) 


L 


(l) 
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where  £0  -  (l/2)m*v02.  This  simple  expression  has  been  found  to  give  very 
good  agreement  to  the  low-field/high-field  interface  ensemble  velocity  when 
evaluated  with  appropriate  vQ  and  kT  parameters. 

Confirmation  of  the  above  model  is  provided  in  Fig.  4  which  shows  the  en¬ 
semble  velocity  for  a  high  field  spike  of  width  3Q0A  beginning  at  x  *  0.  The 
circular  points  are  velocities  calculated  for  x  <  0  when  the  x  •  0  Interface 
is  replaced  by  an  ideal  collecting  contact.  The  close  agreement  between  the 
two  curves  verifies  chat  the  low-field/high-field  boundary  is  acting  like  an 
ideal  collecting  boundary.  The  value  of  <v>+  »  2.5x10?  cm/sec  calculated  from 
(1)  is  in  excellent  agreement  with  the  interface  ensemble  velocity.  It  is  al¬ 
so  noted  in  Fig.  4  chat  the  velocity  to  the  right  of  the  high  field  spike  is 
about  2  x  10?  cm/sec  and  considerably  above  Che  1  kV/cm  steady  state  value 
(Vss)  of  about  0.8  x  10?  cm/sec.  Thus,  the  high  field  spike  has  imparted  a 
considerable  component  of  x  directed  momentum  to  the  ensemble  of  particles  and 
this  persists  for  a  considerable  distance  beyond  the  high  field  spike. 

PLANAR  DOPED  8ARRIER  TRANSISTORS 

The  assumed  field  profile  and  Che  resulting  velocity  profile  for  a  PDB 
transistor  structure  is  shown  in  Fig.  5.  This  particular  structure  has  an 
emitter  barrier  height  of  0.31  volts  and  calculations  are  shown  for  three  col¬ 
lector  barrier  heights  of  0.05,  0.10  and  0.15  volts,  which  give  barrier  height 
differences  of  0.16,  0.11  and  0.06  volts.  These  calculations  are  for  a  base 
width  of  2000A.  Other  calculations  have  been  made  for  base  widths  of  250A  and 
1000 A. 

The  carrier  velocity  within  the  base  region  (0.03  ua  <  *  <  .23  pm)  is 
lower  chan  might  be  expected  from  simple  ballistic  arguments  on  single  car¬ 
riers.  The  velocity  profile  around  the  emitter  field  spike  is  similar  to  chat 
for  a  single  field  spike  aa  shown  in  Fig.  4.  There  is  an  initial  rise  to 
above  2  x  107  ca/see  before  the  high  field  region  followed  by  a  very  narrow 
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peak  la  Che  velocity.  The  velocity  within  the  base  layer  is  lower  than  that 
following  a  single  field  spike  because  of  some  additional  trapping  of  carriers 
between  che  two  energy  barriers  within  the  base  layer. 

The  collector  region  in  the  present  calculations  was  taken  as  18  kV/ca 
and  of  length  5000A.  The  velocity  shows  a  second  peak  in  this  region  (x  > 
0.25  ua)  which  is  followed  by  a  region  of  low  velocity  as  carriers  transfer  to 
the  upper  (L)  valley.  The  velocity  in  che  collector  region  is  found  to 
depend  on  the  height  of  the  collector  barrier  as  shown  by  the  three  curves. 

Soae  carriers  lose  sufficient  energy  within  che  base  layer  to  become  per— 
aanancly  crapped  there.  (Carriers  which  transfer  to  the  L  valley  will  become 
temporarily  crapped  until  they  transfer  back  to  che  r  valley).  In  the  present 
calculations  che  trapped  carriers  constitute  the  base  current  and  were  removed 
from  che  enseable  as  soon  as  they  reached  the  collector  barrier  with  insuffi¬ 
cient  energy  to  traverse  it. 

Fig.  6  shows  calculated  beta  values  as  a  function  of  base  width  for  two 
barrier  height  differences.  As  would  be  expected,  beta  is  a  strong  function 

of  both  base  width  and  barrier  height  difference.  The  solid  curves  are  based 

upon  a  simple  model  which  assumes  that  Ic  -  Ie  exp(-xjj/L)  where  L  is  a  pheno¬ 
menological  parameter  which  was  evaluated  from  che  1000A  base  width  calcula¬ 
tions.  The  Monte  Carlo  calculations  at  250A  and  2000A  are  seea  to  agree  rea¬ 
sonably  well  with  this  simple  model.  Beta  values  above  10  have  been  found  to 
require  base  widths  below  1000A. 

The  average  electron  transit  time  for  the  transistor  structures  shorn  in 
Fig.  5  were  in  the  range  of  5-6  psec  yielding  gain-bandwidth  products  of  ft  - 
30  GHx. 
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Figure  5.  Ensemble  Electron  Velocities  for  the  PDB  Transistor  with 

Collector  Barrier  Heights  of  0.05  volts  (curve  a),  0.1  volts 
(curve  b)  and  0.15  volts  (curve  c)  and  an  Emitter  Barrier  of 
0.21  volts. 
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ABSTRACT 


A  study  of  ion-implanted  MESFET  performance  aa  a  function  of  the 
implantation  energy  and  fluency  and  including  the  effects  of  deep-level 
trapping  state  concentrations  in  the  substrate  has  been  conducted.  Carrier 
concentrations  as  a  function  of  depth  are  determined  through  the  use  of  LSS 
theory  and  a  profiling  model.  An  analytic  device  model,  which  computes  both 
DC  and  RF  characteristics,  is  then  employed  to  predict  MESFET  performances. 

The  study  includes  the  effects  of  depth  dependent  transport  properties  and  has 
indicated  a  number  of  design  rules  for  the  fabrication  of  optimized 
ion-implanted  devices. 


*This  work  was  supported  by  Rockwell  International  Electronics  Research  Center 
and  the  U.S.  Army  Research  Office,  Research  Triangle  Park,  NC  on  contract 
DAAG29-80-K-0080 . 
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INTRODUCTION 

The  performance  of  FETs  fabricated  by  ion-implantation  depends  greatly  on 
carrier  concentrations  and  velocity-field  characteristics  as  a  function  of 
depth  into  the  active  device  layers.  The  presence  of  deep-level  traps  in  the 
semiconductor  contributes  to  the  complexity  of  problems  associated  with  the 
characterization  of  ion-implanted  devices. 

This  paper  presents  the  results  of  a  study  to  determine  the  effects  of 
various  concentration  and  transport  profiles  upon  device  performance.  The 
study  utilizes  theoretical  models  of  both  material  properties  and  device 
characteristics.  This  Information  is  combined  with  experimental  material 
characterization  to  provide  improved  quantitative  accuracy  of  the  model. 

Figure  1  outlines  the  modeling  steps  Involved  in  obtaining  the  results  to 
be  presented.  Initially,  the  determination  of  typical  ranges  for  concentra¬ 
tions  of  deep-level  trapping  states,  Ny(x),  was  made  from  a  novel  measurement 
technique  using  both  differential  capacitance  and  conductance  DLTS  data. 
Shallow-level  donor  concentrations,  Nq(x),  were  then  determined  from  LSS 
theory  as  a  function  of  implantation  parameters.  These  concentrations  as  a 
function  of  depth  were  then  used  in  a  profiling  model  to  determine  the 
free-carrier  profile,  n(x),  for  the  material.  Carrier  transport  properties 
were  also  determined  from  the  trap  and  donor  profiles.  This  was  done  through 
the  use  of  Monte  Carlo  particle  simulations  and  a  model  to  account  for  the 
effects  of  compensation  in  the  semiconductor. 

An  analytic  model  which  utilizes  all  of  this  information  is  then  used  to 
assess  performance  potential.  The  DC  characteristics  and  mall-signal 
S-parameters  along  with  figures-of-merlt  are  computed  by  the  model.  The 
distinction  between  free-carrier  and  donor  profiles,  the  effect  of  deep-level 
traps,  and  the  depth  dependence  of  mobility  and  velocity  are  shown  to  be 
important  considerations  which  have  been  ignored  in  previous  models. 
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THE  DEVICE  MODEL 

A  one-dimensional  analytic  device  model  with  a  small-signal  analysis  has 
been  developed  and  utilized  for  this  work.  The  model  Is  based  upon  the 
principles  presented  by  Pucel  [1]  but  has  been  generalized  for  arbitrary 
doping  profiles  and  Includes  differences  in  the  free-carrier  and  shallow-level 
donor  profiles  and  the  effects  of  deep- levels. 

The  model  assumes  that  the  electron  transport  properties  of  a  material 
can  be  simulated  by  a  two-piece  velocity-field  relationship.  The  two-piece 
approximation  is  defined  from  a  theoretical  velocity-field  characteristic 
determined  by  Monte  Carlo  techniques.  For  electric  fields  less  than  an 
appropriate  saturation  field,  Ea,  the  electron  velocity  Is  described  by  a 
linear  expression, 

v  •  PqE.  (1) 

For  electric  fields  above  Ea  the  electrons  move  at  a  constant,  maximum 
velocity,  va. 

The  low-field  mobility  as  a  function  of  depth  Into  the  channel  is 
determined  from  the  equation  [2], 
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with  k  -  0.025, 

B  -  0.817278, 

C  -  6.252838, 

and  where  %  Is  given  In  m“3. 

Equation  (2)  glvea  the  low-field  nobility  aa  a  function  of  donor  denalty 
and  background  compensation.  To  derive  equation  (2)  the  theoretical  results 
of  Walukiewlcz  at  al.  [3]  were  used  In  conjunction  with  Monte  Carlo 
velocity-field  predictions  [A].  The  Monte  Carlo  results  were  used  to 
determine  the  mobility  as  a  function  of  background  donor  density  with  no  traps 
present,  and  the  Walukiewlcz  values  were  then  normalized  to  the  Monte  Carlo 
numbers.  The  normalized  data  were  finally  curve  fit  to  obtain  equation  (2). 
Very  good  agreement  Is  obtained  when  equation  (2)  is  plotted  against  the 
normalized  Walukiewlcz  values  aa  shown  in  Figure  2. 

The  technique  for  determining  the  maximum  velocity  vB  has  been  previously 
presented  [5,6,7].  The  saturation  velocity  can  be  described  by  the 
expression 

vm  ■  vQ  -  A  log  [(1-y)2  +  By],  (3) 

where 

y  -  [ND(x)/N0]2-5, 

H0  -  1.3  x  1022  a-3, 

k  -  0.0262, 

B  -  0.4, 

and  v0  -  1.40  x  103  m/sec. 

Equation  (3)  gives  va  in  103  m/sec  when  Nj)(x)  is  expressed  in  m~3.  The 
expression  Is  assumed  to  have  the  same  dependence  upon  compensation  ratio  9  as 
equation  (2).  Thus,  the  factor  (l-9)b  is  multiplied  with  equation  (3)  to 
obtain  va  in  the  presence  of  traps.  The  exponent  b  is  defined  In  equation 
(2) ' 
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Equations  (2)  and  (3)  in  conjunction  with  knowledge  of  the  three 
profiles,  n(x)  -  the  free-electron  concentration,  %(x)  -  the  donor 
concentration  and  Nf(x)  -  the  deep-level  concentration,  allow  for  the 
derivation  of  a  device  model  which  includes  the  effects  of  varying  transport 
properties  as  a  function  of  epi-depth.  The  profiles  of  interest  are 
determined  with  a  novel  characterization  technique  [2]  that  combines 
experimental  capacitance-voltage  (C-V)  and  conductance  deep-level  transient 
spectroscopy  (DLTS)  data.  The  experimental  equipment  and  procedure  is 
discussed  elsewhere  [8]. 

For  this  work  an  ion-implanted  one  micron  gate  length  device  [9]  was 
analyzed.  The  experimental  characterization  revealed  the  presence  of  a 
dominant  deep-level  donor  state  0.736  eV  below  the  conduction  band.  The  level 
is  probably  EL2  [10,11]  in  agreement  with  the  findings  of  Martin  et  al.  [12]. 
The  resulting  f ree-carrler ,  shallow-level  donor,  and  deep-level  concentrations 
as  a  function  of  depth  into  the  material  are  shown  in  Figure  3.  Despite  the 
scatter  of  the  data  near  the  tail,  it  is  evident  that  diffusion  of  the 
electrons  has  caused  the  f ree-carrler  and  shallow-level  donor  concentrations 
to  differ  significantly.  The  low-field  mobility  versus  depth  obtained  from 
this  analysis  is  also  shown  in  Figure  3.  In  regions  of  relatively  high  donor 
concentration  the  mobility  varies  inversely  with  doping,  as  expected.  Near 
the  tall  of  the  implanted  region,  however,  the  deep-levels  seriously  degrade 
the  mobility  acting  to  confine  the  active  device  channel. 

The  dc  I-V  characteristics  for  the  device  as  determined  from  the  model 
are  compared  in  Figure  A  with  the  measured  Z-V  curves.  It  should  be  noted 
that  without  the  inclusion  of  the  effects  of  traps  upon  the  carrier  transport 
characteristics  the  excellent  agreement  shown  in  Figure  4  could  not  be 
obtained.  The  traps  have  a  tendency  to  "soften”  the  pinch-off  characteristics 
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of  che  device  (i.e.,  when  crape  are  included  in  che  simulation,  the  slope 
dIE/dVg  is  noc  as  greac  near  pinch-off).  For  all  of  che  devices  sCudied  in 
Chis  work,  this  "softening"  effect  was  required  co  obtain  best  agreement  with 
Che  measured  I-V  characteristics. 

The  model  allows  Che  element  values  for  a  small-signal  equivalent  circuit 
aa  shown  in  Figure  5  co  be  determined.  The  equivalent  circuit  can  be  analyzed 
to  obtain  RF  performance  predictions  as  a  function  of  the  various  Implantation 
parameters  and  material  characteristics. 

RESULTS 

The  results  that  follow  were  obtained  following  the  modeling  steps 
outlined  in  the  Introduction  and  illustrated  in  Figure  1.  Three  parameters 
were  varied  independently.  They  are  1)  implantation  energy,  e,  2)  peak  doping 
density,  Hm-it.  and  3)  trapping  state  density,  %(x).  The  implant  species  was 
assumed  to  be  Si  in  GaAs  and  Che  activation  was  assumed  to  be  100X  for  all 
devices.  The  gate  length  of  the  simulated  devices  is  one  micron  with  a  gate 
width  of  300  microns.  The  trapping  state  density  was  assumed  to  be  constant 
as  a  function  of  depth  for  these  studies.  Note  that  the  peak  doping  density 
can  be  converted  to  a  corresponding  ion  fluency  through  the  simple  relation¬ 
ship 

e  -  (4) 

where  <Jp  is  the  standard  deviation  of  the  projected  range.  The  value  for  <3p 
can  be  obtained  from  LSS  theory  when  the  implant  schedule  is  known. 

Figure  6a  shows  computed  carrier  profiles  for  che  same  75  keV  implant 
with  -  2x10^  cm~3,  but  for  three  different  trapping  state 
concentrations.  The  low-field  mobility  and  maximum  velocity  profiles 
corresponding  to  these  three  cases  are  shown  in  Figure  6b.  It  is  clear  from 
the  figures  that  the  trapping  state  concentration  in  the  material  has  a 
significant  effect  on  both  the  free-carrier  concentration  and  che  transport 
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properties  of  the  implanted  material.  The  equivalent  circuit  element  values 
that  correspond  to  these  cases  are  shown  in  Table  I.  The  deep-level 
concentration  has  the  most  effect  upon  the  device  transconductance  which 
decreases  as  the  trap  concentration  increases  and  the  gate  delay  time  which 
increases  with  the  trap  concentration. 

Table  I 

Ion-Implanted  MESFET 
Equivalent  Circuit  Element  Values 


e  -  75  keV,  Npealt  -  2xl017  cm*  3 
Vda  "  3.5  V,  Id  -  Id88 

_ Deep-Level  Concentration 


Element 

lxlO14  cm-3 

2xl0^3cm-3 

lxl0l®cm“3 

ga(mmho) 

46.4 

45.6 

41.5 

Cgs(pF) 

0.419 

0.417 

0.399 

Rdr(kfl) 

1.22 

1.20 

1.35 

Ri(Q) 

3.337 

3.423 

3.680 

•c(psec) 

8.00 

8.16 

8.41 

fmax(®*) 

66.1 

64.5 

61.8 

The  remaining  parameters  were  constant  with  the  values  C^8  »  0.0577  pF, 
Cgd  -  0.0432  pF,  Rd  *  1.214  Q  and  Rg  ■  1.214  Q. 


For  Che  implant  energy  and  peak  doping  studies,  the  crapping  state 
density  was  left  constant  at  N  -  2x10^3  cm~3.  This  number  was  chosen  to  be  in 
general  agreement  with  Che  results  shown  in  Figure  3.  The  implant  energy  was 
varied  from  50  to  150  keV  while  the  peak  doping  took  values  between  8xl0l$  and 
4x1017  cm"3. 
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The  value  for  f?  Is  computed  from  first  order  considerations  to  be  given  by 

fT  "  2^7  *  <5) 

The  quantity  fm^  Is  obtained  by  noting  the  frequency  at  which  Mason's 
unilateral  gain  (as  predicted  by  the  model)  goes  to  unity  .  Mason's 
collateral  gain  at  a  frequency  of  10  GHz  is  plotted  as  a  function  of  energy  In 
Figure  8.  The  results  clearly  indicate  the  superiority  of  lov-energy  implants 
for  high-frequency  operation. 

The  information  plotted  in  Figure  8  indicates  the  importance  of  Implant 
energy  to  achieve  optimum  high  frequency  performance.  Notice  that  for  a  peak 
doping  of  4xl0l7  cm~3,  a  decrease  of  implant  energy  from  150  to  50  keV  results 
in  better  than  a  92  increase  in  the  unilateral  power  gain,  from  16.2  dB  to 
17.7  dB.  The  superior  high  frequency  performance  of  low  energy  implants  is 
easily  understood  in  terms  of  the  device  physics.  The  available  charge 
carriers  lie  closer  to  the  surface  (gate  contact)  for  low  energy  implants. 
Thus,  the  gate  potential  required  to  alter  the  depletion  width  is  smaller  than 
for  deep  implants  and  the  transconductance  of  the  device  should  rise.  High 
transconductance  is  important  to  the  high  frequency  performance  of  these 
devices.  All  of  the  devices  considered  in  compiling  Figures  7  and  8  were 
compared  at  a  bias  of  Iq  -  10  mA.  The  dashed  line  falling  off  rapidly  at  the 
low-energy  end  of  the  curves  indicate  that  for  energies  lower  than  this, 
IDSS<10  mA. 

The  trapping  state  concentration  in  ion-implanted  semiconductor  material 
can  vary  over  a  wide  range  of  values.  The  quality  of  the  initial 
semi-insulating  substrate,  the  temperature  characteristic  of  the  various 
processing  steps  used,  and  the  type  of  annealing  implemented  all  affect  the 
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kinds  and  relative  concentrations  of  deep-levels  in  the  final  device  [13,14]. 
For  the  trap  study  in  this  work,  an  inplant  energy  of  75  keV  and  a  peak  doping 
density  of  2  x  LO^  ca~3  was  assumed .  The  background  trapping  state 
concentration  was  varied  from  0  to  10^6  cm“^.  p0r  trapping  statfi  concentra¬ 
tions  above  10^6  ca“^,  the  device  would  be  normally  pinched-off. 

Figure  9  illustrates  the  effects  traps  have  on  the  zero  gate  bias 
current,  Iqss>  and  the  pinch-off  potential,  W0o*  The  trapping  state  density 
of  the  figure  ranges  between  10^  and  10^6  cm“^.  Over  this  range  the  zero 
gate  bias  current  varies  between  30  and  18  mA,  while  the  pinch-off  potential 
varies  between  1.62  and  1.38  volts.  This  large  range  of  DC  characteristics 
indicates  that  reproducibility  will  be  a  problem  unless  consistency  of  the 
trapping  state  density  in  the  substrate  can  be  maintained.  It  should  be 
further  noticed  that  for  deep-level  concentrations  below  a  level  of  about  lO1^ 
cm-^,  very  little  change  in  the  characteristics  occurs.  As  the  trapping  state 
concentration  Increases  above  the  low  lO1^  cm”3  ieVel,  however,  the 
current  and  pinch-off  potential  begin  to  fall  rapidly. 

Figure  10  shows  the  pinch-off  "softening'*  effect  mentioned  in  the 
previous  section  for  one  particular  device.  Notice  that  the  figure  plots 
normalized  current  and  voltage.  Thus,  the  trapping  state  density  affects  the 
DC  characteristics  of  the  device  in  a  qualitative  as  well  as  a  quantitative 
way. 

Deep-level  traps  also  have  some  effect  on  the  RF  performance  of  the 
device  as  illustrated  in  Figures  11  and  12.  Figure  11  shows  clearly  the 
decrease  of  f^  and  associated  with  increasing  trapping  state 
concentracions—eapecially  for  deep-level  concentrations  above  the  low 
1015 

cm~3  range.  In  Figure  12  note  Chat  for  low  bias  currents  fx  Increases 
when  few  traps  are  present  while  it  decreases  for  higher  trap  concentrations. 
This  is  easily  explained  in  terms  of  Che  degrading  effects  craps  have  on 
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mobility  and  velocity.  As  the  gate  bias  restricts  current  flow,  a  larger 
fraction  of  the  carriers  are  forced  deeper  into  the  channel.  This  corresponds 
to  the  more  lightly  doped  regions  of  the  device.  If  the  compensation  ratio  is 
fairly  small  (i.e.,  few  traps)  then  from  equations  (1)  and  (3)  the  transport 
properties  are  superior,  and  f?  increases.  If,  on  the  other  hand,  the 
trapping  state  density  is  on  the  same  order  of  magnitude  as  the  shallow-level 
donor  concentration,  then  the  compensation  ratio  approaches  one.  This 
corresponds  to  extreme  degradation  of  mobility  and  velocity  and,  thus,  forces 
ff  to  decrease. 

CONCLUSIONS 

It  has  been  shown  how  deep  level  traps  in  ion-implanted  devices  degrade 
carrier  transport  properties  in  the  semiconductor  material.  The  degradation 
is  more  severe  near  the  tail  of  the  implant  profile  than  near  the  peak.  Thus, 
the  transport  properties  of  the  device  will  be  depth — or  bias — dependent.  A 
modeling  technique  which  accounts  for  this  dependence  has  been  used  to  study 
device  properties  as  a  function  of  fabrication  parameters  and  deep-level 
concentrations.  The  results  indicate  that  low  energy  Implants  should  possess 
superior  high-frequency  properties  and  that  lowering  the  trapping  state 
density  in  the  material  should  improve  device  performance.  For  the  75  keV 
implant  studied  in  this  work,  a  critical  deep-level  concentration  of  about 
1015 

cm~3  was  identified.  For  trapping  state  concentrations  above  this  level 
performance  degradation  becomes  increasingly  severe.  Decreasing  deep-levels 
in  ion-implanted  devices  below  this  critical  trapping  state  concentration 
should  Improve  device  performance.  It  should  be  noted,  however,  that  for  a 
different  implant  schedule  than  the  one  considered  here  (i.e.,  75  keV  with 
Mbit  ■  2xl017  cm )  the  critical  trap  density  may  be  different. 
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Figure  1  Flow-chart  of  modeling  steps  required  to  obtain  study  results. 

Figure  2  Low-field  nobility  as  a  function  of  doping  and  deep-level 

compensation.  The  solid  lines  are  computed  from  equation  (2)  for  0 
■  0.0,  0.2,  0.4,  0.6,  0.8.  The  data  points  are  from  the  normalized 
theoretical  computations  [3]. 

Figure  3  Resulting  concentration  profiles  and  low-field  mobility  profile  for 
one  device. 

Figure  4  Comparison  of  model  predicted  and  measured  I-V  characteristics  for  a 
one  micron  ion-implanted  device. 

Figure  5  The  equivalent  circuit  for  an  FET  used  in  the  analysis. 

Figure  6a  Concentration  profiles  for  a  75  keV  implant  with  ILt  “  2x101?  Cm~^ 
and  with  Hy  -  lxlO14,  2xl015,  lxlO16  cm~3.  The  solid  curve  is  the 
donor  density  while  the  dashed  curves  are  free-carrler  densities 
appropriate  for  the  various  trapping  state  densities. 

Figure  6b  Transport  property  profiles  for  the  three  75  keV  implants  of  Figure 
6a.  The  solid  curves  give  low-field  mobility  as  a  function  of 
depth.  The  dashed  curves  give  maximum  velocity  as  a  function  of 
depth. 

Figure  7a  Predicted  gain-bandwidth  product,  fy,  implantation  energy,  e. 
The  peak  doping,  S„.T .  is  used  as  an  Independent  parameter. 

Figure  7b  Predicted  maximum  frequency  of  oscillation,  vs.  Implantation 

energy,  e.  The  peak  doping,  H,,T.  is  used  as  an  Independent 
parameter . 

Figure  8  Mason's  unilateral  gain,  0,  at  10  GHz  vs.  energy,  e.  trapping  state 
density  is  constant  at  Ny  ■  2x10^  cm-^. 

Figure  9  Predicted  effects  of  varying  trapping  state  concentrations  on  DC 
characteristics . 

Figure  10  normalized  drain  current  predictions  vs.  normalized  gate  voltage 
with  and  without  traps. 

Figure  11  Predicted  gain-bandwidth  product,  fy,  and  maximum  frequency  of 
oscillation,  f..,.  vs.  background  trapping  state  concentration. 

Figure  12  Predicted  gain-bandwidth  product,  fy,  as  a  function  of  normalised 

drain  current.  The  background  trapping  state  concentration  is  used 
as  a  variable  parameter. 
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A  TECHNIQUE  FOR  MODELING  ION- IMPLANTED  CaA»  aESFETS 
IN  THE  PRESENCE  OF  DEEP  LEVELS 


J.  M.  Golio,  G.  N.  Maracas,  D.  Johnson,  R.  J.  Trev 
and  N.  A.  Maanari 
North  Carolina  State  University 
Electrical  and  Computer  Engineering  Department 
Raleigh,  North  Carolina  27650 


INTRODUCTION 


Recent  Interest  in  developing  GaAs  integrated  circuits  for  both 
microwave  [1]  and  high-speed  logic  [2]  applications  has  stimulated  the 
development  of  ion-implanted  MESFETs.  Obtaining  performance  predictions  for 
ion-implanted  FETs,  however,  presents  several  problems  unique  to  these 
devices.  Free-carrler  and  shallow-level  donor  concentrations  will  not  be 
equivalent  within  the  channel  of  an  ion-implanted  device,  and  both 
concentrations  will  vary  as  a  function  of  channel  depth.  Deep-levels  will 
also  exist  within  the  energy  gap  of  the  channel  material.  These  levels  will 
act  to  decrease  the  number  of  carriers  available  for  conduction  and  also  to 
degrade  transport  characteristics.  In  order  to  obtain  accurate  performance 
predictions,  therefore,  it  is  first  necessary  to  quantitatively  characterize 
each  of  these  depth  dependent  material  properties.  In  many  cases  traditional 
characterization  techniques  will  not  provide  sufficient  information  for  an 
accurate  material  characterization. 

This  paper  presents  a  modeling  technique  for  determining  free-carrler, 
shallow-level  donor,  deep-level  and  low-field  mobility  profiles.  The 
technique  makes  use  of  differential  capacitance  (C-V),  conductance  DLTS  and 
channel  mobility  profile  data.  The  effect  of  this  more  accurate  modeling  on 
performance  prediction  is  also  examined. 

DIFFERENTIAL  CAP AC ITANCE 


Typical  free-carrler,  background  donor,  and  deep-level  profiles  for  a 
sample  of  ion-implanted  semiconductor  material  is  shown  in  Figure  I. 
Free-carrler  diffusion  from  highly  doped  to  more  lowly  doped  regions  in  the 
material  causes  the  free-carrler  profile  to  differ  from  that  of  the 
shallow-level  donor  profile.  This  effect  is  most  pronounced  in  shallow 
implants  and  can  lead  to  as  much  as  an  order  of  magnitude  difference  between 
the  two  concentrations  for  Implants  typically  used  in  the  fabrication  of  GaAs 
MESFETs. 

Traditional  C-V  analysis  establishes  the  distribution  of  majority 
carriers  rather  than  the  distribution  of  ionized  impurity  atoms  [3].  The 
actual  impurity  distribution  in  a  measured  sample  can  then  be  computed  from 
the  free-carrler  distribution  [4].  In  the  C-V  measurement  process,  the 
capacitance  is  normally  measured  with  a  small  oscillating  voltage 
superimposed  over  a  fixed  applied  bias.  The  fixed  voltage  bias  translates 
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into  a  depletion  length  Into  the  material  through  the  use  of  the  abrupt 
depletion  approximation.  The  oscillating  voltage  Induces  a  current  due  to  an 
Incremental  charge  being  alternately  covered  and  uncovered  near  the  depletion 
region  edge. 

The  presence  of  deep-levels  in  the  material  adds  complications  to  the 
analysis.  The  traps  allow  the  Incremental  charge  which  Is  measured  to  come 
not  only  from  free-carrlers  In  the  conduction  band,  but  also  from  carriers 
trapped  in  the  deep-level  sites.  The  interpretation  of  C-V  data  under  these 
conditions  has  been  examined  in  detail  by  Klmerllng  [S] .  The  analysis 
considers  the  relationship  between  the  frequency  of  the  oscillating  voltage, 
the  electron  emission  rate  associated  with  the  deep  energy  site,  and  the  rate 
of  change  of  the  dc  bias  as  well  as  the  type  of  trap  site  present.  Following 
Klmerllng,  expressions  to  accurately  interpret  C-V  data  and  obtain  the 
free-carrier  distribution  can  be  developed  provided  the  dominant  trap  site 
can  be  accurately  described  [6,7].  The  effect  that  traps  have  on  the 
interpretation  of  C-V  data  is  shown  in  Figure  2.  The  solid  curve  in  the 
figure  represents  the  computed  free-carrier  concentration  using  traditional 
C-V  analysis  and  assuming  no  deep-levels  are  present  in  the  material  [3] . 

The  dashed  curves  represent  what  the  actual  free-carrier  distribution  in  the 
material  would  be  for  various  constant  deep-acceptor  and  deep-donor  trap 
densities.  In  actual  ion-implanted  material  the  trap  distribution  will  not 
be  constant  with  depth  so  that  more  detailed  analysis  is  necessary. 

DEEP-LEVEL  TRANSIENT  SPECTROSCOPY 

Conductance  DLTS  similar  to  that  proposed  by  Adlerstein  [8]  has  been 
utilized  in  this  work.  In  the  measurement  process  the  source-drain 
conductance  of  an  FET  structure  is  monitored  while  a  gate  voltage  pulse  is 
applied  repetitively.  Band  diagrams  which  correspond  to  times  before,  during 
and  after  the  gate  pulse  are  shown  in  Figure  3.  A  dc  reverse  bias  voltage  has 
the  channel  nearly  pinched-off  before  the  gate  pulse  is  applied.  The 
drain-source  voltage  la  kept  small  to  Insure  that  the  device  remains  in  the 
linear  region.  Thus,  a  small  steady-state  conductance  is  measured.  Uhen  the 
gate  pulse  is  applied,  it  acts  to  decrease  the  depletion  width  and  causes 
empty  electron  traps  formerly  in  the  depletion  region  to  be  filled  (see 
Figure  3b).  Immediately  after  the  pulse,  the  depletion  region  is  deeper  than 
before  the  pulse  since  the  filled  traps  carry  the  added  negative  charge  of 
trapped  electrons.  This  decrease  in  depletion  width  corresponds  to  a 
measured  decrease  in  the  source-drain  conductance.  Because  the  time 
constants  of  deep-levels  are  relatively  long  compared  to  capacitance  time 
constants,  the  total  change  in  conductance  immediately  after  the  pulse  is  not 
difficult  to  determine.  The  more  traditional  DLTS  analysis  which  monitors 
the  rate  of  decay  of  the  trapped  carrier  charge  as  function  of  temperature  is 
used  only  independently  in  this  technique  to  determine  the  trapping  state 
energy  level.  Finally,  by  using  several  gate  pulses  of  different  magnitudes, 
information  about  the  density  of  the  deep  energy  levels  as  a  function  of 
depth  can  be  obtained. 
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Ic  should  be  noted  that  calculation  of  conductance  land  change  of 
conductance  requires  that  low-field  mobility  and  free-cirrier  distributions 
be  known.  For  this  work  mobility  was  calculated  using  the  theoretical 
results  of  Walukiewlcz  et  al.  [9]  In  conjunction  with  Monte  Carlo  [10] 
velocity-field  simulations.  The  data  obtained  in  this  way  was  then  curve  fit 
to  obtain  an  empirical  expression  for  mobility  as  a  function  of  ionized 
impurity  density  and  compensation  ratio.  The  resulting  expression  is  given 
in  reference  [11]  and  plotted  against  normalized  Walukiewlcz  values  in  Figure 
4.  The  figure  illustrates  the  excellent  agreement  between  the  two.  A 
further  check  on  the  mobility  calculation  can  be  made  by  comparing  the 
computed  mobility  with  that  obtained  from  a  magnetotransconductance 
measurement  [ 13 ] . 

MEASUREMENT * APPARATUS 


Figure  7  is  a  schematic  diagram' of  the'  MCNC/NCSU  electrical  properties 
characterization  laboratory.  All  measurements  are  controlled  by  an  LSI  11/23 
computer  which  stores  the  data  and  can  transmit  it  to  larger  mainframe 
computers  and  remote  users.  This  system  will  be  described  more  completely  at 
a  later  time. 

I— V  characteristics  of  devices  are  obtained  with  an  HP4145A 
semiconductor  parameter  analyzer  which  is  connected  to  an  automatic  probe 
station  for  on-chip  measurements.  1  MHz  C-V/GV  data  is  extracted  with  a  PAR 
410  capacitance-voltage  profiler.  This  is  also  used  for  capacitance  and 
conductance  DLTS  measurements.  Fast  A/D  converters  dlgrtize  the  DLTS 
transients  and  store  them  for  later  analysis  [13].  The  Air  Products 
Helitran,  which  is  capable  of  attaining  liquid  Helium  temperatures  (4.2*K) 
can  be  inserted  between  the  poles  of  the  laboratory  electromagnet  for 
experiments  in  magnetic  fields  such  as  mobility  and  magnetoresistance  as  a 
function  of  temperature.  The  cryostat  also  has  quartz  windows  allowing 
optical  excitation  of  the  samples. 

Since  all  measurements  are  completely  automated,  new  experiments  can  be 
performed  simply  by  writing  specialized  software  routines.  This  increases 
the  flexibility  of  the  system  by  increasing  the  number  of  measurements 
possible. 

PROFILE  SOLUTION 

Differential  capacitance  data  allows  for  the  calculation  of  free-carrier 
distribution  if  the  deep-level  distribution  is  known.  Likewise,  DLTS  data 
allows  for  the  calculation  of  deep-level  distributions  provided  the 
free-carrier  distribution  is  known.  The  combination  of  the  two  measurements 
when  coupled  with  appropriate  device  equations  provides  sufficient 
information  for  determining  the  profiles  of  interest.  References  [6]  and  [7] 
outline  a  method  for  obtaining  solutions  to  the  set  of  highly  non-linear 
integral  and  differential  equations  which  must  be  solved  simultaneously.  The 
method  involves  assuming  a  deep-level  distribution  in  the  material  and  then 
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using  Che  C-V  data  Co  compute  Che  DLTS  results.  The  accrual  DLTS  measurement 
data  can  Chen  be  compared  to  Che  predictions  for  Che  assumed  crap  density 
profile.  By  modifying  the  assumed  crap  density  profile  good  agreement 
between  Che  change  in  conductance  predicted  and  that  which  is  measured  can  be 
obtained . 

EXPERIMENTAL  RESULTS 


Measurements  were  performed  on  a  1  pm  gate  length  MESFET  and  a  fat  FET 
(250  z  200  pm)  adjacent  to  it  on  the  same  chip.  The  active  channel  was 
formed  by  implanting  silicon  into  a  Chromium  doped  GaAs  substrate. 

Conductance  DLTS  and  1  MHz  C-V  measurements  were  performed  on  these  devices 
and  a  dominant  deep-level  donor  site  was  identified  at  0.736  eV  below  the 
conduction  band.  This  traj.  is  probably  (EL2) .  The  findings  of  Martin  et  al. 
[12]  also  indicate  that  an  (EL2)  dominance  is  to  be  expected. 

The  measurement  data  was  input  into  a  computer  simulation  which 
calculates  the  desired  profiles  as  discussed  in  references  [6]  and  [7].  For 
a  first  guess  at  the  deep-level  distribution,  constant  density  profiles  were 
assumed.  The  solid  lines  of  Figure  5  show  the  resulting  conductance 
predictions  for  various  trap  densities  along  with  measured  data.  Note  that 
small  values  of  Va-Vp  correspond  to  regions  deep  in  the  active  channel  while 
larger  values  represent  areas  close  to  the  surface.  Using  simple  step 
function  approximations  for  the  trap  density  profile  results  in  considerable 
improvement  in  the  obtained  agreement. 

•m 

Figure  6  illustrates  free-carrier,  ionized  donor  and  deep-level 
concentrations  which  result  from  the  analysis.  Note  that  some  scatter  of  the 
ionized  donor  data  begins  to  occur  deep  in  the  channel.  This  begins  to  occur 
when  the  trap  density  and  ionized  donor  density  are  of  the  same  order  of 
magnitude.  The  uncertainties  in  the  exact  shallow-level  concentration  at 
this  depth  are  not  critical  to  profile  predictions.  This  is  true  since  the 
magnitude  of  all  the  quantities  of  interest  are  small  at  this  depth  when 
compared  to  their  magnitudes  near  the  Implant  peak.  Despite  the  scatter  of 
the  data,  it  is  evident  from  Figure  6  that  the  free-carrier  and  shallow-donor 
concentrations  differ  significantly.  For  lower  energy  (shallow)  implants, 
this  distinction  between  the  two  profiles  is  even  greater. 

The  low-field  mobility  profile  obtained  from  this  analysis  is  also  shown 
in  Figure  6.  The  curve  can  be  compared  to  an  experimentally  obtained 
mobility  and  is  in  good  agreement.  The  low  field  mobility  profile  in  the 
channel  is  obtained  experimentally  by  monitoring  the  magnetic  field 
dependence  of  the  1.6  pm  gate  device  transconductance  [13]. 

EFFECT  ON  DEVICE  PERFORMANCE 


The  effects  which  depth  dependent  concentrations  and  transport 
properties  have  on  device  performance  have  been  studied  through  the  use  of  a 
device  model  which  is  described  in  references  [7]  and  [LI].  The  model  is  a 
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one-dimensional  model  with  small-signal  analysis  which  Accounts  for  depth 
dependent  carrier  concentrations  and  transport  properties. 

The  I— V  predictions  of  the  model  for  a  1.6  micron  ion-implanted  device 
are  compared  in  Figure  8  with  measured  I-V  curves.  The  agreement  is 
excellent.  It  should  be  noted  that  without  the  inclusion  of  depth  dependent 
transport  due  to  traps,  this  agreement  could  not  be  obtained.  The 
deep-levels  have  a  tendency  to  "soften'*  the  pinch-off  characteristics  of  the 
device.  For  all  of  the  devices  studied  in  this  work,  this  "softening"  effect 
was  required  to  obtain  best  agreement  with  measured  I-V  characteristics. 

This  pinch-off  softening  effect  is  also  illustrated  in  Figure  9  for  one 
particular  device.  In  addition  to  the  qualitative  effect  of  Figure  9,  the 
inclusion  of  depth  dependent  transport  properties  also  has  considerable 
quantitative  effect.  Figure  10  shows  the  zero  gate  bias  drain  current  and 
the  pinch-off  potential  as  a  function  of  trapping  state  concentration.  The 
trapping  state  density  of  the  figure  ranges  between  10 1*  an^  iq16  cm~3%  over 
this  range  the  zero  gate  bias  current  varies  between  30  and  18  mA,  while  the 
pinch-off  potential  varies  between  1.62  and  1.38  volts. 

CONCLUSION 


A  modeling  technique  for  determining  free-carrler,  shallow-level  donor, 
deep-level  and  mobility  profiles  of  ion-implanted  devices  has  been 
illustrated.  Experimental  results  obtained  from  C-V,  conductance  DLTS  and 
channel  mobility  profile  measurements  on  a  1.6  (job  device  are  shown  and  have 
been  incorporated  into  the  model.  Agreement  between  the  measured  and 
calculated  device  parameters  is  extremely  good  indicating  the .importance  of 
including  basic  material  parameters  into  a  device  model. 
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Figure  1.  Typical  concentration 
profiles  for  ion-implanted  material. 
Shown  are  the  free-carrier ,  shallow- 
level  donor,  and  deep-level  crapping 
state  profiles. 


Figure  2.  Effect  of  deep-levels  on 
C-V  interpretation.  The  solid  line 
represents  the  computed  free-carrier 
profile  if  no  Crapping  staces  are 
present.  The  dashed  curves  are  re¬ 
sulting  profiles  for  various  constant 
trap  densities. 


MOBILITY  M«2/V«SEC 


88 


Figure  3.  Bend  diagrams  appropriate 
to  conductance  OLTS  measurements. 

The  figures  correspond  to  times 

a)  before  the  pulse  is  applied, 

b)  during  the  pulse,  and  c)  immedi¬ 
ately  after  the  pulse. 


Figure  4.  Low-field  mobility  as  a 
function  of  doping  and  deep-level 
compensation.  The  solid  lines  are 
computed  from  an  empirical  expression 
for  9  -  0.0,  0.2,  0.4,  0.6,  and 
0.8.  The  data  points  are  from 
normalized  theoretical  computations 
[9]. 
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Figure  5.  Change  in  conductance  as  a  function  of  pulsed  voltage.  Data 
points  are  from  conductance  DLTS  measurements.  The  solid  lines  are  model 
predictions  with  NT  -  constant.  The  dashed  line  is  predicted  for  the  three 
level  step  profile  defined  in  the  figure. 
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igure  8.  A  comparison  of  model  predicted  and  measured  I-V  characteristics 
for  a  1.6  micron  ion-implanted  device.  The  lines  are  measured 
I-V  curves.  The  data  points  are  predicted  by  the  model. 
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Figure  9.  Normalised  drain 
current  vs.  normal ized  gate 
bias  with  two  different 
trapping  state  densities. 
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Figure  10.  Zero  gate  bias 
drain  current  and  pinch-off 
potential  vs.  trapping  state 
concentration.  The  implant 
energy  is  75  keV  and  the  pea  If 
doping  density  is  2x10  cm 
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ABSTRACT 


A  new  modeling  technique  for  determining  free-carrler,  shallow-level 
donor,  trap  and  mobility  profiles  of  Ion-implanted  GaAs  MESFETs  Is  presented. 
This  technique  uses  the  results  of  C-V  and  conductance  DLTS  measurements  to 
calibrate  a  theoretical  model.  The  model  eliminates  assumptions  commonly  made 
in  other  simple  characterization  techniques.  A  sample  case  is  presented  which 
indicates  the  importance  of  this  more  rigorous  treatment. 


*This  work  was  partially  supported  by  the  Army  Research  Office,  Durham,  NC  on 
contract  DAAG29— 80-K-0080. 
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NOTATION 


a  epi-layer  thickness 
On  electron  diffusion  coefficient 
e  conduction  band  energy 
D  shallow-level  donor  energy 
y  Fermi  energy  level 

X  deep-level  trapping  state  energy  level 

v  valance  band  energy  level 


Electric  field 


Jn  total  electron  current  density 
k  Boltzmann’s  constant 
L  gate  length 


free-carrler  concentration 


Nc  conduction  band  density  of  states 

Np  ionized  impurity  concentration 

Nt  deep-level  trapping  state  concentration 

q  electronic  charge 

T  temperature 

x„  depletion  width 

Ax  change  in  depletion  width 

yw  Fermi  level  cross  point  of  the  trap  level 

Ay  change  in  the  Fermi  level  cross  point 

Z  gate  width 

dielectric  constant  of  semiconductor 

permittivity  of  free-space 
Mo  low-field  electron  mobility 

Vgi  built  in  potential  of  Schottky-barrler 
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1.  INTRODUCTION 

An  attractive  technology  for  the  development  of  GaAs  integrated  circuits 
for  microwave,  millimeter-wave  and  high-speed  logic  applications  involves  ion- 
implantation  of  conductive  regions  directly  into  high-resistivity,  bulk  grown 
semiconductor  substrates  [1,2].  The  performance  of  FETs  fabricated  in  this 
manner  depends  strongly  upon  carrier  concentration  profiles  and  transport 
characteristics  as  a  function  of  depth  into  the  device  active  layer.  The 
determination  of  these  quantities,  however,  is  complicated  by  the  presence  of 
impurity  compensation  in  the  form  of  deep-levels  within  the  channel  and  by 
free  carrier  diffusion  due  to  non-uniform  impurity  profiles. 

Generally,  C-V  techniques  that  assume  charge  neutrality  conditions 
throughout  the  undepleted  semiconductor  are  utilized  to  analyze  semiconductor 
profiles.  Also,  the  effects  of  deep-levels  are  often  neglected.  These 
assumptions  are  generally  not  valid  for  the  case  of  ion-implanted  material 
typically  used  for  GaAs  MESFETs.  When  diffusion  and  compensation  effects  must 
be  considered  the  transformation  from  bias  dependent  to  depth  dependent 
parameters  is  complicated.  In  general,  it  is  not  possible  to  accurately 
determine  the  various  profiles  of  interest  for  non-uniformly  doped 
semiconductors  solely  from  C-V  measurements. 

In  this  peper  we  present  a  modeling  and  characterization  technique  for 
determining  free  carrier,  shallow-level  donor,  deep-level  and  carrier  mobility 
profiles.  This  technique  is  applicable  for  the  analysis  of  material  with 
arbitrary  deep-level  and  doping  profiles.  This  method  combines  experimental 
differential  capacitance  (C-V)  and  conductance  deep-level  transient 
spectroscopy  (DLTS)  measurements  with  a  theoretical  characterization  modal. 

The  combination  of  the  C-V  and  DLTS  measurements  allows  the  carrier  mobility 
as  a  function  of  depth  to  be  accurately  determined  in  an  indirect  manner. 
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The  profiles  obtained  with  this  method  have  been  used  in  a  theoretical  MESFET 
device  model  to  analyze  the  effects  of  various  profiles  upon  device 
performance  _PT-  Excellent  agreement  between  the  theoretical  and  experimental 
device  performance  was  obtained.  Such  agreement  offers  indirect  evidence  of 
the  accuracy  of  the  characterization  technique. 

2.  GENERAL  DEVICE  EQUATIONS 

Figure  1  shows  typical  f ree-carrler,  background  donor,  and  deep-level 
profiles  for  a  sample  of  ion-implanted  semiconductor  material.  Typically, 
peak  doping  ranges  from  about  9xl0l6ca-3  Co  4xl0^cm~3.  The  deep-level 
concentrations  vary  with  material  growth  technology  and  generally  range  from 
low  10l5ea-3  to  a8  high  as  5xl0^c*“3.  Even  in  the  absence  of  deep  trapping 
centers,  f ree-carrler  diffusion  from  highly  doped  to  more  lowly  doped  regions 
will  cause  the  free-carrler  profile  to  differ  from  that  of  the  ionized  donors 
as  shown  in  Figure  1.  This  effect,  also  known  as  Debye  tailing,  can  be 
significant.  Our  studies  show  that  as  much  as  an  order  of  magnitude 
difference  can  exist  between  these  two  profiles  for  ion-implanted  materials 
typically  used  in  the  fabrication  of  GaAs  MESFETs. 

The  relationship  between  the  free-carrler,  donor,  and  deep- level  profiles 
can  be  obtained  by  combining  the  current  density  equation  for  electrons  with 
Poisson's  equation  in  the  direction  into  the  channel.  The  appropriate  form 
for  the  one-dimensional  current  density  equation  in  n-type  material  is 

dn 

Jn  ■  0  ■  q (pon(x)E(x)  +  Dn—-)  .  (1) 

dx 

To  write  the  charge  balance  equation  appropriate  for  bulk  semiconductor 
material,  the  deep-level  characteristics  must  be  known.  For  the  case  where 
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deep-level  acceptors  act  aa  the  dominant  electron  trap  in  the  semiconductor 

material,  Poisson's  equation  is  written  as  [4] 

il-  {ND(x)-n(x)-NT(x)  >.  (2a) 

dx  «r«b 

When  deep-level  donors  constitute  the  electron  traps  in  the  material,  equation 
(2a)  must  be  modified.  In  undepleted  semiconductor  material,  deep-level  donor 
states  will  normally  lie  below  the  Fermi  level.  Thus,  they  will  be  filled  and 
carry  no  charge.  The  appropriate  form  of  Poisson's  equation  for  this  case  is 
[41 

d£  -  q  {Ni)OO-n(x)  >.  (2b) 

dx  ^*0 

Equation  (1)  and  the  appropriate  form  of  equation  (2)  can  be  combined  to 
express  one  of  the  profiles  in  terms  of  the  other  two.  In  the  deep-acceptor 
trap  case,  the  background  donor  concentration  can  be  written  as 

kTCf^,  d  [~1  dn~| 

Nd(x)  «  n(x)  +  ^(x) - —  — -  (3a) 

q2  dxJji(x)dxJ 

where  the  Einstein  relation  has  been  used.  Likewise,  when  the  traps  consist 
mainly  of  deep-donors,  the  shallow-level  donor  concentration  can  be  written 

kT  Cg  so  d  ["*1  dnl 

ND(x)  -  n(x) - - .  (3b) 

q2  dx  Jji(x)  dxj 

The  relationships  (3a)  and  (3b)  are  identical  to  that  derived  by  Kennedy 
and  O'Brien  [5]  if  trap  centers  were  not  present.  The  last  term  in  both 
equations  (3a)  and  (3b)  represents  the  local  charge  Imbalance  that  will,  in 
general,  exist  within  the  material.  The  Importance  of  this  local  steady-state 
charge  imbalance  in  terms  of  the  device  measurements  will  be  discussed  In  more 
detail  in  the  following  sections. 
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3.  DIFFERENTIAL  CAPACITANCE  MEASUREMENT  (C-V) 

Using  Che  depletion  approximation,  Che  differential  capacitance  (C-V) 
measurement  can  be  considered  as  a  parallel-plate  capacitor  problem.  The 
capacitance,  C,  is  normally  measured  with  a  small  oscillating  voltage,  AV, 
superimposed  over  a  fixed  applied  biaa,  Va.  The  oscillating  voltage  induces  a 
current  due  to  an  incremental  charge,  AQ,  being  alternately  covered  and 
uncovered  near  the  depletion  region  edge,  xw.  Kennedy  et.  al.  [6]  have  shown 
that  in  the  absence  of  deep  states,  traditional  C-V  analysis  estimates  the 
free-carrier  concentration  since  only  the  free-carriers  can  participate  in  the 
differential  change  in  the  device  charge  AQ. 

If  deep  trapping  states  are  present,  then  the  incremental  charge,  AQ,  is 
due  to  free-carriers  in  the  conduction  band  and  electrons  trapped  in  the 
deep-level  sites.  The  interpretation  of  C-V  data  under  these  conditions  is 
dependent  on  the  relationship  between  the  frequency  of  the  oscillating 
voltage,  ui(  AV) ,  the  electron  emission  rate  associated  with  the  deep  energy 
site,  Oq,  and  the  rate  of  change  of  the  dc  bias,  &d(AVa).  Both  deep-level 
donors  and  acceptors  modulated  at  various  speeds  with  respect  to  en  have  been 
examined  by  a  number  of  workers  [7-9].  For  simplicity,  however,  this  paper 
will  examine  only  the  case  where 

u(AV)  »  »  u(  AVa). 

The  band  diagrams  for  this  case  are  shown  in  Figure  2.  These  conditions  are 
chosen  since  they  would  commonly  hold  for  C-V  measurements  on  bulk  GaAs 
substrates  ion-implanted  with  n-type  dopants. 

Kimmerling  [9]  has  shown  that  when  acceptor  type  traps  make  up  the 
deep-levels  the  free-carrier  concentration  under  the  previously  stated 
assumptions  can  be  expressed  as 

Fw  & 

n(x^f)  -  nc_v(xw)  -  »t<7w)  ""  +  l*T<*w) •  (*•) 

XwAx 
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For  donor  type  Craps  Che  corresponding  expression  is 

ywAy 

^(*y)  m  ^c— v(3Cw^~ ^T^yw^ "  •  (^b) 

*w*t 

In  equacions  (4a)  and  (4b),  Che  term  nc_v  is  given  by 

C3  /dC\*1 

Uc-vOO - ( - 1  (5) 

q  %  «o  \dV  . 

and  represencs  Che  cradicional  C-V  expression.  The  Fermi  level  cross-point  yw 
can  be  derived  from  the  band  bending  if  Che  background  donor  concentration  is 
known.  Referring  to  Figure  2a  and  applying  Poisson's  equation  gives 

Ef-€T  -  - - [nd(*w)  -  »D(yc)  -  X(xw)ND(xw)l  (6) 

J 

where  Che  notation 

N(a)  -  J*H(x)dx 

o 

and 

N(a)  ■  /S/tN(z)dxdx 
o  o 

is  adopted.  Knowledge  of  the  valueCp  -  Ef  is  also  required  to  determine  yw 
from  equation  (6).  Normally, £c  -  Cj  is  determined  through  independent 
measurement  [10,13]  and  the  relation 

£c  -  CP  -  kT*n(Nc/n(x))  (7) 

is  then  used  with  this  information  to  determine  tp  -  Zj. 

Equations  (4a),  (4b)  and  (6)  can  also  be  written  in  terms  of  the  quantity 
\(xw)  through  the  relationship 


M*w)  ■  *w-yw 


(8) 


99 


Differentiation  of  equation  (8)  with  respect  to  Xy  and  substitution  into 
equation  (4a)  results  in  the  expression 

y w  ^  ^ 

n(xw)  -  nc_v(xw)-NT(yw) - (1 - )  +  NT(xw).  (9a) 

*v  dxw 


To  obtain  this  expression  small  perturbations  are  assumed  so  that 
Likewise,  equation  (4b)  becomes 


At 

4 


3|Z. 

dx 


yw  dX. 

n(Xv)  -  nc-v(xw)-NT(yw) - (1 - ).  (9b) 

xw  dxw 

Figure  3  shows  a  profile  calculated  from  traditional  C-V  analysis  along 
with  corresponding  free-carrier  profiles  calculated  from  equations  (9a) 
and  (9b)  for  various  constant  deep-level  acceptor  and  donor  trap  densities. 

The  deep  levels  have  a  significant  effect  upon  the  shape  of  the  free  carrier 
profile,  especially  near  the  tail  of  the  implant.  Since  low  noise  devices  are 
typically  biased  so  that  the  conducting  channel  is  in  the  tail  region,  the 
deep-level  concentration  is  expected  to  have  a  significant  Influence  upon 
device  performance. 

4.  CONDUCTANCE  DLTS 

Capacitance  deep-level  transient  spectroscopy  (DLTS)  was  originally 
proposed  by  Lang  [10]  and  described  in  various  forms  by  others  [11-12]  as  a 
technique  to  characterize  traps  in  semiconductors.  Similar  information  can 
also  be  obtained  from  conductance  DLTS  [13],  and  it  is  this  method  that  has 
been  utilized  in  this  work.  The  technique  involves  monitoring  the 
source-drain  conductance  of  an  FET  structure  while  a  repetitive  gate  voltage 
pulse  is  applied.  The  band  diagrams  corresponding  to  times  before,  during  and 
after  the  gate  pulse  are  illustrated  in  Figure  2.  Before  the  pulse,  a  reverse 
bias  Va  has  the  channel  nearly  pinched  off  so  that  only  a  small  steady-state 
conductance,  go,  is  measured.  The  pulse  acts  to  decraase  the  depletion  region 


too 


and  causes  empty  electron  traps  formerly  In  the  depletion  zone  to  be  filled 
(see  Figure  2b).  Immediately  after  the  pulse,  the  depletion  region  is  deeper 
than  before  the  pulse  since  the  filled  traps  carry  the  added  negative  charge 
due  to  the  trapped  electrons.  This  increase  in  depletion  width*— signified  as 
As  in  Figure  2c— causes  a  decrease  in  the  source-drain  conductance,  The 

non-equilibrium  change  in  conductance  decays  as  the  filled  traps  emit 
electrons  back  into  the  conduction  band.  The  decay  process  is  very  slow 
compared  to  gate  capacitance  time  constants,  however,  so  that  Ag  due  to  filled 
traps  immediately  after  the  pulse  is  not  difficult  to  determine.  This  change 
in  conductance  from  the  ataady-atate  value  immediately  after  the  pulse  is  of 
interest  here.  The  actual  rate  of  decay  as  a  function  of  temperature  which 
characterizes  more  traditional  DLTS  analysis  [10-13],  is  used  independently  to 
determine  the  trapping  state  energy  level.  Finally,  by  using  several  gate 
pulses  of  different  magnitudes,  information  about  the  deep  energy  levels  as  a 
function  of  depth  can  be  obtained. 

Under  the  steady-state  condition  of  Figure  3a,  the  conductance  go  can  be 
computed  from 


So 


qZ  a 

—  / 

L  x* 


Mo(x)n(x)dx. 


(10) 


Likewise,  the  change  in  conductance  immediately  after  the  pulse  can  be  written 
from  Figure  2c  as 


qZ  xw+Ax 

Ag  - - /  Mo(x)n(x)dx.  (11) 

L  Xy 

The  conductions  quantities  expressed  by  equations  (10)  and  (11)  represent  only 
the  conductance  beneath  the  gate.  For  short  channel  devices,  the  parasitic 
source  and  drain  resistance-  wi  .  also  contribute  to  the  measured  conductance. 
For  this  work,  the  parasitic  resistances  ware  determined  and  subtracted  from 


the  measured  terminal  conductance. 
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Equations  (10)  and  (11)  require  that  mobility  be  known  as  a  function  of 
donor  density  and  background  compensation.  This  information  was  obtained 
using  the  theoretical  results  of  Walukiewlcz  et  al.  [14]  in  conjunction  with 
Monte  Carlo  velocity-field  predictions.  The  mobility  as  a  function  of 
background  donor  density  with  no  traps  present  was  determined  from  theoretical 
Monte  Carlo  calculations  [IS],  and  the  Walukiewlcz  values  were  then  normalized 
to  the  Monte  Carlo  numbers.  The  normalized  data  was  used  to  obtain  an 
empirical  expression  for  mobility  as  a  function  of  background  donor  density 
and  compensation  ratio.  The  resulting  expression  is 


where 


Mo 


Mmax 

‘log  tr* 
p 

»o 


(l-0)b 


N'  -  %  •  m^, 

D 

Umar  "  8380  (cm2/V*sec) , 


N0  -  23.2553, 
c  -  23.0, 


0  -  »t/Nd, 


(12) 


0.025* (log  N’)2  -  0.817278*(log  N')  +  6.252838 
D  D 

for  No  >  1021 


0.114992 


^  for  Nj,  <  1021  (m-3), 

and  where  Nq  is  given  in  (m~3).  Expression  (12)  is  plotted  against  the 
normalized  Walukiewlcz  values  in  Figure  4.  As  can  be  seen  from  the  plot,  the 
agreement  is  quite  good. 
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Another  restriction  which  applies  to  Che  situations  Illustrated  In 
Figure  2  Is  that  Poisson's  equation  oust  hold.  As  noted  In  section  2,  a 
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non-zero  charge  density  exists  within  the  channel  independent  of  bias  which 
gives  rise  to  a  potential  gradient.  This  potential  is  not  the  quantity  of 
Interest  for  this  analysis.  Rather,  it  is  the  relationship  of  the  change  in 
applied  potential  to  the  depletion  width  which  oust  be  described.  These 
quantities  are  related  through  the  charge  which  can  be  moved — n(x)  and 
NT(x) — and  are  therefore  independent  of  Nq(x) .  Thus,  the  appropriate  form  of 
Poisson's  equation  to  describe  Figure  2a  is 

Va  *  -  {n^)  -  (NT(xw)  -  NT(yw)  -  xwn(xw)  +  X(xw)(NT(xw)  -  NT(yw))}  -  VBI 

«r«b 

(13a) 

for  deep  acceptors,  and 

ya  -  —1—  {n(x^)  +  NT(yw)  *  *w*(*w>  "  Fw^Fw)  >  “  VBI  (13b) 

=r«b 

for  deep  donors.  The  applicable  equation  for  Figure  2b  is  obtained  by 
changing  Va  to  Vp,  xw  to  Xp  and  X^Xy)  to  X.(xp)  in  the  appropriate  form  of 
equation  (13).  Similarly,  from  Figure  2c  for  the  case  of  deep  acceptors 

Vp  *  -  (n(xw+Ax)  -  (Nx(xy+Ax)  -  NT(yp))  -  (xw+Ax)n(xw+Ax)  +  (Xv+Ax-yp) 

®r 

•(Nx(x„+Ax)  -  Nx(yp))>  -VBi,  (14a) 

while  the  case  for  deep  donors  is  described  by 


<1  -  —  _  _ 

Vp  ■  -  (n(xw+Ax)  -  Nx(7p)  “  (xv+As)n(xw+Az)  -  ypNT(yp) }  -  VBI.  (14b) 

«r«b 
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Ic  can  be  noted  from  Figure  2  that  only  the  trap  sites  within  a  distance 
of  Sv-Xlxt,)  from  the  surface  can  alternately  be  filled  and  emptied  by  the 
pulse  Vp.  For  the  case  studied  here,  this  corresponds  to  a  depth  of  about  0.1 
micron.  For  distances  into  the  device  greater  than  this,  the  trap 
concentration  enters  into  the  computations  though  the  free-carrier 
concentration  as  specified  by  equation  (9),  and  through  carrier  mobility  as 
expressed  in  equation  (12). 

5.  METHOD  OF  S0LPTI0N 

The  differential  capacitance  data,  conductance  DLTS  data,  and  measurement 
of  the  trapping  state  energy  level  when  coupled  with  the  equations  of  the 
preceding  sections  provide  sufficient  information  for  determining  the  profiles 
of  Interest.  Certain  problems  exist,  however,  in  utilizing  the  data  and 
equations  together. 

Equation  (10)  requires  that  a  value  for  the  maximum  epllayer  depth  be 
known.  With  ion-implanted  profiles,  however,  this  is  a  difficult  quantity  to 
determine.  For  the  low-energy  (70  -  140  keV)  Si  implants  studied  here,  it  was 
found  that  a  good  value  for  epi-thlckness  was  on  the  order  of  0.3  to  0.35 
micron.  This  corresponds  to  a  depth  into  the  material  where  free-carrier 
density  reaches  a  level  a  little  more  than  two  orders  of  magnitude  below  the 
peak  free-carrier  density  in  the  structure. 

The  limitation  of  the  C-V  measurement  technique  to  obtain  data  near  the 
surface  of  the  sample  creates  another  difficulty  in  the  analysis.  Equations 
(3),  (9),  (13)  and  (14)  all  require  knowledge  of  the  free-carrier 
concentration  throughout  the  active  layer.  This  difficulty  is  overcome  by 
requiring  that  equations  (10)  and  (13)  hold  simultaneously.  In  fact,  if  the 


L 


C-V  data  were  available  throughout  the  active  layer,  then  either  equation  (10) 
or  (13)  would  be  redundant  and  the  problem  would  be  overspecified.  In 
collecting  conductance  DLTS  data,  both  the  reverse  bias,  Va,  and  the 
steady-state  conductance,  go,  are  measured.  Also,  from  C-V  analysis  and 
equation  (4),  a  value  for  free-carrier  concentration  is  obtained  from  some 
minimum  obtainable  depletion  width,  x^q,  to  the  maximum  epi-layer  extent,  a. 
For  values  of  x  between  the  material  surface  and  x^q,  the  free-carrier 
profile  can  only  be  estimated.  The  initial  estimate  of  the  free-carrier 
concentration  near  the  surface  can  be  made  by  assuming  that  the  profile  should 
be  approximately  gausslan  in  this  region  with  the  slope  of  the  profile  flat  at 
the  surface.  Once  a  profile  has  been  established  over  the  entire  active 
region,  equation  (13)  can  be  used  to  calculate  the  depletion  width,  xw, 
corresponding  to  applied  bias,  Va.  Likewise,  equation  (10)  can  be  used  to 
compute  x„  for  a  steady-state  conductance,  go .  Any  discrepancy  between  the 
two  values  for  xw  obtained  by  these  equations  must  be  due  to  inaccuracies  in 
the  built-in  potential  used  or  in  the  free-carrier  concentration  estimate 
between  the  surface  and  x^  n .  Adjustments  can  then  be  made  to  either  of  these 
quantities  until  the  two  values  for  depletion  width  agree. 

It  should  be  noted  that  this  method  produces  only  an  equivalent  built-in 
potential/  free-carrier  profile  pair  near  the  surface.  However,  since  under 
normal  device  operation  of  an  ion-implanted  MESFET  this  region  would  always  be 
depleted,  this  information  is  usually  adequate. 

To  solve  the  equations  of  the  preceding  sections  simultaneously  in  a 
rigorous  fashion  is  difficult.  Equations  (11)  and  (14)  must  hold  for  each 
pulse  used  in  the  conductance  DLTS  measurement.  Thus,  if  twenty  different 
pulse  heights  were  used ,  there  would  be  over  forty  integral  and  differential 
equations  that  must  be  solved  simultaneously. 
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An  alternative  to  this  method  is  to  assume  that  the  background  trap 
profile  is  known,  and  to  compute  the  resulting  conductance  OLTS  data*  The 
actual  DLTS  measurement  data  can  then  be  compared  to  the  predictions  for  the 
assumed  trap  profile.  By  making  modifications  on  the  assumed  trap  profile 
good  agreement  between  the  change  in  conductance  predicted  and  that  which  is 
measured  can  be  obtained. 

This  latter  method  was  chosen  for  the  work  presented  here.  The  system  of 
equations  was  solved  iteratively:  each  equation  being  solved  individually  to 
update  only  one  quantity  at  a  time.  To  begin,  the  trap  concentration  was 
assumed  equal  to  some  constant  value  throughout  the  device.  The  depletion 
width,  la,  corresponding  to  the  steady-state  bias,  Va ,  was  determined  using 
equations  (10)  and  (13)  as  outlined  above.  Equation  (3)  was  used  to  update 
the  shallow-level  donor  profile.  The  function  W*v)  was  then  obtained  from 
equation  (6)  and  the  free-carrler  profile  determined  through  the  use  of 
equation  (9).  Normally,  three  to  six  iterations  through  the  equations  were 
required  before  convergence  was  achieved. 

Once  the  free-carrier  and  shallow-level  donor  profiles  were  determined, 
equation  (14)  was  used  to  compute  the  depletion  width  change.  Ax,  as  a 
function  of  pulse  bias,  Vp.  This  calculation  requires  the  mobility  to  be 
determined  as  a  function  of  depth  through  use  of  equation  (12).  This 
information  is  then  easily  utilized  with  equation  (11)  to  predict  conductance 
DLTS  data.  Finally,  the  results  were  compared  to  the  DLTS  data  actually 
measured  on  the  device.  At  this  point,  adjustments  could  be  made  to  the  trap 
profile  originally  assumed,  and  the  process  repeated  until  good  agreement  was 


obtained 
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6.  RESULTS 

A  1.6  micron  gate  length  MESFET  and  a  differential  capacitance  test 
pattern  were  fabricated  on  a  Silicon  implanted  GaAs  substrate  [16]. 
Conductance  DLTS  and  C-V  measurements  were  performed  on  these  devices  and  a 
dominant  deep-level  donor  state  (electron  trap)  was  identified  at  0.736  eV 
below  the  conduction  band.  Although  positive  identification  of  this  trap  is 
difficult,  the  level  is  probably  EL2  [17,18].  The  dominance  of  this  level 
agrees  with  the  findings  of  Martin  et.  al.  [19].  The  C-V  and  DLTS  data  was 


used  as  input  to  a  computer  simulation  which  calculates  the  desired  profiles 
as  outlined  in  the  preceding  sections. 

As  a  first  guess,  the  trap  concentration  was  assumed  constant.  The  solid 
lines  of  Figure  5  give  the  resulting  conductance  DLTS  predictions  for  various 
trap  densities  along  with  the  measured  data.  Notice  that  small  values  of 
Va-V p  represent  areas  deep  in  the  device  while  larger  values  represent  areas 
close  to  the  surface.  The  use  of  simple  step  function  approximations  for  the 
trap  profile  results  in  improved  agreement  between  the  model  and  experimental 
results.  A  step-profile  defined  by 


^*3.0x10*5  cm ”5  o.O  <  x  <  0.040  micron 

%  *  <  2.5x10*5  cm“3  0.04  <  x  <  0.055  micron 

U. 0x10*5  cm“3  0.055  <  x  <  a 


was  used  to  obtain  the  results  shown  by  the  dashed  line  of  Figure  5.  Here  the 
agreement  between  model  predictions  and  measurements  is  excellent. 

The  resulting  f ree-carrler ,  shallow-level  donor,  and  deep-level  trap 
concentrations  as  a  function  of  depth  into  the  material  are  shown  in  Figure  6. 
Notice  that  deep  into  the  channel  there  is  some  scatter  of  the  shallow-level 
donor  data.  This  begins  to  occur  when  the  trap  concentration  and  the 
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shallow-level  donor  concentrations  are  of  the  same  order  of  magnitude.  The 
uncertainties  in  the  exact  shallow-level  concentration  at  this  depth  into  the 
channel  are  not  critical  to  the  profile  predictions*  This  is  true  since  the 
aagnitude  of  all  the  quantities  of  interest  are  small  at  this  depth  when 
compared  to  their  magnitudes  near  the  implantation  peak.  Despite  the  scatter 
of  the  data,  it  is  evident  from  Figure  6  that  diffusion  of  the  electrons  has 
affected  free-carrier  and  shallow-level  donor  concentrations.  For  lower 
energy  (shallow)  implants,  this  distinction  between  the  two  profiles  is  even 
greater. 

The  low-field  mobility  profile  obtained  from  this  analysis  is  also  shown 
in  Figure  6.  It  should  be  emphasized  that  the  mobility  versus  depth  profile 
is  determined  indirectly  with  the  use  of  equation  (12).  Since  this  technique 
requires  only  the  C-V  and  DLTS  measurements  it  allows  the  various  profiles  to 
be  determined  in  an  efficient  manner. 

The  mobility  profile  shown  in  Figure  6  is  in  good  agreement  with  the 
results  presented  by  Das  and  Kim  [20],  Hobgood  et  al.  [21]  and  Jay  and  Halils 
[22].  In  addition,  this  technique  has  been  used  to  determine  profiles  for  use 
in  a  theoretical  device  model  to  study  the  effects  of  the  various  profiles 
upon  device  performance  [3].  Excellent  agreement  between  the  theoretical 
predictions  and  experimental  device  performance  was  obtained  [3]. 

From  Figure  4  it  should  be  noted  chat  when  the  deep-level  density  is  less 
than  about  10l5cm-3>  the  C-V  data  can  be  Interpreted  in  the  traditional  manner 
and  the  Ny(x)  terms  in  the  equations  can  be  neglected.  Under  these  special 
conditions  n(x)  is  also  known  and  the  entire  analysis  is  slmplfled.  Even  for 
this  special  case,  however,  the  shallow-level  donor  and  free  electron 
concentrations  wet  be  treated  as  separate  functions  in  order  to  accurately 
determine  the  carrier  mobility  profile. 
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7.  CONCLUSIONS 

A  modeling  technique  for  determining  the  free-carrler,  shallow-level 
donor,  deep-level  and  mobility  profiles  of  ion-implanted  GaAs  material 
suitable  for  MESFET  fabrication  has  been  presented.  The  technique  involves 
non-destructive  electrical  C-V  and  conductance  DLTS  measurements  coupled  with 
the  solution  of  the  applicable  device  equations.  The  appropriate  equations 
have  been  developed  and  a  solution  method  presented.  Parameter  profiles  as  a 
function  of  depth  suitable  for  device  characterization  and  modeling  are  easily 
determined.  Deep-levels  and  free-carrler  diffusion  are  shown  to  have 
significant  effects  upon  free-carrler  transport  characteristics.  The  use  of 
this  technique  on  a  sample  case  has  been  presented. 
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Figure  5 


Typical  concentration  prof ilea  for  ion-implanted  Material*  Shown 
are  Che  free-carrler,  shallow-level  donor,  and  deep-level  crapping 
scace  profiles.  The  disclncdon  beCween  free-carrler  and 
shallow-level  donor  concencradons  will  exisc  even  wichouc  Che 
presence  of  deep-levels. 

Band  diagrams  appropriate  Co  conductance  DLTS  measurements •  The 
flguree  correspond  to  tines  a)  before  Che  pulse  ie  applied,  b) 
during  Che  pulse,  and  c)  limed lately  after  Che  pulse. 

Effect  of  deep-levels  on  C-V  interpretation.  The  solid  line 
represents  Che  conpuCed  free-carrler  profile  if  no  crapping  states 
are  presenc.  The  dashed  curves  are  resulting  profiles  for  various 
constant  trap  densities. 

Low-field  nobility  as  a  function  of  doping  and  deep-level 
compensation.  The  solid  lines  are  computed  from  equation  (12)  for 
0  “  0.0,  0.2,  0.4,  0.6,  and  0.8.  The  data  points  are  from 
normalized  theoretical  computations  [ 14 ) . 

H 

Change  in  conductance  as  a  function  of  pulsed  voltage.  Data  points 
are  from  conductance  DLTS  measurements.  The  solid  lines  are  model 
predictions  with  %  “  constant.  The  dashed  line  is  predicted  for 
the  three  level  step  profile  defined  in  the  figure. 


Figure  6  Resulting  concentration  profiles  and  low-field  mobility  profiles. 
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A  New  Method  of  Tolerance  Enhancement  for 
Filters  and  Amplifier  Matching  Networks* 

A.  N.  Riddle  and  R.  J.  Trew 
Department  of  Electrical  and  Computer  Engineering 
North  Carolina  State  University 
Raleigh,  NC  27650 
(919)  737-2336 

Abstract 

A  new  method  for  increasing  the  tolerance  of  passive  networks  to  complex 
load  variations  is  presented.  The  method  is  based  upon  utilization  of  the 
angular  relationship  between  a  network’s  reflection  coefficient  and  the 
corresponding  changes  with  component  variations.  A  new  filter  polynomial  that 
produces  a  response  which  has  a  greater  tolerance  to  lo^s,  component 
variations  and  load  reactance  variations  than  Butterworth,  Chebys'hev  or 
Elliptic  structures  has  been  developed  using  this  method.  Examples  using  this 
new  filter  for  tolerance  enhancement  of  matching  filters  and  FET  feedback 
amplifiers  are  presented. 


*Thls  work  was  supported  by  the  Army  Research  Office,  Durham  NC  on  contract 
DAAG29-80-K-0080  and  by  the  Office'  of  Naval  Research  through  an  ONR  Fellowship 
to  Mr.  Riddle. 
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Introduction 

One  method  for  the  development  of  low-cost  microwave  Integrated  circuits 
involves  the  design  of  passive  networks  that  are  tolerant  of  component  and 
load  impedance  variations.  Such  networks  when  used,  for  example,  in  amplifier 
circuits  would  produce  essentially  uniform  terminal  characteristics  while 
allowing  for  typical  device— to-device  parameter  variations.  While  this 
technique  has  considerable  potential  for  use  In  designing  low  cost  MICs,  very 
little  work  has  been  reported  on  appropriate  design  procedures.  A  new  method 
for  determining  suitable  matching  networks  is  presented  in  this  paper.  In 
particular,  a  solution  to  the  approximation  problem  for  deriving  filters  with 
greater  tolerances  to  load  parasitica,  loss  and  element  tolerances  is 
presented.  A  polynomial  la  developed  for  synthesizing  a  network  with  less 
sensitivity  to  parameter  variations  than  Butterworth,  Chebyshev  or  Elliptic 
structures.  Examples  of  circuits  utilizing  this  new  filter  polynomial  are 
presented  and  compared  with  conventional  filter  designs'.  The  method  presented 
in  this  paper  is  Intended  to  be  used  with  the  synthesis  procedures  discussed 
elsewhere  [1,2].  Only  lumped  element  prototypes  are  considered  so  that  the 
effects  of  the  filtering  function,  rather  than  a  particular  realization 
method,,  may  be  studied. 

The  first  detailed  consideration  of  ideal  responses  for  matching  filters 
was  presented  by  Fano  [3].  Although  Fano  demonstrated  that  a  low  ripple 
Chebyshev  response  approximating  a  constant  mismatch  was  superior  to  a  large 
ripple  Chebyshev  response  with  the  same  peak  mismatch,  he  did  not  explore 
filter,  responses  other  than  Butterworth,  Chebyshev  and  Elliptic  types.  The 
filter  polynomial  developed  in  this  paper  results  in  a  response  which 
maximizes  the  tolerance  of  the  filter  to  load  reactance  variations  such  as  in 
the  input  capacity  of  an  FET.  This  response  is  derived  by  realizing  that  the 
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tolerance  of  e  filter  to  an  eleaent  variation  is  maximized  when  that  variation 
causes  a  change  in  Su(f)  which  is  perpendicular  to  the  direction  of  Sn(f) 
[4].  Therefore,  a  filtering  polynomial  is  developed  which  causes  the  change 
in  Sii(f)  with  load  reactance  variations  to  be  perpendicular  to  the  direction 
of  Sn(f)  at  the  peaks  of  Su(f).  When  the  angle  factor  (defined  as  the  angle 
between  the  direction  of  Sj_i(f)  and  the  change  in  Sxi(f)  due  to  a  variation  in 
any  component)  is  an  odd  multiple  of  90  degrees  the  response  is  maximally 
tolerant  to  variations  in  that  component  [4].  Consideration  of  the  angle 
factor  leads  to  a  filter  response  which  has  a  relatively  flat  mismatch  over 
most  of  the  band  and  a  small  dip  in  Sj^(f)  at  the  band  edge.  Such  a 
characteristic  has  improved  performance  over  traditional  filters. 

Theory 

The  transformation  from  the  load  ialttance  plane  to  the  input  reflection 
coefficient  plane  and  the  reverse  process  are  both  bilinear  transformations. 
Therefore,  since  changes  in  the  real  and  imaginary  parts  of  the  load  are 
orthogonal,  resulting  changes  in  the  input  reflection  coefficient  will  be 
orthogonal.  This  property  can  be  exploited  to  increase  the  tolerance  of 
filter  networks  to  reactive  variations  in  the  load.  If  the  filter  network  is 
characterized  with  an  ABCD  matrix  and  a  load  impedance  of  the  form 

ZL(f)  -  R  +  j  X  <f) 

is  assumed,  the  change  in  a  reflection  coefficient  due  to  load  reactance 
variations  can  be  written  as 


w- 


ill 


(1) 


If  the  load  is  described  as  an  admittance  the  corresponding  equation  is 

W- 

where  TL(f)  -  G  +  JB(f). 


120 

Since  equations  (1)  and  (2)  define  a  derivative  from  only  one  circuit 
evaluation,  they  provide  efficient  tolerance  analysis. 

Figure  1  shows  the  center  of  the  reflection  coefficient  plane  with  a 
circle  enscriblng  a  maxima  tolerable  reflection,  rM-,.  for  a  circuit  such  as 
an  amplifier  input  matching  network.  A  point,  Sn(f0),  is  located  inside  this 
circle  and  represents  a  point  on  the  curve  of  input  reflection  versus 
frequency .  Increases  and  decreases  in  the  device  reactance,  typically  the 
input  capacitance,  will  move  the  point  at  Su(f0).  This  differential  movement 
is  given  exactly  by  equation  (1)  for  any  network  described  by  S-parameters 
with  load  impedance  Z^. 

The  alternating  line  of  Figure  I  shows  the  maximum  symmetric  tolerance 
about  the  point  at  Sn(f0).  The  dashed  line  of  Figure  1  shows  an  asymmetric 
tolerance  which  for  practical  purposes  should  be  limited  to  its 
symmetric  variation.  This  shows  how  variations  in  the  device  input 
capacitance  will  have  a  minimum  tolerance  if  they  cause  radial  excursions  in 
SU(f0)  and  a  maximum  tolerance  if  they  cause  excursions  of  Su(f0)  along  the 
alternating  line  of  Figure  1.  The  tolerance  in  the  device  input  capacitance 
would  be  that  percentage  change  which  caused  Sn(f0)  to  lie  on  the  I*,,,,., 
circle. 


The  change  in  Sn(f)  with  frequency  is  given  by  the  equation. 


where  A,  B,  C,  and  I)  are  the  ABCD  matrix  parameters  for  the  matching  network. 
Equation  (3.)  shows  that  for  all  low  Q  or  resistive  matching  networks  the 
absolute  direction  of  3Su(l)/81  i>  exactly  the  same  as  that  of  5Sn(f )/ax(f ), 
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depending  only  on  Che  sign  of  OX/ Of.  For  these  networks  reective  variations 
in  the  load  have  their  maximum  tolerance  when  the  magnitude  of  Su(f)  is  a 
constant  over  the  passband,  i.e.  the  matching  filter  has  a  flat  mismatched 
response.  While  practical  networks  are  not  resistive  and  rarely  low  Q, 
equation  (3)  leads  to  the  conclusion  chat  matching  networks  should  provide  a 
flat  mismatch  at  least  over  the  center  of  the  band  and  should  not  have  large 
ripples.  For  this  reason  a  Butterworth  filter  has  a  midband  tolerance  to  load 
reactance  variations  greater  than  a  Chebyshev  filter  even  when  the  Butterworth 
response  has  a  flat  mismatch  greater  than  the  peak  of  the  Chebyshev  (Figure  2 
and  Table  II).  The  superior  performance  of  the  Butterworth  response  in  its 
midband  is  due  to  its  90  degree  angle  factor. 

A  classical  approximation  for  the  tolerance  of  a  filter  to  load 
variations  is  given  by  [5,6] 

Z  classical 
Sll 

where  Sg  is  the  classical  sensitivity  defined  as: 


Note  that  equation  (1)  may  be  used  to  compute  the  classical  sensitivity 
without  the  need  for  a  derivative. 

This  approximation  gives  the  minimum  bound  on  the  load  reactance 
tolerance  as  the  tolerance  goes  to  zero.  A  better  approximation  to  the  load 
reactance  tolerance  uses  the  angular  relationship  between  Sn(f)  and  its 
change  with  respect  to  a  load  reactance  to  give 

^min 

Z  corrected  tolerance  (f)  *  100  •  |sn(f)|  | (5) 


tolerance  (f)  “ 


fas  -  , 

|Sii(f)|  I  Vu| 
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where 

dain  »  minimum  {S^  co»  (4  -  9)  i  -  [Sn  sin  ( ♦  -  9) P  > 

«nd  SU  -  |su(f)|  ,  ♦  -  <  SU<f),  9  -  <  . 

•  The  difference  between  9  end  $  it  the  angle  factor,  and  maximum  tolerances 
occur  at  odd  multiples  of  90* .  This  tolerance  measure  will  always  be  greater 
than  or  equal  to  the  classical  first  order  tolerance  approximation  because  the 
classical  tolerance  uses  the  minimum  radial  distance  from  Su(f0)  to  the  rHX 
circle,  as  seen  by  the  dashed  line  in  Figure  1,  while  the  corrected  tolerance 
includes  the  angle  of  response  variation.  It  will  be  demonstrated  in  the 
results  section  that  this  corrected  tolerance  gives  a  better  approximation  to 
the  tolerance  magnitude  and  the  frequency  where  the  circuit  first  falls  out  of 
specification. 

Because  the  Butterworth  filter  has  optimal  midband  tolerances,  and  poor 
tolerances  at  the  band  edges,  a  more  tolerant  filter  can  be  derived  from  a 

j 

response  which  has  a  wide,  flat  mismatch  at  the  midband  and  small  dips  in 
Sll(f)  coward  the  band  edges.  It  is  important  that  the  dips  in  Su(f)  do  not 
go  to  zero  so  that  they  cause  only  a  small  ripple  in  the  response.  An 
increase  in  band  edge  tolerances  results.  This  new  response  appears  as  a 
mismatched  second  order  Chebyshev  filter.  The  filtering  polynomial,  however, 
is  not  a  Chebyshev  polynomial.  The  new  filter  polynomial  has  two  ripples  for 
all  orders  and  exists  only  for  third  and  higher  order  filters.  A  Chebyshev 
polynomial,  however,  may  be  considered  as  a  degenerate  case  for  a  second  order 
filter.  The  filtering  polynomial  has  a  larger  tolerance  to  load  variations 
over  its  passband  than  a  Butterworth,  Chebyshev,  or  a  mismatched  version  of 
either.  For  a  third  order  filter  the  polynomial  is: 


r3«Cp)  ■c£P3  +  (£P2+0P+1» 


(6) 
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where  a  is  the  distance  to  the  Su(f)  zeros  (matched  esse)  in  the  normalized 
filter, 

0  - 

end  p  •  «  +  jw.  The  characteristics  of  this  filter  are: 

Faa(0)  -  1.0, 

^nctCJ®)  "  0.0, 

and  FaoKJ)  “  1*0  for  the  matched  case. 

This  may  be  used  to  construct  a  fourth  order  filter  with  the  same  a  as  a 
third  order,  and  higher  orders  may  be  obtained  from  the  recursion  formula: 

*(n+l)a  *  <P2  +  l)F(n-l)a  +  P*na.  (8) 

where  a  is  the  order  of  the  filter* 

Results 

The  network  of  Figure  3  was  used  to  examine  the  tolerances  of 
Buttervorth,  Chebyshev  and  the  new  filter  to  variations  in  a  series  load 
capacitance.  In  order  to  give  a  fair  comparison  between  the  filter  types, 
each  filter  is  designed  to  have  a  bandwidth  such  that 

I  ‘  w 

where  K  is  the  same  for  all  filters.  All  filters,  therefore,  have  the  same 
ability  to  absorb  complex  load  reactances. 

The  maximum  percentage  tolerance  of  a  series  or  parallel  load  capacitance 
assuming  that  Sn(f0)  -  0,  the  matching  network  is  lossless,  the  network  is 
reciprocal  and  the  percentage  change  in  the  capacitance  is  small  is  calculated 

*A  fourth  order  polynomial  is: 

*4«<p>  *  fr  ♦  irjr  p3  +  p2  +  p  +  1 


fe-i 

T-t?- 
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with  the  expression, 

l  tolerance  -  M  .  ^~T  ■  (10) 

where  Qf  is  the  device  input  Q  at  £0.  This  expression  establishes  an  upper 
limit  for  circuit  tolerancing  and  aids  in  device  selection. 

In  the  network  shown  in  Figure  3  the  load  is  a  narrowband  model  of  the 
input  impedance  of  a  250 wide  FET  used  over  the  3.7  to  4.2  GHz  range.  Lf ,  Cf 
and  R  are  1.4  nH,  .3  pF  and  8  ohms  respectively.  Note  that  ideal  responses 
were  obtained  by  neutralizing  the  load  reactances  with  negative  components  and 
using  an  ideal  transformer.  In  practice,  Youla's  broadband  equalization 
theory  [7]  would  be  used  to  neutralize  the  load  reactances  and  the  ideal 
transformer  could  be  absorbed  into  Che  filter  by  Norton's  transformation  or 
the  use  cf  lmittance  inverters.  Figures  3  and  4  show  plots  of'  third  order 
Chebyshev  and  Buttervorth  filters  with  a  normal  and  2Z  high  FET  capacitance. 
These  plots  are  only  from  3.75  to  4  GHz  for  clarity,  and  have  the  normal 
capacitance,  .3  pF,  as  the  circled  points  and  the  high  value  as  triangular 
points.  These  filters  were  designed  for  peak  reflections  less  than  .26  from 
3.7  to  4.2  GHz.  Note  how  the  Chebyshev  filter  violates  its  specifications  in 
the  midband  while  the  mismatched  Buttervorth  filter  fails  at  the  band  edges. 
The  Buttervorth  filter  also  had  a  higher  K  than  the  Chebyshev  which  shows  that 
for  a  given  value  of  K  a  Chebyshev  response  is  a  better  choice  than  a 
Buttervorth  response. 

The  Chebyshev  filter  is  less  tolerant  because  the  reactance  variation 
causes  a  change  in  Sn(f)  that  is  not  quite  tangent  to  Sn(f)  at  the  ripple 
peaks  (the  angle  factor  is  not  90  degrees),  and  the  change  in  Su(f)  due  to 
reactive  load  variations  becomes  tangent  to  the  curve  of  Su(f)  at  higher 
frequencies.  This  means  that  the  rapid  magnitude  change  in  Su(f)  allows 
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reactive  variations  in  the  load  to  cause  nearly  radial  variations  in  Sn(f  ) 
thereby  degrading  the  tolerance.  Chebyshev  (and  Elliptic)  responses  are  only 
optimal  when  both  passband  and  stopband  responses  have  tight  specifications. 

As  can  be  seen  fron  Figures  4  and  5,  the  mismatched  Butterworth  response  has  a 
greater  tolerance  in  the  midband  even  though  its  smallest  reflection  is  equal 
to  the  peak  reflection  of  the  Chebyshev.  The  Butterworth  filter's  tolerance 
is  reduced  at  the  band  edges  because  variations  in  the  load  reactance  cause 
changes  in  Sii(f)  which  begin  to  be  tangent  to  the  Su(f)  vector.  Since  the 
Butterworth  filter  has  nearly  ideal  tolerance  properties  in  its  midband, 
modifying  its  band  edge  response  will  lead  toward  optimally  tolerant  filters. 

As  shown  in  equation  (1),  the  design  of  load  tolerant  filters  in  lossless 
reciprocal  structures  Involves  only  the  manipulation  of  the  S2i(f)  vector. 
Reciprocal  structures  must  rely  on  angular  relationships  or  loss  for  tolerant 
design,  as  seen  via  the  S2i(f)  and  Si2<i)  terms  in  equations  (1)  and  (3).  The 
Si2(f)  term  in  equation  (1)  shows  how  circulators  and  isolators  allow 
tolerance  to  load  variations  in  nonreciprocal  networks.  A  second  point  of 
clarification  concerns  the  straight  lines  drawn  in  Figure  1.  The  bilinear 
nature  of  networks  causes  most  variations  to  have  a  finite  curvature.  This 
causes  the  calculations  in  equations  (4),  (5)  and  (10),  which  are  based  on  the 
straight  lines  of  Figure  1,  to  be  approximate. 

Varloua  bandpass  networks  represented  by  Figure  3  were  characterized  by 
their  tolerances  to  load  reactance  variations,  element  variations  and  loss 
over  the  3.7  to  4.2  GHz  range.  Table  I  gives  the  element  values  for  a 
Chebyshev  filter  of  .5  GHz  bandwidth  and  ■  .26,  a  mismatched  Chebyshev 

filter  of  .5617  GHz  bandwidth,  pux  -  .23  and  p^m  “  >115,  and  a  mismatched 
version  of  the  filter  of  equation  (6)  with  .53  GHz  bandwidth,  a  *  .78,  * 

.23  and  -  .115.  All  of  the  filters  have  a  center  frequency  of  3.942  GHz. 
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▲11  of  the  above  filters  also  have  a  K,  as  dsfioad  by  aquation  (10),  of 
approximately  7.4  x  10^.  This  was  done  to  ensure  a  physical  basis  of 
coaparison  since  K  relates  directly  to  the  load  reactance  which  nay  be 
absorbed  by  the  filter  [1].  For  lossy  filters,  the  conductances  and 
resistances  of  Table  I  appear  in  parallel  or  series  with  their  respective 
resonant  circuits  shown  in  Figure  3.  ▲  large  bandwidth  Butterworth  filter 

with  large  ml  snatch,  f%in  -  .278,  was  included  to  show  how  it  is  the 
coabinatlon  of  the  angle  factor  and  the  reflection  coefficient  magnitude  that 
deteralne  the  tolerance  to  load  variations.  As  shown  in  Table  II,  this 
Butterworth  filter  has  a  larger  tolerance  than  the  Chebyshev  filter,  yet  its 
reflection  coefficient  is  always  greater  than  the  Chebyshev.  However,  the 
Butterworth  filter  shown  here  is  not  couparabla  to  the  other  filters,  because 
of  its  bandwidth  of  2  GHs  which  makes  its  K  far  too  large  for  comparison  (see 
Figure  2). 

Table  II  gives  the  tolerances  of  the  various  filter  types  to  load 
reactance  variations.  The  additional  tolerances  in  Table  II  are  for  the  first 
and  second  least  tolerant  filter  elenents,  as  labeled.  The  tolerance 
information  is  given  for  a  maximum  $n(f)  -  .333  -  rM«  with  only  one 
component  being  varied.  This  table  shows  that  the  mismatched  Chebyshev  and 
the  new  filter  have  up  to  a  fifty  percent  Increase  in  device  tolerance  and  a 
three  hundred  percent  increase  in  element  tolerances  over  a  Chebyshev  filter. 
Furthermore,  the  new  filter  shows  a  greater  Increase  in  device  tolerance  end 
most  sensitive  element  tolerance  than  the  mismatched  Chebyshev.  These 
Increases  in  tolerance  can  neks  circuits  much  easier  to  tune  and  even 
eliminate  tuning  in  some  cases.  The  mismatched  Butterworth  response  is 
tabulated  only  to  illustrate  the  effect  of  the  angle  factor  on  circuit 
performance.  This  wide  bandwidth  filter  has  a  nearly  constant  reflection 
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coefficient  of  p  ■  0.278  across  the  3.7  to  4.2  GHz  band ,  and  an  angle  factor 
around  90  degrees.  Although  the  reflection  coefficient  of  this  Buttervorth 
response  is  higher  than  the  peak  reflection  coefficient  of  the  Chebyshev 
response,  the  Buttervorth  filter  has  a  larger  tolerance  to  load  variations. 

The  larger  tolerance  is  due  to  the  90  degree  angle  factor  of  the  broadband 
Buttervort.h  filter,  and  so  demonstrates  the  Importance  of  the  angle  factor  in 
tolerant  filter  design.  The  improvement  of  the  mismatched  Chebyshev  response 
over  the  normal  Chebyshev  response  may  be  attributed  to  the  reduction  of  its 
maximum  reflection  coefficient,  while  the  improvement  of  the  new  filter  over 
the  mismatched  Chebyshev  is  because  the  new  filter  has  an  angle  factor  nearer 
to  90*  at  the  peak  of  Sii(f).  Thus,  the  new  filter  is  the  most  tolerant 
design.  The  ideal  tolerance  of  equation  (9)  is  included  in  the  table  as  an 
approximate  upper  limit  for  the  device  tolerance. 

Table  III  compares  the  tolerances  of  the  new  filter  to  load  reactance 
variations  using  the  first  order  classical  [equation  (4)]  and  angle  factor 
corrected  [equation  (5)]  methods.  The  classical  method  gives  a  minimum  device 
tolerance  of  1.3Z  at  3.94  GHz ,  while  the  corrected  method  gives  a  minimum 
device  tolerance  of  1.7Z  at  3.85  and  4.05  GHz .  The  actual  device  tolerances 
are  1.83Z  at  3.85  GHz  and  -1.8Z  at  3.7  and  4.05  GHz  simultaneously.  This 
demonstrates  that  the  corrected  response  gives  a  better  measure  of  the  device 
tolerance  and  the  frequency  where  the  specification  is  violated. 

It  is  well  known  that  monolithic  microwave  circuits  are  lossy  compared  to 
hybrid  circuits.  Therefore  it  is  important  to  noise  figure  and  circuit 
performance  that  the  filters  be  affected  as  little  as  possible  by  loss. 

Figures  6  through  9  demonstrate  the  effects  of  uniform  element  Qs  of  25  on 
the  three  filter  types  discussed.  In  both  reflection  and  transmission,  the 
new  filter  is  least  affected  fay  loss.  The  new  filter  gives  .7  and  2.5  db  less 
loss  at  the  band  edges  and  .4  and  1.3  db  less  loss  at  the  band  center  than  the 
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mismatched  and  noraal  Chebyshev  filter a,  respectively.  Note  that 
pradiatortion  would  help  the  band  edge  response  but  hurt  the  band  center 
response.  The  end  result  is  that  the  new  filter  would  give  a  better  noise 
figure  than  either  of  the  Chebyshev  filters.  Table  IV  coapares  the  tolerancea 
of  the  lossy  filters.  Again,  the  new  filter  gives  an  overall  higher  tolerance 
to  both  eleaent  and  device  (load)  variations. 

Table  V  gives  the  tolerance  results  for  a  broadband  FET  feedback 
aaplifler.  This  aaplifier  used  a  750  p  wide  FET  of  Cf  -  .6pF  and  gg|  “  .065 
ahos  with  an  I£  matching  filter  (C|_,  L2)  and  a  feedback  resistor  (Rpg)  of  180 
ohas,  as  shown  in  Figure  10.  The  FET  model  included  source  inductance  and 
.  other  pertinent  components.  Two  different  matching  filters  were  used  with 
.'this  amplifier,  and  have  responses  shown  in  Figure  11.  The  first  design  was 
“tuned"  to  provide  zero  reflection  near  7  GHz ,  and  results  in  a  familiar 
humped  response.  The  second  response  was  designed  to  give  a  finite  but 
flatter  VSWR  across  the  entire  band.  'Oils  latter  design  is  less 
straightforward,  but  results  in  a  much  greater  tolerance  to  both  the  FET  and 
the  circuit  elements. 

Conclusion 

The  consideration  of  the  angle  between  the  Su(f)  vector  and  the  change 
in  Su(f)  due  to  load  reactance  variations  has  been  used  to  design  a  more 
tolerant  filter  characteristic.  This  angle  factor  also  explains  why  Chebyshev 
filters  have  a  low  tolerance  for  a  given  mismatch.  A  corrected  first  order 

approximation  to  load  reactance  tolerance  was  given  which  has  greater  accuracy 

r 

with  no  additional  circuit  calculations.  While  the  filter  was  designed 
specifically  for  a  greater  tolerance  to  load  reactance  variations,  the  filter 
showed  increased  tolerances  to  eleaent  variations  and  the  effects  of  loss. 
Since  Insensitivity  to  loss  and  the  greatest  possible  tolerances  are  essential 
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Co  designing  lower  cost  end  easily  tuned  microwave  circuits,  the  design 
techniques  described  above  should  aid  the  production  of  all  microwave 
circuits,  especially  monolithic  circuits. 

The  basic  concepts  of  this  paper  are  given  in  equations  (1)  and  (5). 

These  nay  be  generalized  to  state  that  the  sensitivity  of  any  element  in  a 
filter  is  proportional  to  the  product  of  the  amount  of  power  delivered  to  it 
from  the  source  and  the  amount  of  power  it  may  deliver  back  to  the  source;  and 
that  since  circuit  respones  are  vectors,  it  is  useful  to  view  the  problem  of 
maximizing  a  tolerance  as  that  of  maximizing  the  distance  the  response  must 
travel  before  it  violates  a  specification.  These  concepts  are  best  satisfied 
in  reciprocal  filtering  networks  by  using  a  mismatched  response  and  making  the 
angular  relationship  between  the  response  and  its  perturbation  due  to  any  one 
element  variation  equal  to  an  odd  multiple  of  90  degrees.. 
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TABLE  1 

FILTER  ELEMENT  VALUES 
f0  -  3.942  GHz 


1  Type 

BW 

(GHz) 

n 

Ci(pF) 

Lx(nH) 

GiCmQ) 

L2(nR)  C2(pF)  R2(Q) 

C3(pF)  L3(nH)  G3(mQ) 

|:heby shev 

.5 

2.5 

8.729 

.1867 

8.648 

18.11 

.09 

17.94 

8.729 

.1867 

8.64 

VChebyshev 
■(mismatched ) 

.5617 

2.228 

8.585 

.1899 

8.51 

13.99 

.1165 

13.86 

7.005 

.2327 

6.94 

New  filter 
|CF3.78 

1 mismatched)* 

.53 

1.978 

8.594 

.1897 

8.52 

12.49 

.131 

12.37 

5.138 

.3173 

5.09 

fButterworth** 
|( mismatched) 

2.0 

1.879 

.5735 

2.842 

N. 

.447 

3.64 

.213 

7.64 

I  *Thls  denotes 

a  third  order  filter 

with  reflection 

zeros 

at  .78 

of  the 

filter 

bandwidth. 

l**The  Butterworth  response  has  coo  large  a  bandwidth  and  is  used  only  for  illustrating  the 
importance  of  the  angle  factor. 


TABLE  II 

BANDPASS  MATCHING  FILTER  TOLERANCES 


Circuit  Parameters 


Allowed  Element  Tolerances  (I) 
for  Tmax  “  >333  and  3. 7-4.2  GHz  range 
(see  Figure  3) 


^Filter 

BH(GHz) 

Pmax 

Device 

Cf 

1 

2 

-Chebyshev 

.5 

.26 

i.i2f  ~ 

.642 

l2*  C2 

'  1.96*  CTT  Lt 

*Cheby  shev 

.5617 

.23 

1.432 

2.022 

L2>  C2 

3.392  Ci,  Li 

j'.New  Filter 

I 

.53 

.23 

1.82 

2.92 

Cl.  L2 

3.052  L2,  C2 

Butterworth* 

2. 

.278 

2.22 

|  Ideal  (eqn  10) 

3.92 

to  illustrate  how  the  angle  factor  allows  a  relatively  high  tolerance  even  when  the  reflection 
coefficient  is  large. 
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TABLE  III 

LOAD  REACTANCE  TOLERANCE  OF  LOSSLESS  FILTERS 


.  t 

Classical 
Tol  (X) 

(eqn  3) 
Freq  (GHz) 

Corrected 
Tol  (X) 

(eqn  4) 

Freq  (GHz) 

Actual 

Tol  (X)  Freq  (GHz) 

Chebyshev 

■mm 

■Km 

^753 

Chebyshev 

(mismatched) 

3.8 

n 

3.825 

1.43 

3.825 

New  Filter 
(mismatched) 

/ 

B 

3.94 

IB 

3.85 

1.83 

3.85 

TABLE  IV 

N. 

ELEMENT  TOLERANCES  OF  LOSSY  FILTERS 

element  Q  *  25 

Allowed  Element  Tolerances  (Z) 
for  Imax  "  *333  and  3. 7-4.2  GHz  range 
(see  Figure  3  and  Table  I) 


Circuit  Parameters 

Device 

Circuit 

Cf 

1 

2 

Filter 

BW(GHz)  PO.X 

. 

Chebyshev. 

.5 

78) 

~rsr 

CL.  L1 

l7H“ 

l2»  c2 

Chebyshev 

.5617 

.23 

1.6X 

•  9X 

Cl,  I* 

2.2X 

L2,  C2 

New  Filter 

.53 

.23 

1.8X 

.9X 

Cl,  Li 

2.4X 

L2,  c2 

< 


Figure  1.  A  Saleh  chart  containing  the  oaxlaua  tolerable  reflection,  raax, 
and  the  input  reflection  coefficient  at  fQ,  Su(f0).  The  dashed 
radial  line  represents  the  direction  of  miniaua  absolute  variation 
in  Sn(fQ).  The  alternating  line  perpendicular  to  the  Su(f0) 
vector  represents  the  direction  of  oaxlaua  absolute  tolerance  since 
the  distance  between  Sj_i(f0)  and  the  circle  Is  aaxiaized. 

Figure  2.  A  plot  of  VSWR.  versus  frequency  for  a  third  order  Chebyshev  filter 
(solid  line)  with  a  peak  reflection  of  .258,  a  mismatched  Chebyshev 
filter  (dashed  line)  with  maximum  and  minimum  reflection  of  .23  and 
.115  respectively,  and  a  large  bandwidth  mismatched  Butterworth 
filter  (alternating  line)  with  a  midband  reflection  coefficient  of 
.278. 

Figure  3.  The  third  order  lumped  filter  network  used  in  examining  the 

tolerances  of  filter  responses.  Inductors  and  capacitors  one 
through  three  represent  lossless  filter  elements.  &,  Cf  and  Lf 
represent  a  narrowband  model  of  an  FET  input  impedance.  The  ideal 
transformer  and  negative  elements  are  used  to  equalize  the  filter 
for  the  complex  FET  impedance  without  disturbing  the  filter 
components.  This  enables  filter  and  device  component  varlationa  to 
be  analyzed. 

Figure  4.  Input  reflection  plot  of  a  narrowband  mismatched  Butterworth  filter 
versus  frequency.  The  circles  represent  the  true  response,  and  the 
triangles  represent  the  response  with  Cf  designed  to  cover  3.7  to 
4.2  GHz,  but  only  3.75  to  4  GHz  is  shown  for  clarity.  The 
circle  is  for  a  reflection  coefficient  of  .333.  Note  how  the 
Butterworth  filter  fails  at  the  band  edge  and  has  ideal  tolerance 
characteristics  in  its  aidband. 

Figure  5.  Input  reflection  plot -of  the  Chebyshev  filter  in  Figure  2.  Only 

3.75  to  4  GHz  of  the  3.7  to  4.2  GHz  pasaband  is  shown  for  clarity. 
Again,  the  fma*  circle  is  at  a  reflection  of  .333,  the  circles 
represent  the  true  response,  and  the  triangles  represent  the 
response  for  a  2Z  increase  in  Cf.  Note  how  the  Chebyshev  response 
is  moved  in  a  nearly  radial  direction  at  the  response  peak.  This 
radial  variation  causes  the  Chebyshev  to  have  a  poor  tolerance  to 
load  variations. 

Figure  6.  VSWR  versus  frequency  plot  for  three  of  the  filters  described  in 
Tables  I  through  IV.  The  solid  line  is  the  Chebyshev  with  a  peak 
reflection  of  .258.  The  dashed  line  is  a  mismatched  Chebyshev  with 
reflection  maxima  and  minima  of  .23  and  .115  respectively.  The 
'alternating  line  is  the  new  filter  with  reflection  maxima  and 
mlnlna  of  .23  and  .115  respectively.  All  of  the  filters  have  the 
same  return-loss  bandwidth  product,  which  means  they  all  have  the 
same  ability  to  equalise  a  complex  load. 
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Figure  7.  VSUE  versus  frequency  plot  for  the  filters  described  In  Figure  6 

having  lossy  elements  with  Q's  of  25.  Note  the  narrow  bandwidth  of 
the  Chebyshev  (solid  line)  and  the  flat  response  of  the  new  filter 
(alternating  line). 

Figure  8.  Gain  versus  frequency  plot  for  the  lossless  filters  described  In 
Figure  6.  Note  the  Chebyshev  (solid  line),  mismatched  Chebyshev 
(dashed  line),  and  mismatched  F3.7g  (alternating  line)  responses. 
The  new  filter  has  a  double  peaked  response  for  all  orders. 

Figure  9.  Gain  versus  frequency  plot  for  the  filters  described  in  Figure  6 
having  lossy  elements  with  Q's  of  25.  The  Chebyshev  response 
(solid  line)  is  severely  rounded  and  has  the  greatest  loss.  The 
mismatched  Chebyshev  response  (dashed  line)  has  less  loss,  but  the 
lowest  loss  is  in  the  F3.7g  response  (alternating  line).  The  F3.7g 
response  is  relatively  flat  and  has  one  to  two  db  less  loss  than 
the  Chebyshev  response.  This  allows  less  tuning  and  a  lower  noise 
figure. 

Figure  10.  Model  of  an  FET  feedback  amplifier. 

Figure  11.  VSWR  versus  frequency  plot  for  the  Vet  feedback  amplifier  with  a 
tuned  input  matching  network  (solid  line),  and  a  flat  response 
matching  network  (dashed  line). 


IMPEDANCE  COORDINATES— 50-OHM  characteristic  impedance 
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A  Theory  of  Even  and  Odd  Order  Low-Pass  Impedance  Matching  Filters* 


A.  N.  Riddle  and  R.  J.  Trew 
Electrical  and  Computer  Engineering  Department 
North  Carolina  State  University 
Raleigh,  North  Carolina  27650 
(919)  737-2336 


ABSTRACT 

A  theory  of  impedance  matching  with  low-pass  filters  is  presented  which 
includes  both  even  and  odd  order  networks*  It  is  shown  that  the  odd  order 
networks  contain  an  insensitive  element  which  should  prove  useful  in 
minimizing  the  tuning  (and  thus  the  cost)  of  microwave  amplifiers*  A  fifth 
order  Chebyshev  network  is  used  as  an  example  in  which  the  last  shunt 
capacitor  in  the  matching,  network  is  varied  ±502  while  maintaining  a  VSWR  < 
1.5:1  over  nearly  an  octave* 


*This  work  was  supported  by  the  Army  Research  Office  under  contract 
No.  DAAG29-80-K-0080 . 
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A  Theory  of  Even  and  Odd  Order  Low-Pass  Impedance  Matching  Filters* 

A.  N.  Riddle  and  R.  J.  Trew 
Electrical  and  Computer  Engineering  Department 
North  Carolina  State  University 
Raleigh,  North  Carolina  27650 
(919)  737-2336 

INTRODUCTION 

Low-pass  LC  filter  networks  have  formed  the  basis  of  many  useful 
impedance  matching  networks  [1-5].  This  is  because  low-pass  networks  resemble 
the  topology  of  most  devices  (especially  bipolar  transistors),  low-pass 
networks  allow  bias  current  to  pass  to  the  device,  and  low-pass  networks  are 
less  sensitive  to  component  variations  than  band-pass  networks.  Various 
impedance-transforming  low-pass  networks  (matching  filters)  have  been  designed 
(2-4)  using  a  low-pass  filter  prototype  with  a  transformation  which  results  in 
an  even  order  low-pass  matching  filter  with  twice  the  number  of  components  as 
the  original  filter.  Tables  of  these  networks  have  proved  useful  in  amplifier 
design  [6].  It  should  be  emphasized  that  the  networks  described  in  [2-4]  are 
true  Impedance  matching  devices  which  are  capable  of  zero  reflection 
coefficients  and  unity  gains  at  specific  frequencies.  These  networks 
accomplish  this  by  virtue  of  the  frequency  transformation  used  in  their 
filtering  function.  The  filtering  function  is  typically  of  a  Butterworth, 
Chebyshev  or  Elliptic  function  form.  The  frequency  transformation  maps  u  to  a 
function  of  to*  which  creates  an  even  order  filter  with  an  even  number  of 
components.  The  resulting  filtering  function  has  twice  the  order  of  the 
filtering  function  before  application  of  the  frequency  transformation.  Vihile 
odd-order  matching  filters  [7]  are  at  least  as  useful  as  even-order  matching 
filters,  a  derivation  of  the  odd-order  filtering  function  has  not  been 
reported.  This  paper  presents  an  extension  to  the  theory  of  low-pass 
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impedance  matching  filters  by  including  a  derivation  of  the  filtering  function 
for  filters  with  even  and  odd  numbers  of  components. 

Since  definitive  treatments  of  broadband  equalization  networks  already 
exist  [8,9],  this  paper  will  not  discuss  those  methods.  Rather,  this  paper 
describes  the  derivation  of  the  impedance  matching  polynomials  to  be  used  in 
Youla's  method  of  equalization,  A  simple  method  of  reactance  equalization  is 
described  in  this  paper  to  illustrate  the  usefulness  of  odd-order,  low-pass 
matching  filters.  The  filters  described  in  this  paper  should  not  be  confused 
with  those  in  [8]  since  the  latter  achieved  impedance  matching  with  ideal 
transformers  and  this  paper  describes  low-pass  filters  capable  of  impedance 
transformation . 

A  main  point  in  this  paper  is  the  demonstration  of  how  the  symmetry  of 
odd-order  networks  allows  a  certain  component  to  vary  with  little  influence  on 
the  impedance  matching  properties  of  the  network.  The  insensitive  component 
is  the’  reactance  in  series  with  the  larger  terminating  resistance,  or  the 
reactance  in  parallel  with  the  smaller  terminating  resistance.  This  means 
that  the  shunt  capacity  of  a  bipolar  transistor  input  impedance  may  be  used  as 
part  of  the  matching  network,  and  yet  device-to-device  variations  will  cause 
negligible  change  in  circuit  performance.  This  property  results  because  the 
first  component  in  series  with  a  low  load  impedance  (with  respect  to  the 
generator  impedance)  has  a  greater  Influence  on  impedance  matching  than  the 
component  in  parallel  with  the  load  impedance.  Therefore,  such  circuits  have 
an  Inherent  insensitivity  to  device  variations.  This  inherent  insensitivity 
is  very  important  in  the  development  of  low-cost  microwave  amplifiers  using 
on-chip  matching  networks. 
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Traditionally,  impedance  matching  networks  are  divided  into  classes 
which  do  [1]  and  do  not  [2]  use  ideal  transformers.  This  paper  considers 
the  insertion- loss  theory  of  low-pass  impedance  matching  networks  without 
ideal  transformers.  Two  classes  of  polynomials  are  presented  in  order  to 
demonstrate  the  use  of  this  theory  in  both  narrow-  and  broad-band  matching 
networks. 

Throughout  this  paper,  network  frequency  and  impedance  scaling  nota¬ 
tion  [7]  is  used  so  that  Ri  *  1  fi  and  =  1  Hz. 

THEORY 

The  following  method  for  synthesizing  matching  networks  is  based  on 
what  is  loosely  termed  "insertion-loss"  theory.  This  theory  uses  the  ideas 
of  a  lossless  network  and  the  conservation  of  energy.  Consider  a  generator, 
lossless  network,  and  a  load  resistance  as  shown  in  Fig.  1.  It  is  assumed 
that  the  generator  is  capable  of  delivering  one  unit  of  power,  and  it  follows 
that  the  sum  of  the  power  reflected  from  the  network  and  the  power  transmitted 
to  the  load  must  equal  unity,  the  power  available  from  the  generator.  In  S- 
parame ter  notation [10] , 


|sh<p)I2  +  |s12(P)(2  -  l. 


(i) 


where  generally  p  -  o  +  ju  and  S^2(p)  *  S21(p)  (from  reciprocity).  Su(p) 
and  S21(p)  are  represented  by  the  polynomial  functions  D(p)  and  H(p)  as 
shown  in  equations  (2)  and  (3)j 


Sn<p) 


Usl 

H(p) 


(2) 


s22(p) 


1 

/e  H(p) 


(3) 
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The  polynomial  which  describes  Che  characteristic  function  of  the 
matching  filter,  D(p) ,  must  allow  the  transfer  function  of  the  network, 
S2i(p),  to  have  a  loss  at  p~0  which  is  dependent  upon  the  generator  and 
load  resistances,  and  zero  loss  at  some  specified  frequencies.  Using 
S-parameter  notation,  and  p-jui, 

S21<P)S21<-P)  -  - i i - ,  (4) 

eD(p)DC-p)  +  1  eH(p)H(-p) 

where,  with  reference  to  Fig.  1, 

CR!-R2)2 

e  “  ^R1R2  <5) 

and 

D(p)  -  (1  +  bp)M(p2),  (6) 

where  M(p2)  satisfies  certain  conditions,  and  b  is  a  parameter.  The 
conditions  on  M(p2)  are  that 

M(0)  -  1,  (7) 

and  for  some  given  u>i,  the  matching  frequencies, 

M(-u^)  -  0.  (8) 

The  parameter  b  is  a  real  number  which  is  approximately  proportional 
to  the  square  of  the  Q  of  the  load  resistance  and  the  component  nearest  the 
load.  An  empirical  expression  for  b  in  the  case  of  a  fifth  order  Chebyshev 
network,  which  is  given  in  the  synthesis  example,  is 

b  -  [ _ - \*  •  (9) 

\  .42(1  +  .^225.)  ! 

*2 


Note  that  for  b-0,  the  Cn  in  Fig.  1  becomes  zero  (i.e.,  an  open  circuit), 
the  network  is  of  even  order  and  the  expressions  reduce  to  those  described 
by  existing  theory  [8]. 


where 
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Then, 


where 


■  Che  upper  passbend  edge,  and 
u>L  "  the  lower  passband  edge. 


H(-w2) 


Cg(w) 

V»,> 


n  odd 
n  even  ’ 


(14) 


(15) 


and 


w-o 


SYNTHESIS 

It  will  prove  useful  to  give  an  example  of  matching  filter  synthesis. 
This  example  will  consist  of  a  fifth  order  Chebyshev  matching  filter.  A 
fifth  order  network  is  chosen  so  that  the  zeros  of  the  Chebyshev  function 
will  be  evident  and  the  insensitivity  of  the  element  nearest  the  load  can 
be  demonstrated.  Note  that  impedance  and  frequency  scaling  have  been  used. 

Consider  the  matching  filter  of  Fig.  2  terminated  with 

Rl  -  1.0  n 

R2  -  *3333  0  . 

Via  equation  (5)  it  is  found  that 


e  -  .3333 


154 


Let  us  take  <d0  *  1  and  choose  a  bandwidth  of  approximately  31Z. 

Therefore  it  is  chosen  that 

A  -  .3115, 

since  A  is  a  paremeter  roughly  equal  to  the  fractional  bandwidth.  Once 
A  and  uq  are  known,  both  and  may  be  found  by  using  equations  (12) 
and  (13),  so, 

uu  •  1.145 
UL  -  .8295  . 

Since  n  -  5  is  odd,  K  -  2  from  equation  (15),  and  so  the  second  order 
Chebyshev  polynomial  found  in  P-2  1  may  be  combined  with  equation  (11)  to  yield. 


M(-u2) 


C2Cw) 

c2<Wz> 


2W2-! 


a£_  2j/  -  4^2  +  -1)  • 

"  2«4  A2  A2  a2 

M(-  2)  -  — i -  (20.615  w*  -  41.229  <u2  +  19.615) 

19.615 


(16) 


(17) 

(18) 


and 


Then,  letting  b-1. 


M(p2)  -  1.051  p4  +  2.102  p2  +  1. 


D(p)  -  (l+p)M(p2) 

-  1.051  p5  +  1.051  p4  +  2.102  p3  +  2.102  p2  +  p  +  1.  (19) 

The  matching  filter  is  synthesized  from  this  characteristic  function 
by  using  the  relations, 

S21(P)S21<-P)  “  eD(p)D(_p)+1  cH(p)H(-p) 


1  -  8^^(p)S^^(”P) 
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These  relations  reduce  to 


Sii(p) 


DM 

H(p) 


(20) 


where  the  quadrantal  symmetry  of  the  zeros  of  eD(p)D(-p)+l  may  be  used  to 
efficiently  derive  H(p)  within  a  constant  multiplier.  This  constant  multi¬ 
plier  la  found  by  recognizing  that  since  these  are  low-pass  filters. 


so. 


lim 

p— 


Sn(p)  - 


lim 


H(p) 


lim  P(p) 
P"**  H(p) 

_  +  lim 
"  ~  P-** 


-  *  1. 


D(p) 


(21) 

(22) 


with  the  positive  sign  generally  used. 

The  network  elements  are  then  derived  by  a  continued  fraction 
expansion  of 


1  +  Su(p) 

Y(P>  "  -  . 

1  -  Sji(p)  (23) 

Substituting  equation  (19)  into  equations  (4),  (1),  and  (20)  gives  the 
result  of  equation  (23)  as 

y(p)  .  2.1  p5+3.87  p*+7.45  p3+8.58  p2+5.77  p+3  .  (24) 

1.76  p*+3.24  p3+4.38  p2  +  3.77  p+1 

A  continued  fraction  expansion  of  equation  (24)  gives  the  elements  of 
Fig.  2  as 

Cx  -  1.19  F 
L2  -  .79  H 
C3  -  2.94  F 
L4  -  .463  H 
end  C5  •  1.64  F» 
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It  should  be  recognized  that  if  C5  vas  the  parasitic  capacity  of  a 
load  impedance,  then  equation  (9)  could  have  been  used  to  calculate  b  so 
that  the  synthesized  network  would  achieve  a  C5  of  the  desired  value. 

In  the  previous  example. 


C 


C5r2  ^2  .  / (1.64) (. 3333)2  _  ^ 

.42(1+  ^1)  .42(1+  ) 

Rl  .3333' 


03 


which  is  so  close  to  the  chosen  b-1  as  to  give  an  excellent  approximation. 


SENSITIVITY 


In  order  to  demonstrate  the  insensitivity  of  the  matching  characteristics 
of  this  network  to  variations  in  C5,  where  may  be  taken  as  the  parasitic 
parallel  capacity  associated  with  the  load,  the  VSWR  of  the  network  was 
calculated  as  Cj  was  varied  *  50Z. 

The  solid,  dashed,  and  alternating  lines  in  Fig.  3  indicate  the  network's 
input  VSWR  for  C5  *  1.64  F,  Cj,«  1.5  (1.64)F  and  C5  *  .5  (1.64)F,  respectively. 
Note  that  the  variation  of  C5  shifts  the  center  frequency  of  the  network,  and 
increased  the  input  VSWR  (Fig.  3);  however,  the  network  obtains  a  VSWR  <1.5 
for  almost  an  octave  in  all  cases. 

The  insensitivity  of  these  networks  is  a  general  property  which  can  be 
deduced  from  the  Chain  (ABCD)  matrix  representation  [13] .  This  representation 
shows  that  the  use  of  different  resistors  in  a  doubly  terminated  network  will 
cause  the  transfer  function  to  be  less  dependent  on  the  last  element  at  one 
end  of  the  filter.  For  the  odd-order  filter  with  shunt  elements  at  the  input 
end  output,  the  element  nearest  the  low  resistance  termination  is  least  sensi¬ 
tive.  Thus,  the  odd-order  filter  is  the  only  network  which  is  capable  of  e 
desensitized  shunt  capacitor  on  the  low  lspedaace'  side  of  the  network. 
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This  desensitization  may  be  heuristically  explained  by  noting  that 
the  basic  matching  circuit  is  in  the  form  of  an  L-network.  A  match  between 
two  unequal  resistances  is  mainly  dependent  upon  the  first  element  shunting 
the  larger  resistance,  and  the  first  element  in  series  with  the  lower 
resistance.  Because  of  this,  elements  in  shunt  with  the  lower  resistance 
and  in  series  with  the  higher  resistance  play  a  lesser  role  in  the  matching 
circuit. 


Conclusions 

A  general  form  for  the  characteristic  function  of  a  low-pass  matching 
filter  has  been  introduced.  This  characteristic  function  allows  both  even 
and  odd-order  networks  to  be  derived.  Furthermore,  it  has  been  shown  that 
the  odd-order  networks  can  minimize  the  sensitivity  of  the  matching  filter 
to  variations  in  the  parallel  capacity  of  a  low  impedance  load.  This 
property  should  be  useful  for  monolithic  circuits  and  the  internal  switching 
networks  of  packaged  transistors.  In  the  Interest  of  generality,  two 
different  classes  of  polynomials  have  been  given  which  satisfy  the  require¬ 
ments  of  the  characteristic  function  for  narrow-band  and  broad-band 
applications. 

It  should  be  noted  that  the  parameter  b  is  roughly  proportional  to 
the  square  of  the  Q  of  the  last  element  and  the  termination  in  the  matching 
filter.  Therefore,  as  b  is  reduced,  the  bandwidth  of  the  filter  will 
Increase  slightly  as  the  odd  order  filter  degenerates  to  an  even  order 
filter. 

This  analysis  is  presented  to  further  the  understanding  of  impedance 
matching  networks.  Since  impedance  matching  is  sine  qua  none  in  micro¬ 
wave  amplifiers  end  antennas,  the  microwave  engineer  should  have  as  many 
tools  as  is  possible.- 
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Fifth  order  low-pass  matching  filter  with  element  values 
given  in  the  text. 
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